moredq calculates useful inertia tensor perturbation quantities
for an elastic moon with constant q.
ss  time (julian date)
cm  3rd principal moment of inertia of the Moon
didk,didq    partial of delta i wrt k2 or k2/q
diddk,diddq    corresponding partials of delta i dot (time deriv)
didkm,didqm    portion of above depending on moon mass
note that the calculations are explicitly dependent
only on the moon mass, and not at all on the earth,
except insofar as Mmoon=Mass(3)*Mass(10)
analytic expressions for elastic perturbations supplied by dh eckh

calculate sensitivity of delta I and delta I dot to change in the
rotation state, to be used in calculating the partial derivatives
of the angular acceleration
if iparm.gt.1 then also calculate sensitivity to orbit state

Calculate partials of elasticity/dissipation corrections wrt parameters
Information passed in common:
Mmc,Dmmc= Moon mass divided by C, and partials thereof
Mass= array of mass ratios
Output:
ddidp  = partial of inertia tensor correction
ddiddp = partial of inertia tensor rate correction