subr psrphs - j.f.chandler - 1984 jun
compute pulsar pulse phase observables
in this routine the delay from the arrival (free-space) at
the solar-system barycenter to the observation is converted to
pulse phase.