예제 #1
0
def test_on_mills():
    mills=959804082*1e3
    j2k_ott = marstime.j2000_offset_tt(marstime.julian_tt(marstime.julian(mills)))
    assert within_error(j2k_ott, 151.344, 0.001)
    ls = marstime.Mars_Ls(j2k_ott)
    cy = marstime.Clancy_Year(j2k_ott)
    my = marstime.Mars_Year(j2k_ott)
    assert within_error(ls, 0.0352, 0.001)
    assert within_error(cy, 24, 0.5)
    assert within_error(my, 25, 0.5)
예제 #2
0
def test_math_vs_numpy():
    try:
        import numpy as np
        use_numpy=True
    except:
        use_numpy=False
    
    if use_numpy:
        mil = 1073137591000
        jdut = marstime.julian(mil)
        marstime.use_numpy = False
        tt_math = marstime.utc_to_tt_offset(jdut)

        marstime.use_numpy = True
        tt_numpy = marstime.utc_to_tt_offset(jdut)
        assert within_error(tt_math,tt_numpy,1e-3)

        j2k = marstime.j2000_offset_tt(marstime.julian_tt(jdut))
        marstime.use_numpy = False
        mym = marstime.Mars_Year(j2k)
        marstime.use_numpy = True
        myn = marstime.Mars_Year(j2k)
        assert within_error(mym,myn,0.5)
#----        
        
        mil=0
        jdut = marstime.julian(mil)
        marstime.use_numpy = False
        tt_math = marstime.utc_to_tt_offset(jdut)

        marstime.use_numpy = True
        tt_numpy = marstime.utc_to_tt_offset(jdut)
        assert within_error(tt_math,tt_numpy,1e-4)

        j2k = marstime.j2000_offset_tt(marstime.julian_tt(jdut))
        marstime.use_numpy = False
        mym = marstime.Mars_Year(j2k)
        marstime.use_numpy = True
        myn = marstime.Mars_Year(j2k)
        assert within_error(mym,myn,0.5)
예제 #3
0
def test_j2000_offset_tt():
    assert marstime.j2000_offset_tt() > 2456203.-2451545.
    assert marstime.j2000_offset_tt(marstime.j2000_epoch())==0
    assert marstime.j2000_offset_tt(2451645)==100
예제 #4
0
def test_spirit_landing():
    mil = 1073137591000
    longitude=184.702
    latitude=-14.460
    jdut = marstime.julian(mil)
    assert within_error(jdut,2453008.07397,1e-3)
    
    tt_utc = marstime.utc_to_tt_offset(jdut)
    assert within_error(tt_utc,64.184,1e-3)
    
    jday_tt = marstime.julian_tt(jdut)
    assert within_error(jday_tt, 2453008.07471,1e-3)
    
    j2000_ott = marstime.j2000_offset_tt(jday_tt)
    assert within_error(j2000_ott, 1463.07471, 1e-3)

    m = marstime.Mars_Mean_Anomaly(j2000_ott)
    assert within_error(m, 66.06851,1e3)

    alpha = marstime.FMS_Angle(j2000_ott)
    assert within_error(alpha, 317.09363,1e-4)

    pbs = marstime.alpha_perturbs(j2000_ott)
    assert within_error(pbs, 0.01614, 1e-3)

    v_m = marstime.equation_of_center(j2000_ott)
    assert within_error(v_m, 10.22959, 1e-4)

    ls = marstime.Mars_Ls(j2000_ott)
    assert within_error(ls, 327.32322, 1e-4)

    eot = marstime.equation_of_time(j2000_ott)
    assert within_error(eot, -12.77557, 1e-4)

    mtc = marstime.Coordinated_Mars_Time(j2000_ott)
    assert within_error(mtc,13.16542, 1e-4)

    lmst = marstime.Local_Mean_Solar_Time(longitude,j2000_ott)
    assert within_error(lmst, 0.85196, 1e-4)

    ltst = marstime.Local_True_Solar_Time(longitude,j2000_ott)
    assert within_error(ltst, 0.00025, 1e-5)

    subsol = marstime.subsolar_longitude(j2000_ott)
    assert within_error(subsol, 4.70575,1e-4)

    dec = marstime.solar_declination(ls)
    assert within_error(dec, -13.42105,1e-2)

    rm = marstime.heliocentric_distance(j2000_ott)
    assert within_error(rm,1.47767, 1e-4)

    im = marstime.heliocentric_longitude(j2000_ott)
    assert within_error(im, 52.37469, 1e-4)

    bm = marstime.heliocentric_latitude(j2000_ott)
    assert within_error(bm, 0.08962, 1e-4)

#This is where we deviate significantly from the original marstime algorithm.
    sz = marstime.solar_zenith(longitude, latitude,j2000_ott)
    assert within_error(sz, 151.93895, 0.2)

    sa = marstime.solar_azimuth(longitude, latitude,j2000_ott)
    assert within_error(sa, 179.99225,1e-2)
예제 #5
0
def test_midnight_crossing():
    mil = 947116800000
    longitude=0.
    latitude=0.
    jdut = marstime.julian(mil)
    assert within_error(jdut,2451549.5,1e-3)
    
    tt_utc = marstime.utc_to_tt_offset(jdut)
    assert within_error(tt_utc,64.184,1e-3)
    
    jday_tt = marstime.julian_tt(jdut)
    assert within_error(jday_tt, 2451549.50074,1e-3)
    
    j2000_ott = marstime.j2000_offset_tt(jday_tt)
    assert within_error(j2000_ott, 4.50074, 1e-3)

    m = marstime.Mars_Mean_Anomaly(j2000_ott)
    assert within_error(m, 21.74548,1e3)

    alpha = marstime.FMS_Angle(j2000_ott)
    assert within_error(alpha, 272.74486,1e-4)

    pbs = marstime.alpha_perturbs(j2000_ott)
    assert within_error(pbs, 0.00142, 1e-3)

    v_m = marstime.equation_of_center(j2000_ott)
    assert within_error(v_m, 4.44191, 1e-4)

    ls = marstime.Mars_Ls(j2000_ott)
    assert within_error(ls, 277.18677, 1e-4)

    eot = marstime.equation_of_time(j2000_ott)
    assert within_error(eot, -5.18764, 1e-4)

    mtc = marstime.Coordinated_Mars_Time(j2000_ott)
    assert within_error(mtc,23.99431, 1e-4)

    lmst = marstime.Local_Mean_Solar_Time(longitude,j2000_ott)
    assert within_error(lmst, 23.99431, 1e-4)

    ltst = marstime.Local_True_Solar_Time(longitude,j2000_ott)
    assert within_error(ltst, 23.64847, 1e-4)

    subsol = marstime.subsolar_longitude(j2000_ott)
    assert within_error(subsol, 174.72703,1e-4)

    dec = marstime.solar_declination(ls)
    assert within_error(dec, -25.22838,1e-4)

    rm = marstime.heliocentric_distance(j2000_ott)
    assert within_error(rm, 1.39358, 1e-4)

    im = marstime.heliocentric_longitude(j2000_ott)
    assert within_error(im, 2.26270, 1e-4)

    bm = marstime.heliocentric_latitude(j2000_ott)
    assert within_error(bm, -1.35959, 1e-4)

    sz = marstime.solar_zenith(longitude, latitude,j2000_ott)
    assert within_error(sz, 154.26182, 1e-4)

    sa = marstime.solar_azimuth(longitude, latitude,j2000_ott)
    assert within_error(sa, 191.03687,1e-4)