see subroutine bodred of input link for explanation of
mcond(1-30) or mrcond(1-30)   called cond(1-6),con(1-24) in bodred
mcon1(1-12) or mrcon1(1-12)   called con1(1-12)          in bodred
epsm(1-6) or epsmr(1-6)       called eps(1-6)            in bodred
jdm1,jdmn0 (or jdmr0),jdm2    called jd1,jd0,jd2         in bodred
imn or ilib                   called itape               in bodred
km(1-100) or kmr(1-100)       called k(1-100)            in bodred
kim(1-n) or kir(1-n)          called ki(1-n)             in bodred
intm or intmr                 called int                 in bodred
mplnt                         called nplnt               in bodred
mcentr                        called ncentr              in bodred

ipar= number of partial derivatives on output tape + 1
iparm for orbit, iparmr for rotation
if ipar=1, no partials on tape
iparm and iparmr calculated in subroutine morset from
kim(1-n) and kir(1-n), respectively