setup initial condition for equations of motion and
equations for partial derivatives w.r.t. initial osculating
eliptic orbital elements
or initial position and velocity (depending on icnd)
*       start=1300

metric parameters beta or gamma not equal to one give
principle of equivalence violation (unless ratio of
gravitational to total energy is zero), as do non-zero
deltas or deltap

initialize central body tape reading
kkp(84)=0  ignore central body tape during integration
kkp(84)=1  read central body tape for positions, velocities only
kkp(84)=2  read central body tape for partials as well