reference compare

geospace-code · Oct 31, 2018 · 457dd39 · 457dd39
1 parent 25884b2
commit 457dd39
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/tests/test_kepler.py b/tests/test_kepler.py
@@ -5,6 +5,7 @@
 import xarray
 import georinex as gr
 from datetime import datetime, timedelta
+import numpy as np
 
 TGPS0 = datetime(1980, 1, 6)
 
@@ -28,6 +29,10 @@
     'M0': 2.24295542,
 }
 
+xref = -5.67841101e6
+yref = -2.49239629e7
+zref = 7.05651887e6
+
 time = TGPS0 + timedelta(weeks=910, seconds=4.03272930e5)
 
 dat = xarray.Dataset(sv, attrs={'svtype': 'G'},
@@ -36,3 +41,6 @@
 x, y, z = gr.keplerian2ecef(dat)
 
 print(x, y, z)
+
+magerr = np.sqrt((x-xref)**2 + (y-yref)**2 + (z-zref)**2)
+print('error magnitude [meters]', magerr)