for k=1,...,numspt we have
spot(k)  = 4 character observed spot name
nsplnt(k)= planet number of body on which spot lies (if.ge.0)
nsplnt(k)= -1,-2,-3 for various types of earth spot coordinates
(changed to 3 before return to calling program)
nsplnt(k)= -4 spot is a star or quasar
nsplnt(k)= -5 spot is a pulsar
spcord(1,k) =radius of spot relative to center of body (kilometers)
spcord(2,k) =east longitude of spot (degrees).  on moon or planet
longitude is angle in right hand coordinate system,
whereas on earth longitude is angle in left hand
coordinate system for observing sites.
for earth spots it is right handed coordinate system
spcord(3,k) =latitude of spot (degrees).
lspcrd(i,k) =0 spcord(i,j) not adjusted (i=1,3)
lspcrd(i,k) =1 spcord(i,j) adjusted in least squares analysis (i=1,3)
kspt        =flag indicating type of coordinates
if nsplnt=3 and kspt<0, same logic as for nsplnt<0
otherwise if kspt=-1, convert cylindrical to spherical
otherwise if kspt=-4, convert cartesian to spherical
otherwise if kspt=-5 spot is moving and its position
relative to planet must be supplied in the obslib,
in which case, spcord(1) is set to -5 on return

v. spot coordinates           (sptred)
for each spot
card 1a
columns
1- 4  spot name                                         (1a4)
5- 7  planet number                                     (i3)
8- 8  not used                                          (blank)
9-24  first  spot coordinate                            (f16.9)
or  9-26   (f18.11) ##
25-40  second spot coordinate                            (f16.9)
or 27-44   (f18.11) ##
41-56  third  spot coordinate                            (f16.9)
or 45-62   (f18.11) ##
57-58  not used in old format                            (blank)
64 velocity flag: if '6', another card follows  (1x,a1) ##
59-64  not used                                          (blank)
or 65-70 ##
65-70  lsp=1,0  adjust or not the spot coordinates       (3i2)
or 71-76   (3i2) ##
71-72  kspt  type of coordinates                         (i2)
or 77-78   (i2) ##

card 1b - implemented only in new format
note: velocities for pulsars are implemented separately, not here
1- 4  not used, but must be non-blank
5- 8  not used
9-26  upward spot velocity                            (f18.11) ##
27-44  westward spot velocity                          (f18.11) ##
45-62  northward spot velocity                         (f18.11) ##
63-70  reference epoch                                 (f8.0)   ##
71-76  ls= 1  adjust each of three spot velocities     (3i2)    ##
ls= 0  do not adjust
79-80  should be '##', but will be ignored

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

input:       l      = postion of spot data in the final array
(integer)
npl    = body identification number to select
(integer*4)