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()
