cappallo/ash   august 1970  subroutine shppar
modified june,1978 by r.b. goldstein and j. chandler
shppar calculates partial derivatives of height(lslope=1) or
slope(lslope=2) of a planetary surface wrt spherical harmonic
coefficients.
it is assumed that /legend/ contains leg. poly. for this co-lat.