m.ash   march 1974   subroutine efgout
write output efg satellite coordinates for cloudcroft 48 inch
optical observatory (+ others)
binary efg output for least squares fit of polynomials in efgfit
option to output angles instead of e,f,g
called from mdeldp. entry points efg1,efg2 called from radar.

jct(30) = 0 no e,f,g (+edot,fdot,gdot) tape output
earth centered rotating coordinates. e in true equator of
date towards greenwich meridian, g to north, f completes
right hand system
jct(30).ne.0 such tape output on data set iabs(jct(30))
unless jct(40) says there is angle output instead
.lt.0 no frintout of e,f,g
.gt.0 printout of e,f,g at end of dummy observation series
jct(31) = space defense center satellite number for e,f,g output
jct(32) = 0 no binary e,f,g output in subroutine efgout for fitting
polynomials in subroutine efgfit
jct(32).ne.0 such binary output on data set iabs(jct(32)) if
jct(30).ne.0
.gt.0 jct(30) output with jct(32) output
.lt.0 jct(30) output is supressed

jct(40) = 0 no angle output for optical observatories
jct(40).ne.0 such tape output on data set iabs(jct(30))
in format compatable with data general nova computer
jct(40).gt.0 printout of angles, data set jct(40) used as intermediate
ebcdic buffer
jct(40).le.0 no printout of angles
jct(41) = 0 angle output referred to true equinox and equator of date
for jct(40) output
jct(41).gt.0 angle output referred to mean equinox and equator of
jct(41). for eample, jct(41)=76 denotes output referred
to mean equinox and equator of 1976.0
jct(42) = 0 no refraction correction made in computing jct(40) angles
jct(42) = 1 refraction correction made in computing jct(42) angles
jct(43) =   file number for jct(40) angle output (0 is file 1, 1 is
file 2, etc, for data general nova control word at end of
514 character record)
jct(44) = 0 no punch card output for haystack
jct(44).gt.0 80 character card image output for haystack on data set
jct(44) (punch if 7) if jct(40).gt.0

jct(46) = 0 jct(40) print angle output has range in meters and
range rate in cm/min
jct(46) = 1 jct(40) print angle output has one way time delay in
seconds and one way doppler shift in hz

jct(47) = 0 usual jct(40) print and tape output
jct(47) > 0 jct(40) tape output has satellite radius instead
of doppler rate, print output has special les-8/9
format with 6 frequency doppler
jct(47) = 8 les-8 output
jct(47) = 9 les-9 output

jct(48) = 0 right ascension and hour angle rates are seconds
of time per second in jct(40) print output
jct(48) = 1 right ascension and hour angle rates are seconds
of arc per second in jct(40) print output

xplsc     = position and velocity of earth satellite, lunar
orbiter or spot on moon relative to center of earth
referred to mean equinox and equator of the
reference epoch at retarded
time (units are light seconds and light seconds per
second) (calculated in subroutine radmn)