def get_baryv(ra, dec, mjd, T, obs="GB"): """ get_baryv(ra, dec, mjd, T): Determine the average barycentric velocity towards 'ra', 'dec' during an observation from 'obs'. The RA and DEC are in the standard string format (i.e. 'hh:mm:ss.ssss' and 'dd:mm:ss.ssss'). 'T' is in sec and 'mjd' is (of course) in MJD. """ tts = pu.span(mjd, mjd + T / 86400.0, 100) nn = len(tts) bts = numpy.zeros(nn, dtype=numpy.float64) vel = numpy.zeros(nn, dtype=numpy.float64) presto.barycenter(tts, bts, vel, nn, ra, dec, obs, "DE200") return vel.mean()
def get_baryv(ra, dec, mjd, T, obs="GB"): """ get_baryv(ra, dec, mjd, T): Determine the average barycentric velocity towards 'ra', 'dec' during an observation from 'obs'. The RA and DEC are in the standard string format (i.e. 'hh:mm:ss.ssss' and 'dd:mm:ss.ssss'). 'T' is in sec and 'mjd' is (of course) in MJD. """ tts = pu.span(mjd, mjd+T/86400.0, 100) nn = len(tts) bts = numpy.zeros(nn, dtype=numpy.float64) vel = numpy.zeros(nn, dtype=numpy.float64) presto.barycenter(tts, bts, vel, nn, ra, dec, obs, "DE200") return vel.mean()