m.e.ash   aug 1970    subroutine stmrdr
main program for processing delay and doppler for signal
sent from a ground station to a satellite over to another
satellite and then down to another ground station

mocpar = 0 compar called in midst of least squares iteration
read all three data sets iobcon, iobs, iabs1 which become in
subroutines of compar iiobcn, iiobs, iiabs1
mocpar = 1 compar called at end of least squares iteration to
calculate dummy observations
read only iobcon which become in subroutines of compar iiobcn

nprec  = 0 nutation-precession determined midway between
receiving and sending times
nprec  = 1 nutation-precession determined at both receiving and
sending times

ndprec = 0 derivative of nutation-precession matrix not used in
determining velocities
ndprec = 1 derivative of nutation-precession matrix used in
determining velocities

neatm  =-1 effect of earth atmosphere on radar signal is
ignored
neatm  = 0 effect of earth atmosphere on time delay only is
included
neatm  = 1 effect of earth atmosphere on time delay and doppler
shift is included

savb(i,1), i=1,numsav,  saved precession-nutation,etc. from input
observation library tape for radar, optical
and first part of transit observation
savb(i,2) i=1,numsav   these saved quantities from iiobs

resltb(1,j) = time delay for radar observation (seconds)
right ascension for optical observation (sec.of time)
time that first limb of planet cuts first limb of sun
for transit observation (first part) (seconds)
resltb(2,j) = doppler shift for radar observation (cycles/sec)
declination for optical observatoon (sec.of arc)
time that second limb of planet cuts first limb of sun
for transit observation (first part) (seconds)

set up observation data set quantities
niobc,niobs,niabs1  indicate status of data sets  iiobc,
iiobs, iiabs1
-1  not to be read
0  ready to be read
1  already read
if  iiobc, iiobs, or iiabs1  are initially zero,  data set
is not read in any case.