def ut2jde(y, m, d):
    ''' convert UT to JDE, compensate the delta T '''

    jd = g2jd(y, m, d)
    deltat = deltaT(y, m) / 86400.0
    return jd + deltat
 convert UT to JDE, compensate the delta T 
aa.py