def lightabbr_high(jd):
    '''compute light abberation based on A & A p156
    the error will be less than 0.001"

    '''
    t = (jd - J2000) / 365250.0

    # variation of the Sun's longitude
    var_lon = 3548.330
    for r in LIGHTABBR_TABLE:
        var_lon += t ** r[0] * r[1] * sin(r[2] + r[3] * t)

    t = (jd - J2000) / 36525.0
    t2 = t * t
    t3 = t * t2

    # mean anomaly of the Sun
    M = (357.52910 + 35999.0503 * t - 0.0001559 * t2
                                    - 0.00000048 * t3) * DEG2RAD
    # the eccentricity of Earth's orbit
    e = 0.016708617 - 0.000042037 * t - 0.0000001236 * t2
    # Sun's equation of center
    C = ((1.9146 - 0.004817 * t - 0.000014 * t2) * sin(M)
         + (0.019993 - 0.000101 * t) * sin(2 * M)
         + 0.00029 * sin(3 * M)) * DEG2RAD
    # true anomaly
    v = M + C
    # Sun's distance from the Earth, in AU
    R = (1.000001018 * (1 - e * e)) / (1 + e * cos(v))

    res = -0.005775518 * R * var_lon * ASEC2RAD

    return res
compute light abberation based on A & A p156
    the error will be less than 0.001"

    
aa.py