r. goldstein     april, 1976
interpolations are now done by calls to mntrp,emtrp,
pltrp, and sbtrp.
polynomial extrapolations are now done by dguess

oct. 1976..large sections of code replaced by calls to
subroutines propco, radrel. and phadop
june 1978.. converted to implicit real*8 , interface with plcms
and streamlined radrel

xslcns(1-6,j) = position and velocity (not used) of the center
of mass of the solar system relative to the sun
at receive time (j=1), reflection time (j=2)
and send time (j=3)

correct for propagation effects by calling propco
set up ipct as follows:
-2: first leg of phase delay doppler
-1: second leg of phase delay doppler
1: time delay

there are three possible doppler observables:
lopler=-1  observable is phase
lopler= 1  observable is instantaneous frequency
loppler=0 has been removed from the code and has been gone
since before 1975.