def test_solar_zenith_and_elevation():
j2day=151.2737 #ls=0
x = marstime.subsolar_longitude(j2day)
assert within_error(x,114.113,1e-3)
assert within_error(marstime.solar_zenith(x,0,j2day),0,1e-4)
assert within_error(marstime.solar_zenith(x,15,j2day),15,1e-4)
assert within_error(marstime.solar_zenith(x+15,0,j2day),15,1e-4)
assert within_error(marstime.solar_elevation(x,45,j2day)-
marstime.solar_zenith(x,45,j2day)
,0,1e-4)
j2day=349.8778 #ls=90
x = marstime.subsolar_longitude(j2day)
assert within_error(x,232.7006,1e-3)
assert within_error(marstime.solar_zenith(x,0,j2day),25.441,1e-3)
assert within_error(marstime.solar_elevation(x,0,j2day),64.5581,1e-4)
try:
a= marstime.solar_zenith(0,100,0)
assert False
except ValueError:
assert True
except:
assert False
try:
a= marstime.solar_zenith(0,-100,0)
assert False
except ValueError:
assert True
except:
assert False
import marstime
import datetime
import scipy
import scipy.optimize
import argparse
import test_simple_julian_offset
if __name__ == "__main__":
parser = argparse.ArgumentParser(description="Testsuite for marstime library")
parser.add_argument("date",type=str, nargs="?",
help="Earth date in ISO format YYYY/MM/DD")
parser.add_argument("time",type=str, nargs="?",
help="Earth time in ISO format HH:ii:ss")
parser.add_argument("longitude",type=float,
help="East Longitude")
parser.add_argument("latitude",type=float,
help="East Longitude")
args = parser.parse_args()
dt = args.date + " " + args.time
jdate = test_simple_julian_offset.simple_julian_offset(datetime.datetime.strptime(dt, "%Y/%m/%d %H:%M:%S"))
west_longitude = marstime.east_to_west(args.longitude)
print marstime.solar_zenith(west_longitude, args.latitude, jdate)
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)
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)
def test_solar_zenith():
j2day=151.2737 #ls=0
x = marstime.subsolar_longitude(j2day)
assert within_error(x,114.113,1e-3)
assert within_error(marstime.solar_zenith(x,0,j2day),0,1e-4)
