in0  - fortran unit for editing input stream
nstop- cumulative error count
jdpad- adjustment, if any, for all jd1-jd2 integration ranges
meaning of kall:
-1 - just initialize
0 - initialize and read in all nmlst2 namelists
1 - already initialized, just read in one namelist

jdem0   = initial time for earth-moon barycenter integration
jdmn0   = initial time for moon integration
jder0   = initial time for earth rotation integration
jdmr0   = initial time for moon rotation integration
jdpl0(j)= initial time for planet nplnt(j) integration, j=1,numpln
initial time is julian day number, with numerical
integration commencing from midnight beginning of day of
calendar date corresponding to julian day number. if 0,
no integration. if initial conditions are adjusted, there
will be integration on second and subsequent least squares
analysis iterations because initial time will be set to
initial time on tape during comparison of theory and
observation. if negative, checkpoint restart mode.
see explanation of jdbdy0 above or explanation of jd0 in
subroutine bodred where these quantiies are read.

Aplnt(j)=8 character planet name, j=1,u_mxpl=16
Aplnt(-3)=' EMBARY '
Aplnt(-2)='  MOON  '
Aplnt(-1)=' EROTAT '
Aplnt(0) =' MROTAT '
Aplnt(17)='  MOON  '
Aplnt(18)='  SUN   '
earth, moon, sun names are setup in block data routine for main link.

nplnt(j)= planet number, j=1,numpln
1 mercury
2 venus
4 mars
5 jupiter
6 saturn
7 uranus
8 neptune
9 pluto
11-30  other natural planets,asteroids,satellites
31,... artificial space probes
-1,-2,-4,...,-9,-11,... rotation for planet iabs(nplnt)
3,10,-3,-10 do not exist in nplnt(j),j=1,numpln, but instead
are built into nplnt(j),j=-3,0,  imagined
to be associated with earth-moon barycenter, moon,
earth rotation, moon rotation
0 does not exist in nplnt(j),j=1,numpln, imagined to be
associated with sun

npcent(j)=central body for planet nplnt(j), j=1,numpln
-1 sun with integration in satellite-probe link
0  sun with integration in planet link if nplnt.le.30
3  earth (not earth-moon barycenter)
10 moon
other numbers as above

planet shape for nshp gt 0
special coding needed to transfer shape coefficients
and their l-vectors from overlay arrays in
/scoef/ into /scoef4/

nshp=2, local shape model-altitude grid
must transfer real*4 coefs from /scoef/ to
real*10 slots in /scoef4/ using array transf(_,4,___)
to maintain allignment of different planet models
to allow one planet to use r*4 grid model, while another
planet used r*8 fourier or sph. harmonic model