Skip to content

Reading.py module mis-calculates ordinates from absolutes and baselines

The method get_ordinates() in /geomagio/residual/Reading.py mis-calculates e_ord and h_ord. Here is what they should be:

e_ord = h_abs * np.sin(np.radians(d_ord))
h_ord = np.sqrt(h_abs**2 - e_ord**2) - h_bas

This small change will break several unit tests in adjusted_tests.py, so update these as well.

Edited by Erin (Josh) Rigler