def test_apxq2g_lowprecision():
fa.loadapxsh(os.path.join(os.path.dirname(apexpy.__file__), 'apexsh.dat'),
2000)
glat, glon, error = fa.apxq2g(60, 15, 100, -1)
assert_allclose(glat, 51.00891876220703)
assert_allclose(glon, -66.11973571777344)
assert_allclose(error, -9999.0)
def test_g2q2d():
fa.loadapxsh(os.path.join(os.path.dirname(apexpy.__file__), 'apexsh.dat'),
2000)
for lat in [0, 30, 60, 89]:
for lon in [-179, -90, 0, 90, 179]:
qlat, qlon, _, _, _ = fa.apxg2q(lat, lon, 100, 0)
glat, glon, _ = fa.apxq2g(qlat, qlon, 100, 1e-10)
assert_allclose(glat, lat, atol=0.01)
assert_allclose(glon, lon, atol=0.01)
def test_apxq2g():
""" Test fortran quasi-dipole to geographic
"""
fa.loadapxsh(os.path.join(os.path.dirname(apexpy.__file__), 'apexsh.dat'),
2000)
glat, glon, error = fa.apxq2g(60, 15, 100, 1e-2)
assert_allclose(
[glat, glon, error],
[50.97946548461914, -66.16902923583984, 0.00010020843910751864],
atol=1e-6)
def test_apxg2q():
"""Test fortran apex geographic to quasi-dipole
"""
fa.loadapxsh(os.path.join(os.path.dirname(apexpy.__file__), 'apexsh.dat'),
2000)
qlat, qlon, f1, f2, f = fa.apxg2q(60, 15, 100, 1)
assert_allclose(qlat, 56.531288146972656)
assert_allclose(qlon, 94.1068344116211)
assert_allclose(f1, [1.07978308, 0.10027108], rtol=1e-6)
assert_allclose(f2, [-0.24546318, 0.90718889], rtol=1e-6)
assert_allclose(f, 1.0041800737380981)
def test_set_epoch():
A = Apex(date=2000.2, refh=300)
assert_allclose(A.year, 2000.2)
ret_2000_2_py = A._geo2apex(60, 15, 100)
A.set_epoch(2000.8)
assert_allclose(A.year, 2000.8)
ret_2000_8_py = A._geo2apex(60, 15, 100)
assert ret_2000_2_py != ret_2000_8_py
fa.loadapxsh(A.datafile, 2000.2)
ret_2000_2_apex = fa.apxg2all(60, 15, 100, 300, 0)[2:4]
fa.loadapxsh(A.datafile, 2000.8)
ret_2000_8_apex = fa.apxg2all(60, 15, 100, 300, 0)[2:4]
assert ret_2000_2_apex != ret_2000_8_apex
assert_allclose(ret_2000_2_py, ret_2000_2_apex)
assert_allclose(ret_2000_8_py, ret_2000_8_apex)
def test_apxg2all():
fa.loadapxsh(os.path.join(os.path.dirname(apexpy.__file__), 'apexsh.dat'),
2000)
qlat, qlon, mlat, mlon, f1, f2, f, d1, d2, d3, d, e1, e2, e3 = fa.apxg2all(
60, 15, 100, 300, 1)
assert_allclose(qlat, 56.531288146972656)
assert_allclose(qlon, 94.1068344116211)
assert_allclose(mlat, 55.94841766357422)
assert_allclose(mlon, 94.1068344116211)
assert_allclose(f1, [1.07978308, 0.10027108], rtol=1e-6)
assert_allclose(f2, [-0.24546318, 0.90718889], rtol=1e-6)
assert_allclose(f, 1.0041800737380981)
assert_allclose(d1, [0.94957006, 0.2569305, 0.09049489], rtol=1e-6)
assert_allclose(d2, [0.10011058, -1.07805467, -0.33892459], rtol=1e-6)
assert_allclose(d3, [0.00865366, 0.27327025, -0.86666465], rtol=1e-6)
assert_allclose(d, 1.1003910303115845)
assert_allclose(e1, [1.02692986, 0.08382936, 0.03668636], rtol=1e-6)
assert_allclose(e2, [0.24740215, -0.82374191, -0.25726584], rtol=1e-6)
assert_allclose(e3, [0.01047837, 0.33089212, -1.04940987], rtol=1e-6)
def test_set_epoch():
apex_out = Apex(date=2000.2, refh=300)
assert_allclose(apex_out.year, 2000.2)
ret_2000_2_py = apex_out._geo2apex(60, 15, 100)
apex_out.set_epoch(2000.8)
assert_allclose(apex_out.year, 2000.8)
ret_2000_8_py = apex_out._geo2apex(60, 15, 100)
assert ret_2000_2_py != ret_2000_8_py
fa.loadapxsh(apex_out.datafile, 2000.2)
ret_2000_2_apex = fa.apxg2all(60, 15, 100, 300, 0)[2:4]
fa.loadapxsh(apex_out.datafile, 2000.8)
ret_2000_8_apex = fa.apxg2all(60, 15, 100, 300, 0)[2:4]
assert ret_2000_2_apex != ret_2000_8_apex
assert_allclose(ret_2000_2_py, ret_2000_2_apex)
assert_allclose(ret_2000_8_py, ret_2000_8_apex)
import os
import pytest
from numpy.testing import assert_allclose
import apexpy
from apexpy import fortranapex as fa
##############################################################################
# NOTE: the test results (numbers) were obtained by running the code that is #
# tested, therefore the tests below only check that nothing changes when #
# refactoring etc., and not if the results are actually correct #
##############################################################################
fa.loadapxsh(os.path.join(os.path.dirname(apexpy.__file__), 'apexsh.dat'), 2000)
def test_apxg2q():
qlat, qlon, f1, f2, f = fa.apxg2q(60, 15, 100, 1)
assert_allclose(qlat, 56.531288146972656)
assert_allclose(qlon, 94.1068344116211)
assert_allclose(f1, [1.07978308, 0.10027108], rtol=1e-6)
assert_allclose(f2, [-0.24546318, 0.90718889], rtol=1e-6)
assert_allclose(f, 1.0041800737380981)
def test_apxg2all():
qlat, qlon, mlat, mlon, f1, f2, f, d1, d2, d3, d, e1, e2, e3 = fa.apxg2all(60, 15, 100, 300, 1)
assert_allclose(qlat, 56.531288146972656)
import os
import pytest
from numpy.testing import assert_allclose
import apexpy
from apexpy import fortranapex as fa
##############################################################################
# NOTE: the test results (numbers) were obtained by running the code that is #
# tested, therefore the tests below only check that nothing changes when #
# refactoring etc., and not if the results are actually correct #
##############################################################################
fa.loadapxsh(os.path.join(os.path.dirname(apexpy.__file__), 'apexsh.dat'),
2000)
def test_apxg2q():
qlat, qlon, f1, f2, f = fa.apxg2q(60, 15, 100, 1)
assert_allclose(qlat, 56.531288146972656)
assert_allclose(qlon, 94.1068344116211)
assert_allclose(f1, [1.07978308, 0.10027108], rtol=1e-6)
assert_allclose(f2, [-0.24546318, 0.90718889], rtol=1e-6)
assert_allclose(f, 1.0041800737380981)
def test_apxg2all():
qlat, qlon, mlat, mlon, f1, f2, f, d1, d2, d3, d, e1, e2, e3 = fa.apxg2all(
60, 15, 100, 300, 1)
