(a purist would say 'segment' instead of 'link'. however,
this program on the ibm 370 computer, now 120000 cards long,
grew out of a 6000 card long program on the ibm 7094
computer. the overlay terminology in the 7094 operating
system used 'link', whereas the terminology in os/360 is
'segment'. the word 'link' is imbedded at a number of
places in the program in printout formats.)

gauss = gaussian gravitational constant
pi    = ratio between the circumference and diameter of a
circle
twopi = 2*pi
pitwo = pi/2
convd = pi/180
convds= convd/3600
convhs= twopi/86400

common       routine(s)             stand-alone routine(s)
------       ---------------------  ----------------------
anctrl       anset
bdctrl       prmred
bdydta       prtcrd
bernum       prdobs
dtparm       prmred
empcnd       bdyred                 addmoon
eqnphs       cmpar2
ethhar       bdyred
france       bdyred
ktrap        prdobs
lcntrl       prmred
leon         comset                 libang
lfix         partl
ltrap        cmpar1, comset         plastest
mascon       prmred
mnsprt       cmpar2                 plastest
monhar       bdyred
mtrap        cmpar1
namtim                              addmoon, intrpp, orbdif, plastest
pmoon
number       cmpar1, comset         plastest
plndta                              pmoon
plnhar       bdyred
psrstf       bdyred
scoef4       bdyred
skymap       cmpar2
sitcrd       cmpar1                 plastest
smlbdy       bdyred
stats        comset                 plastest
tidal        fermtr
tidal        radar
wrkcom       bodred
xprcom       morset, sbset

21    fn      plantlnk  perform earth-moon barycenter and planet-
asteroid numerical integrations (motion and
partial derivatives), write results on
magnetic tape or disk (one data set per
body).

26    sbfn    satpblnk  perform satellite-probe numerical
integrations (motion and partial
derivatives), write results on magnetic tape
or disk (one data set per body).

40    compar  extralnk  process observations of sun, moon, planets,
asteroids, satellites, artificial space
probes, stars
(compute observed minus theoretical values
and partial derivatives of observations,
write results on magnetic tape or disk)

77    nrmfrm  analzlnk  restore saved solution of normal equations
or form normal equations from subset of
saved normal equations. solve by inverting
matrix with double precision arithmetic
using fact that matrix is symmetric.
(extended precision intermediate
operations if desired)