Added a fix for the timeobj nan bug
@@ -227,6 +227,7 @@ def bkjd_to_time(bkjd, timecorr, timslice, bjdref=2454833.):
bjd = bkjd + bjdref
jd = bjd - timecorr
jd += (0.25 + 0.62 * (5 - timslice)) / 86400.
+ jd[~np.isfinite(jd)] = 0
return Time(jd, format='jd')
