kscrd(k)    =-1 cylindrical observing site coordinates relative to
center of mass of earth
kscrd(k)    = 0 spherical observing site coordinates relative to
center of mass of earth
kscrd(k)    = n.gt.0, spherical observing site coordinates relative
to center of mass of body with planet number n
(not programmed yet, program can only handle earth based
observing sites for now)
all coordinates are in body fixed reference systems

for spherical site coordinates
scord(1,k)  = radius in kilometers
scord(2,k)  = west longitude in degrees
(west of greenwich for geocentric)
scord(3,k)  = north latitude in degrees

for cylindrical site coordinates
scord(1,k)  = radius in x-y plane in kilometers
scord(2,k)  = west longitude in degrees
(west of greenwich for geocentric)
scord(3,k)  = z distance in kilometers

for cartesian site coordinates (for conversion to cylindrical)
x direction is the 0,0 of latitude and longitude
scord(1,k)  = x distance in kilometers
scord(2,k)  = y distance in kilometers
scord(3,k)  = z distance in kilometers

for all sites
scord(4,k)  = upward velocity (mm/yr)
scord(5,k)  = westward velocity (mm/yr)
scord(6,k)  = northward velocity (mm/yr)
t0site(k)   = reference epoch for site coordinates: julian day number

iv observing sites             (sitred)
for each site
card 1a
columns
1- 8  site name (only first 4 char used operationally)  (1a8)
9-24  first  site coordinate                            (f16.9)
or  9-26   (f18.11) ##
25-40  second site coordinate                            (f16.9)
or 27-44   (f18.11) ##
41-56  third  site coordinate                            (f16.9)
or 45-62   (f18.11) ##
58  velocity flag: if '6', then another card follows  (1x,a1)
or    64   (1x,a1) ##
59-64  not used
or 65-70 ##
65-70  ls= 1  adjust each of three site coordinates      (3i2)
ls= 0  do not adjust each of three site coordinates
or 71-76   (3i2) ##
71-72  ks  type of input site coordinates                (i2)
or 77-78   (i2) ##
ks=  n  (n.gt.0)  spherical coords. on body  n
ks=  0            spherical coords. on earth
ks= -1            cylindrical coords. on earth
ks= -2            geodetic (spheroidal) coords. on earth
with radius and flattening read in
or setup before sitred called
ks= -3            geodetic (spheroidal) coords. on earth
same as ks=-2 except
radius and flattening individually input
ks= -4            XYZ coords. on earth (convert to sph)
also applies to velocities, if any, which
will then be converted to topocentric
if the second site card is omitted, then the velocity is assumed to
be zero, and the components are not to be adjusted
## if columns 79-80 are '##", use the alternate (new) format

card 1b
1- 8  not used, but must be non-blank
9-24  upward site velocity                              (f16.9)
or  9-26   (f18.11) ##
25-40  westward site velocity                            (f16.9)
or 27-44   (f18.11) ##
41-56  northward site velocity                           (f16.9)
or 45-62   (f18.11) ##
57-64  reference epoch                                   (f8.0)
or 63-70   (f8.0) ##
65-70  ls= 1  adjust each of three site coordinates      (3i2)
ls= 0  do not adjust each of three site coordinates
or 71-76   (3i2) ##

ks =positive spherical coordinates on body ks
ks =  0   geocentric spherical coordinates on earth
ks = -1   cylindrical coordinates on earth
ks = -2   geodetic (spheroidal) coordinates on earth with
radius and flattening as input to program as whole
ks = -3 same as ks=-2 except that radius and flattening
to be specially read in