|
@@ 129-146 (lines=18) @@
|
| 126 |
|
|
| 127 |
|
del code, trace_lat, trace_lon, trace_r |
| 128 |
|
|
| 129 |
|
def test_convert_A2G_TRACE(self): |
| 130 |
|
"""Test convert from magnetic to geodetic coordinates using trace""" |
| 131 |
|
code = aacgmv2._aacgmv2.A2G + aacgmv2._aacgmv2.TRACE |
| 132 |
|
trace_lat = [30.7644, 50.3958] |
| 133 |
|
trace_lon = [-94.1809, -77.8019] |
| 134 |
|
trace_r = [1133.6277, 305.7156] |
| 135 |
|
|
| 136 |
|
for i,dargs in enumerate(self.date_args): |
| 137 |
|
aacgmv2._aacgmv2.set_datetime(*dargs) |
| 138 |
|
(self.mlat, self.mlon, |
| 139 |
|
self.rshell) = aacgmv2._aacgmv2.convert(self.lat_in[i], |
| 140 |
|
self.lon_in[i], |
| 141 |
|
self.alt_in[i], code) |
| 142 |
|
np.testing.assert_almost_equal(self.mlat, trace_lat[i], decimal=4) |
| 143 |
|
np.testing.assert_almost_equal(self.mlon, trace_lon[i], decimal=4) |
| 144 |
|
np.testing.assert_almost_equal(self.rshell, trace_r[i], decimal=4) |
| 145 |
|
|
| 146 |
|
del code, trace_lat, trace_lon, trace_r |
| 147 |
|
|
| 148 |
|
def test_convert_high_denied(self): |
| 149 |
|
"""Test for failure when converting to high altitude geodetic to |
|
@@ 110-127 (lines=18) @@
|
| 107 |
|
|
| 108 |
|
del lat_comp, lon_comp, r_comp |
| 109 |
|
|
| 110 |
|
def test_convert_G2A_TRACE(self): |
| 111 |
|
"""Test convert from geodetic to magnetic coordinates using trace""" |
| 112 |
|
code = aacgmv2._aacgmv2.G2A + aacgmv2._aacgmv2.TRACE |
| 113 |
|
trace_lat = [48.1948, 58.1633] |
| 114 |
|
trace_lon = [57.7588, 81.0756] |
| 115 |
|
trace_r = [1.1775, 1.0457] |
| 116 |
|
|
| 117 |
|
for i,dargs in enumerate(self.date_args): |
| 118 |
|
aacgmv2._aacgmv2.set_datetime(*dargs) |
| 119 |
|
(self.mlat, self.mlon, |
| 120 |
|
self.rshell) = aacgmv2._aacgmv2.convert(self.lat_in[i], |
| 121 |
|
self.lon_in[i], |
| 122 |
|
self.alt_in[i], code) |
| 123 |
|
np.testing.assert_almost_equal(self.mlat, trace_lat[i], decimal=4) |
| 124 |
|
np.testing.assert_almost_equal(self.mlon, trace_lon[i], decimal=4) |
| 125 |
|
np.testing.assert_almost_equal(self.rshell, trace_r[i], decimal=4) |
| 126 |
|
|
| 127 |
|
del code, trace_lat, trace_lon, trace_r |
| 128 |
|
|
| 129 |
|
def test_convert_A2G_TRACE(self): |
| 130 |
|
"""Test convert from magnetic to geodetic coordinates using trace""" |
aburrell /
aacgmv2
