Return an astropy Time value as MJD in longdouble
The returned value is accurate to within a nanosecond, while the precision of long double MJDs (near the present) is roughly 0.7 ns.