def test_0(self):
dat = astrodate.AstroDate()
dat.set_longitude(self.lng0)
dat.set_zone_correction(self.zc0)
dat.set_daylight_savings(self.dst0)
dat.set_with_tuple(self.lct0)
d = dat.to_lct(self.zc0_1, self.dst0_1)
t = dat.get_tuple()
self.assertTupleEqual(t, self.lct0)
t = d.get_tuple()
self.assertTupleEqual(t, self.lct0_1)
d = dat.to_lct(self.zc0_2, self.dst0_2)
t = dat.get_tuple()
self.assertTupleEqual(t, self.lct0)
t = d.get_tuple()
self.assertTupleEqual(t, self.lct0_2)
dat.to_utc()
t = dat.get_tuple()
self.assertTupleEqual(t, self.utc0)
dat.to_lct()
t = dat.get_tuple()
self.assertTupleEqual(t, self.lct0)
def test_to_date_tuple_from_julian(self):
dat = astrodate.AstroDate()
dat.set_with_julian(0.0)
d = dat.get_tuple()
self.assertTupleEqual(d, (-4712, 1, 1, 12, 0, 0, 'utc'))
dat.set_with_julian(2299160.5)
d = dat.get_tuple()
self.assertTupleEqual(d, (1582, 10, 15, 0, 0, 0, 'utc'))
dat.set_with_julian(2448251.3125)
d = dat.get_tuple()
self.assertTupleEqual(d, (1990, 12, 25, 19, 30, 0, 'utc'))
dat.set_with_julian(2448988.5)
d = dat.get_tuple()
self.assertTupleEqual(d, (1993, 1, 1, 0, 0, 0, 'utc'))
dat.set_with_julian(2448989.0)
d = dat.get_tuple()
self.assertTupleEqual(d, (1993, 1, 1, 12, 0, 0, 'utc'))
dat.set_with_julian(2449078.5)
d = dat.get_tuple()
self.assertTupleEqual(d, (1993, 4, 1, 0, 0, 0, 'utc'))
dat.set_with_julian(2451545.0)
d = dat.get_tuple()
self.assertTupleEqual(d, (2000, 1, 1, 12, 0, 0, 'utc'))
def test_ref2(self):
dat = astrodate.AstroDate()
dat.set_longitude(self.lng_ref2)
dat.set_with_tuple(self.utc_ref2)
dat.to_lst()
t = dat.get_tuple()
self.assertTupleEqual(t, self.lst_ref2)
def calculate_mean_obliquity(self, epochTD=None):
if epochTD is not None:
d = astrodate.AstroDate().alloc_with_epochTD(epochTD)
else:
if self.date is None:
raise ValueError, "Date (TD) is required!"
self.date.to_tdt()
d = self.date
n = earth.Nutation()
n.calculate_with_julianTD(d.get_julian())
return n.get_mean_obliquity()
def test_equatorial_to_ecliptic(self):
d = astrodate.AstroDate().alloc_with_julian(astrodate.J1950, "tdt")
c = AstroCoord()
c.set_date(d)
c.set_with_tuple((9, 34, 53.6, 19, 32, 14.2, u'equ|hrs|ra'))
# print(c.get_pretty_coordinate())
c.to_ecliptic()
# print(c.get_pretty_coordinate())
t = c.get_tuple()
self.assertTupleEqual(t, (139, 41, 9.29, 4, 52, 22.00, u'ecl'))
def test_equatorial_to_galactic(self):
d = astrodate.AstroDate().alloc_with_julian(astrodate.J1950, "tdt")
c = AstroCoord()
c.set_date(d)
c.set_with_tuple((10, 21, 0.0, 10, 3, 11, u'equ|hrs|ra'))
# print(c.get_pretty_coordinate())
c.to_galactic()
# print(c.get_pretty_coordinate())
t = c.get_tuple()
self.assertTupleEqual(t, (232, 14, 52.51, 51, 7, 20.16, u'gal'))
def test_ecliptic_to_equatorial(self):
d = astrodate.AstroDate().alloc_with_julian(astrodate.J1950, "tdt")
c = AstroCoord()
c.set_date(d)
c.set_with_tuple((139, 41, 10, 4, 52, 31, u'ecl'))
# print(c.get_pretty_coordinate())
c.to_equatorial(1980.0)
# print(c.get_pretty_coordinate())
t = c.get_tuple()
self.assertTupleEqual(t, (9, 34, 53.58, 19, 32, 14.18, u'equ|hrs|ra'))
def test_horizon_to_equatorial(self):
c = AstroCoord()
c.set_latitude(52.0)
c.set_with_tuple((19, 20, 3.64, 283, 16, 15.7, u'hor'))
# print(c.get_pretty_coordinate())
c.to_equatorial()
# print(c.get_pretty_coordinate())
t = c.get_tuple()
self.assertTupleEqual(t, (5, 51, 44.00, 23, 13, 10.0, 'equ|hrs|ha'))
# print
d = astrodate.AstroDate()
d.set_with_tuple((2017, 1, 1, 0, 0, 0.0, u'lst'))
c.set_latitude(31.9583)
c.set_date(d)
c.set_with_tuple((77, 21, 40.56, 39, 33, 35.22, u'hor'))
# print(c.get_pretty_coordinate())
c.to_equatorial()
# print(c.get_pretty_coordinate())
t = c.get_tuple()
self.assertTupleEqual(t, (23, 17, 15.7, 41, 16, 9.02, 'equ|hrs|ha'))
c.to_right_ascension()
# print(c.get_pretty_coordinate())
c.to_degrees()
# print(c.get_pretty_coordinate())
t = c.get_tuple()
self.assertTupleEqual(t, (10, 41, 4.50, 41, 16, 9.02, 'equ|deg|ra'))
# print
d.set_longitude(-1.9166667)
d.set_with_tuple((1998, 8, 10, 23, 10, 00, u'utc'))
d.to_lst()
c.set_latitude(52.5)
c.set_date(d)
c.set_with_tuple((49, 10, 7.93, 269, 8, 48.01, u'hor'))
# print(c.get_pretty_coordinate())
c.to_equatorial()
# print(c.get_pretty_coordinate())
t = c.get_tuple()
self.assertTupleEqual(t, (3, 37, 31.93, 36, 28, 0.0, 'equ|hrs|ha'))
c.to_right_ascension()
# print(c.get_pretty_coordinate())
t = c.get_tuple()
self.assertTupleEqual(t, (16, 41, 42.0, 36, 28, 0.0, 'equ|hrs|ra'))
def test_nutation(self):
d = astrodate.AstroDate().alloc_with_tuple(
(1987, 4, 10, 0, 0, 0.0, u'td'))
jde = d.get_julian()
n = earth.Nutation()
n.calculate_with_julianTD(jde)
nutation_in_longitude = n.get_nutation_in_longitude()
nutation_in_obliquity = n.get_nutation_in_obliquity()
mean_obliquity = n.get_mean_obliquity()
true_obliquity = n.get_true_obliquity()
self.assertEqual(nutation_in_longitude, -0.001052332550403307)
self.assertEqual(nutation_in_obliquity, 0.002622947155040129)
self.assertEqual(mean_obliquity, 23.440946290957324)
self.assertEqual(true_obliquity, 23.443569238112364)
def test_equatorial_mode(self):
d = astrodate.AstroDate()
d.set_longitude(-64.0)
d.set_with_tuple((1980, 4, 22, 14, 36, 51.67, u'utc'))
d.to_lst()
c = AstroCoord()
c.set_date(d)
c.set_with_tuple((18, 32, 21.0, -8, 12, 5.9, u'equ|hrs|ra'))
# print(Rigel_Equ_Hrs_RA)
c.to_hour_angle()
# print(c.get_pretty_coordinate())
t = c.get_tuple()
self.assertTupleEqual(t, (5, 51, 44.23, -8, 12, 5.9, u'equ|hrs|ha'))
c.to_degrees()
# print(c.get_pretty_coordinate())
t = c.get_tuple()
self.assertTupleEqual(t, (87, 56, 3.45, -8, 12, 5.9, u'equ|deg|ha'))
def test_to_julian_from_date_tuple(self):
dat = astrodate.AstroDate()
dat.set_with_tuple((-4712, 1, 1, 12, 0, 0))
jd = dat.get_julian()
self.assertEquals(jd, 0.0)
dat.set_with_tuple((1582, 10, 15))
jd = dat.get_julian()
self.assertEquals(jd, 2299160.5)
dat.set_with_tuple((1980, 2, 17))
jd = dat.get_julian()
self.assertEquals(jd, 2444286.5)
dat.set_with_tuple((1985, 2, 17))
jd = dat.get_julian()
self.assertEquals(jd, 2446113.5)
dat.set_with_tuple((1990, 12, 25, 19, 30, 0))
jd = dat.get_julian()
self.assertEquals(jd, 2448251.3125)
dat.set_with_tuple((1993, 1, 1))
jd = dat.get_julian()
self.assertEquals(jd, 2448988.5)
dat.set_with_tuple((1993, 1, 1, 12, 0, 0))
jd = dat.get_julian()
self.assertEquals(jd, 2448989.0)
dat.set_with_tuple((1993, 4, 1))
jd = dat.get_julian()
self.assertEquals(jd, 2449078.5)
dat.set_with_tuple((2000, 1, 1, 12, 0, 0))
jd = dat.get_julian()
self.assertEquals(jd, 2451545.0)
def test_crn(self):
d = astrodate.AstroDate().alloc_with_tuple((1975, 1, 27, 0, 0, 0.0, u'lct'))
crn = sun.calculate_crn_with_dateTD(d)
self.assertEqual(crn, 1624)
def test_2(self):
dat = astrodate.AstroDate()
dat.set_longitude(self.lng2)
dat.set_zone_correction(self.zc2)
dat.set_daylight_savings(self.dst2)
dat.set_with_tuple(self.lct2)
dat.to_lct()
t = dat.get_tuple()
self.assertTupleEqual(t, self.lct2)
dat.to_utc()
t = dat.get_tuple()
self.assertTupleEqual(t, self.utc2)
dat.to_lct()
t = dat.get_tuple()
self.assertTupleEqual(t, self.lct2)
dat.set_with_tuple(self.utc2)
dat.to_gst()
t = dat.get_tuple()
self.assertTupleEqual(t, self.gst2)
dat.to_utc()
t = dat.get_tuple()
self.assertTupleEqual(t, self.utc2)
dat.set_with_tuple(self.gst2)
dat.to_lst()
t = dat.get_tuple()
self.assertTupleEqual(t, self.lst2)
dat.to_gst()
t = dat.get_tuple()
self.assertTupleEqual(t, self.gst2)
dat.set_with_tuple(self.utc2)
dat.to_tdt()
t = dat.get_tuple()
self.assertTupleEqual(t, self.tdt2)
dat.to_utc()
t = dat.get_tuple()
self.assertTupleEqual(t, self.utc2)
dat.set_with_tuple(self.lct2)
dat.to_lst()
t = dat.get_tuple()
self.assertTupleEqual(t, self.lst2)
dat.to_lct()
t = dat.get_tuple()
self.assertTupleEqual(t, self.lct2)
dat.set_with_tuple(self.lct2)
dat.to_tdt()
t = dat.get_tuple()
self.assertTupleEqual(t, self.tdt2)
dat.to_lct()
t = dat.get_tuple()
self.assertTupleEqual(t, self.lct2)
def test_1(self):
dat = astrodate.AstroDate()
dat.set_longitude(self.lng1)
dat.set_zone_correction(self.zc1)
dat.set_daylight_savings(self.dst1)
dat.set_with_tuple(self.lct1)
d = dat.to_lct(self.zc1_1, self.dst1_1)
t = dat.get_tuple()
self.assertTupleEqual(t, self.lct1)
t = d.get_tuple()
self.assertTupleEqual(t, self.lct1_1)
d = dat.to_lct(self.zc1_2, self.dst1_2)
t = dat.get_tuple()
self.assertTupleEqual(t, self.lct1)
t = d.get_tuple()
self.assertTupleEqual(t, self.lct1_2)
dat.to_utc()
t = dat.get_tuple()
self.assertTupleEqual(t, self.utc1)
dat.to_lct()
t = dat.get_tuple()
self.assertTupleEqual(t, self.lct1)
dat.set_with_tuple(self.utc1)
dat.to_gst()
t = dat.get_tuple()
self.assertTupleEqual(t, self.gst1)
dat.to_utc()
t = dat.get_tuple()
self.assertTupleEqual(t, self.utc1)
dat.set_with_tuple(self.gst1)
dat.to_lst()
t = dat.get_tuple()
self.assertTupleEqual(t, self.lst1)
dat.to_gst()
t = dat.get_tuple()
self.assertTupleEqual(t, self.gst1)
dat.set_with_tuple(self.utc1)
dat.to_tdt()
t = dat.get_tuple()
self.assertTupleEqual(t, self.tdt1)
dat.to_utc()
t = dat.get_tuple()
self.assertTupleEqual(t, self.utc1)
dat.set_with_tuple(self.lct1)
dat.to_lst()
t = dat.get_tuple()
self.assertTupleEqual(t, self.lst1)
dat.to_lct()
t = dat.get_tuple()
self.assertTupleEqual(t, self.lct1)
dat.set_with_tuple(self.lct1)
dat.to_tdt()
t = dat.get_tuple()
self.assertTupleEqual(t, self.tdt1)
dat.to_lct()
t = dat.get_tuple()
self.assertTupleEqual(t, self.lct1)
if __name__ == "__main__":
jde = 2448908.5
p = calculate_position_with_julianTD(jde)
print("Latitude0: {}".format(p.get_latitude0()))
print("Latitude: {}".format(p.get_latitude()))
print("Apparent Longitude: {}".format(p.get_apparent_longitude()))
print("True Longitude0: {}".format(p.get_true_longitude0()))
print("True Longitude: {}".format(p.get_true_longitude()))
print("Radius: {}".format(p.get_earth_sun_radius()))
print
year = 2017
ad = astrodate.AstroDate().alloc_now_utc()
jdve = p.calculate_season_start(year, VERNAL_EQUINOX)
ad.set_with_julian(jdve, astrodate.TIME_MODE_TDT)
ve = ad.get_pretty_string()
print("Vernal Equinox: {} ({})".format(ve, jdve))
jdss = p.calculate_season_start(year, SUMMER_SOLSTICE)
ad.set_with_julian(jdss, astrodate.TIME_MODE_TDT)
ss = ad.get_pretty_string()
print("Summer Solstice: {} ({})".format(ss, jdss))
jdae = p.calculate_season_start(year, AUTUMNAL_EQUINOX)
ad.set_with_julian(jdae, astrodate.TIME_MODE_TDT)
ae = ad.get_pretty_string()
print("Autumnal Equinox: {} ({})".format(ae, jdae))
