m.e. ash  feb 1968  subroutine optic
ash/becker     oct 1967    subroutine optic
computation of theoretical value of right ascension and
declination for a meridian (transit) circle or photographic
observation of the sun, moon or a planet

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

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.