'Turn a time to MJD, returning integer and fractional parts.
Parameters
jd1 : float, array, or Time
first part of two-part JD, or Time object
jd2 : float or array, optional
second part of two-part JD.
Default is 0., ignored if jd1 is `~astropy.time.Time`.
Returns
mjd : float or array
integer part of MJD
utc : float or array
fractional part of MJD'

def mjd_utc(self, jd1, jd2=0.0):

    try:
        (jd1, jd2) = (jd1.utc.jd1, jd1.utc.jd2)
    except Exception:
        pass
    mjd = np.floor(((jd1 - MJD_ZERO) + jd2))
    utc = ((jd1 - (MJD_ZERO + mjd)) + jd2)
    return (mjd, utc)
