r.king and r.cappallo   july 1977   subroutine morhar
evaluation of the effect of the gravity fields of the earth and
moon in the right side of the equations for the lunar orbit and
rotation.  called by morfn, this routine incorporates (at present)
m.slade's coding from monfn for the orbital effect and king and
cappallo's coding from mrtfn for the rotation effect. the two sets
of code should be standardized soon and the effects of 4th degree
forces and torques added.

ncall= -2 setup once per iteration of a given step for partial
= -1 setup once per iteration of a given step for motion
=  0 setup once per step
=  k evaluate acceleration for equation k.gt.0
j    = pointer into integration output array

ecor  = earth relative to sun (re is distance)
mcor  = moon relative to sun (rm is distance)
mecor = moon relative to earth (rem is distance)
mcor etc. are equatorial coordinates
smcor etc. are selenodetic coordinates

at this point, JH is the order of the cosine tesseral term,
II is the index counting from C21, and
DCNM is the partial of the specified coefficient w.r.t. the
parameter of interest, i.e., either 1 (for the coefficent itself)
or the appropriate partial (for beta, gamma, or J2 indirect terms)