mored calculates and applys corrections to the lunar rotational
equations of motion to model elasticity and dissipation. entry
moredp calculates the rhs of the variational eqs. wrt:
k2(moon) = mrcond(12), and
k2 * lunar elastic response lag = mrcond(13).
see rjc's thesis for relevant formulae.

call subroutines to calculate perturbations to inertia tensor
and partials wrt k2 and a dissipation parameter.
MOREDT does calculation for kt (constant-time-lag, kmr83=0 or 1)
MOREDQ does calculation for k/q (constant-q model, kmr83=2)
if partial derivatives are needed, these subroutines are called
again to return the partials wrt the rotational state