See https://github.com/wbalmer/backtrack/issues/16#issuecomment-1864149058 for the analysis. Because the conversion is so easy and only has to be done once while reading the observations, we should convert the JD_UTC's of our datapoints to JD_TT so the targets are evaluated at the right time.
if jd_utc defines the epochs of our observed datapoints the needed times are
astropy.time.Time(jd_utc,scale="utc",format="jd").tt.jd
If needed all the jd_tt's can be converted back to jd_utc for plotting / printing purposes. Not sure if the conversions is fast enough to keep everything in UTC until we call app_star. For the default plots you will not see the ~68 seconds difference.
See https://github.com/wbalmer/backtrack/issues/16#issuecomment-1864149058 for the analysis. Because the conversion is so easy and only has to be done once while reading the observations, we should convert the JD_UTC's of our datapoints to JD_TT so the targets are evaluated at the right time. if jd_utc defines the epochs of our observed datapoints the needed times are
astropy.time.Time(jd_utc,scale="utc",format="jd").tt.jd
If needed all the jd_tt's can be converted back to jd_utc for plotting / printing purposes. Not sure if the conversions is fast enough to keep everything in UTC until we call app_star. For the default plots you will not see the ~68 seconds difference.