|
@@ 108-125 (lines=18) @@
|
| 105 |
|
|
| 106 |
|
del code, trace_lat, trace_lon, trace_r |
| 107 |
|
|
| 108 |
|
def test_convert_A2G_TRACE(self): |
| 109 |
|
"""Test convert from magnetic to geodetic coordinates using trace""" |
| 110 |
|
code = aacgmv2._aacgmv2.A2G + aacgmv2._aacgmv2.TRACE |
| 111 |
|
trace_lat = [30.7644, 50.3958] |
| 112 |
|
trace_lon = [-94.1809, -77.8019] |
| 113 |
|
trace_r = [1133.6277, 305.7156] |
| 114 |
|
|
| 115 |
|
for i,dargs in enumerate(self.date_args): |
| 116 |
|
aacgmv2._aacgmv2.set_datetime(*dargs) |
| 117 |
|
(self.mlat, self.mlon, |
| 118 |
|
self.rshell) = aacgmv2._aacgmv2.convert(self.lat_in[i], |
| 119 |
|
self.lon_in[i], |
| 120 |
|
self.alt_in[i], code) |
| 121 |
|
np.testing.assert_almost_equal(self.mlat, trace_lat[i], decimal=4) |
| 122 |
|
np.testing.assert_almost_equal(self.mlon, trace_lon[i], decimal=4) |
| 123 |
|
np.testing.assert_almost_equal(self.rshell, trace_r[i], decimal=4) |
| 124 |
|
|
| 125 |
|
del code, trace_lat, trace_lon, trace_r |
| 126 |
|
|
| 127 |
|
def test_convert_high_denied(self): |
| 128 |
|
"""Test for failure when converting to high altitude geodetic to |
|
@@ 89-106 (lines=18) @@
|
| 86 |
|
|
| 87 |
|
del lat_comp, lon_comp, r_comp |
| 88 |
|
|
| 89 |
|
def test_convert_G2A_TRACE(self): |
| 90 |
|
"""Test convert from geodetic to magnetic coordinates using trace""" |
| 91 |
|
code = aacgmv2._aacgmv2.G2A + aacgmv2._aacgmv2.TRACE |
| 92 |
|
trace_lat = [48.1948, 58.1633] |
| 93 |
|
trace_lon = [57.7588, 81.0756] |
| 94 |
|
trace_r = [1.1775, 1.0457] |
| 95 |
|
|
| 96 |
|
for i,dargs in enumerate(self.date_args): |
| 97 |
|
aacgmv2._aacgmv2.set_datetime(*dargs) |
| 98 |
|
(self.mlat, self.mlon, |
| 99 |
|
self.rshell) = aacgmv2._aacgmv2.convert(self.lat_in[i], |
| 100 |
|
self.lon_in[i], |
| 101 |
|
self.alt_in[i], code) |
| 102 |
|
np.testing.assert_almost_equal(self.mlat, trace_lat[i], decimal=4) |
| 103 |
|
np.testing.assert_almost_equal(self.mlon, trace_lon[i], decimal=4) |
| 104 |
|
np.testing.assert_almost_equal(self.rshell, trace_r[i], decimal=4) |
| 105 |
|
|
| 106 |
|
del code, trace_lat, trace_lon, trace_r |
| 107 |
|
|
| 108 |
|
def test_convert_A2G_TRACE(self): |
| 109 |
|
"""Test convert from magnetic to geodetic coordinates using trace""" |
aburrell /
aacgmv2
