m.e.ash   june 1969    subroutine hpartl
calculate partial derivative w.r.t. gravitational potential
harmonic coefficients for artificial space probe observations
or observations of the moon

kicki - calling routine: 1=radar, 2=optic, 3=trnsit, 4=interf
length- size of L vector
lhar  - input L vector
mhar  - L vector from old observation library
nqlnti- object whose integration is to be searched for partials
ntype - code for type of harmonics: 31=zonal, 41=tesseral cosine,
51=tesseral sine, 46+56 not implemented
nstrti- starting index into ki vector
mti   - flag indicating presence of old obslib if >0