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

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.

special sun-moon mode
xsun(1-3)=position of satellite relative to sun with aberration
corrections
xsun(7)=r**2
xsun(8)=r
xsun(9)=pitch angle or azimuth positive from yaw to roll
xsun(10)=roll angle or elevation.
these angles in ideal satellite reference frame pitch to north,
roll along satellite velocity vector
xsun(11)=fraction of moon's disk illuminated by sun
as seen from satellite
xsun(12)=moon phase angle as seen from satellite
xsun(13)=moon angular diameter as seen from satellite