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

ntmdly =-1 classical determination of time delay
ntmdly = 0 special relativistic determination of time delay
(same as classical because deviation of earth
velocity from its average is small)
ntmdly = 1 general relativistic determination of time delay

ndop   =-1 classical determination of doppler shift
ndop   = 0 special relativistic determination of doppler shift
(same as classical because deviation of earth
velocity from its average is small)
ndop   = 1 general relativistic determination of doppler shift

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.