r. b. goldstein  september 1977
input:  st in days from con1(1)
output: a,b,inc,psi and trig functions of inc,psi
algorithim:    calls rotab to calculate a,b at beginning and
end of a "pass". all values of all returned
parameters within the pass are obtained from
linear interpolation.