Right now model.get_barycentric_toas() returns an astropy Quantity with units of days.
I think it might be preferable to have it return a Time() object with scale="tdb".
For one thing, I tried to use format_toa_line() to write out those barycentric times but it caused an exception. I think very few places use this function so far, so I think it should not be too big a deal to change. Also, this would keep the timescale along with the MJD value.
Right now
model.get_barycentric_toas()
returns an astropy Quantity with units of days. I think it might be preferable to have it return aTime()
object with scale="tdb". For one thing, I tried to useformat_toa_line()
to write out those barycentric times but it caused an exception. I think very few places use this function so far, so I think it should not be too big a deal to change. Also, this would keep the timescale along with the MJD value.