ash/forni  november 1966    subroutine plnhed
printout of data page for planet numerical integration
write first two records of output ephemeris tape
updated feb. 1980   kcl: ptcon added

see subroutine bodred of input link for explanation of
cond(1-6)       called cond(1-6)   in bodred
con(1-24)       called con(1-24)   in bodred
con1(1-12)      called con1(1-12)  in bodred
tcon(1-30)      called tcon(1-30)  in bodred
epsp(1-6)       called eps(1-6)    in bodred
name            called name        in bodred
jdp1,jdp0,jdp2  called jd1,jd0,jd2 in bodred
jplnt           called itape       in bodred
kp(1-100)       called k(1-100)    in bodred
kkp(1-100)      called kk(1-100)   in bodred
ki(1-n)         called ki(1-n)     in bodred
intp            called int         in bodred
intp1           called int1        in bodred
intp2           called int2        in bodred
nplnt           called nplnt       in bodred
ncentr          called ncentr      in bodred
ihrp,iminp,secp called ihr,imin,secin bodred
icnd            called icnd        in bodred

iparp= number of partial derivatives on output tape + 1
if iparp=1, no partials on tape
iparp calculated in subroutine setup from values of kp(1-30)
the meaning of the partials j=2,iparp is given by kp(1-30).

klan = 0  for earth-moon barycenter integration
klan = 1,...u_mxpl  for planet integration,indicating that record klan+4
of iplcon data set and arrays with second dimension
klan from /namtim/ and /empcnd/ were used in putting
data into /petuna/