test nav3 header add

Author: scivision

README.md

@@ -71,16 +71,16 @@ python -m pytest
 
 ## Usage
 
-The simplest command-line use is through the top-level `ReadRinex` script.
+The simplest command-line use is through the top-level `python -m georinex.read` script.
 Normally you'd use the `-p` option with single files to plot, if not converting.
 
 * Read single RINEX3 or RINEX 2 Obs or Nav file:
   ```sh
-  ReadRinex myrinex.XXx
+  python -m georinex.read myrinex.XXx
   ```
 * Read NetCDF converted RINEX data:
   ```sh
-  ReadRinex myrinex.nc
+  python -m georinex.read myrinex.nc
   ```
 * Batch convert RINEX to NetCDF4 / HDF5 (this example for RINEX 2 OBS):
   ```sh
@@ -215,7 +215,7 @@ obs = gr.load('myfile.o', use='E')
 ```
 
 loads only Galileo data by the parameter E.
-`ReadRinex` allow this to be specified as the -use command line parameter.
+`python -m georinex.read` allow this to be specified as the -use command line parameter.
 
 If however you want to do this after loading all the data anyway, you can make a Boolean indexer
 

src/georinex/nav3.py

@@ -7,8 +7,6 @@
 import math
 from datetime import datetime
 
-
-#
 from .rio import opener, rinexinfo
 from .common import rinex_string_to_float
 

src/georinex/tests/test_nav3.py

@@ -5,10 +5,16 @@
 from datetime import datetime
 import georinex as gr
 
-#
+
 R = Path(__file__).parent / "data"
 
 
+def test_nav3header():
+    hdr = gr.rinexheader(R / "demo.17n")
+    assert hdr['IONOSPHERIC CORR']['GPSA'] == approx([1.1176e-08, -1.4901e-08, -5.9605e-08, 1.1921e-07])
+    assert hdr['TIME SYSTEM CORR']['GPUT'] == approx([-3.7252902985e-09, -1.065814104e-14, 61440, 1976])
+
+
 def test_time():
     times = gr.gettime(R / "VILL00ESP_R_20181700000_01D_MN.rnx.gz")
 
@@ -193,10 +199,10 @@ def test_large_all(sv, size):
 )
 def test_ref(rfn, ncfn):
     """
-    python ReadRinex.py tests/data/galileo3.15n -o r3galileo.nc
-    python ReadRinex.py tests/data/demo.17n -o r3gps.nc
-    python ReadRinex.py tests/data/qzss3.14n -o r3qzss.nc
-    python ReadRinex.py tests/data/demo3.10n -o r3sbas.nc
+    python -m georinex.read tests/data/galileo3.15n -o r3galileo.nc
+    python -m georinex.read tests/data/demo.17n -o r3gps.nc
+    python -m georinex.read tests/data/qzss3.14n -o r3qzss.nc
+    python -m georinex.read tests/data/demo3.10n -o r3sbas.nc
     """
     pytest.importorskip("netCDF4")
 

src/georinex/tests/test_obs2.py

@@ -193,8 +193,8 @@ def test_one_sv():
 @pytest.mark.parametrize("use", (None, {"G", "R", "S"}))
 def test_all_systems(tmp_path, use):
     """
-    ./ReadRinex.py tests/demo.10o -o r2all.nc
-    ./ReadRinex.py tests/demo.10n -o r2all.nc
+    python -m georinex.read tests/demo.10o -o r2all.nc
+    python -m georinex.read tests/demo.10n -o r2all.nc
     """
     pytest.importorskip("netCDF4")
 
@@ -222,7 +222,7 @@ def test_all_systems(tmp_path, use):
 
 @pytest.mark.parametrize("use", ("G", ["G"]))
 def test_one_system(use):
-    """./ReadRinex.py tests/demo.10o -u G -o r2G.nc"""
+    """python -m georinex.read tests/demo.10o -u G -o r2G.nc"""
     pytest.importorskip("netCDF4")
 
     truth = xarray.open_dataset(R / "r2G.nc", group="OBS")
@@ -233,7 +233,7 @@ def test_one_system(use):
 
 
 def test_multi_system():
-    """./ReadRinex.py tests/demo.10o -u G R -o r2GR.nc"""
+    """python -m georinex.read tests/demo.10o -u G R -o r2GR.nc"""
     pytest.importorskip("netCDF4")
 
     truth = xarray.open_dataset(R / "r2GR.nc", group="OBS")
@@ -245,7 +245,7 @@ def test_multi_system():
 
 def test_all_indicators():
     """
-    ./ReadRinex.py tests/demo.10o -useindicators  -o r2all_indicators.nc
+    python -m georinex.read tests/demo.10o -useindicators  -o r2all_indicators.nc
     """
     pytest.importorskip("netCDF4")
 
@@ -258,7 +258,7 @@ def test_all_indicators():
 
 def test_meas_indicators():
     """
-    ./ReadRinex.py tests/demo.10o -useindicators -m C1 -o r2_C1_indicators.nc
+    python -m georinex.read tests/demo.10o -useindicators -m C1 -o r2_C1_indicators.nc
     """
     pytest.importorskip("netCDF4")
 

src/georinex/tests/test_obs3.py

@@ -141,7 +141,7 @@ def test_bad_system():
 @pytest.mark.parametrize("use", ("G", ["G"]))
 def test_one_system(use):
     """
-    ./ReadRinex.py -q tests/demo3.10o  -u G -o r3G.nc
+    python -m georinex.read -q tests/demo3.10o  -u G -o r3G.nc
     """
     pytest.importorskip("netCDF4")
 
@@ -159,7 +159,7 @@ def test_one_system(use):
 
 def test_multi_system():
     """
-    ./ReadRinex.py -q tests/demo3.10o  -u G R -o r3GR.nc
+    python -m georinex.read -q tests/demo3.10o  -u G R -o r3GR.nc
     """
     pytest.importorskip("netCDF4")
 
@@ -173,7 +173,7 @@ def test_multi_system():
 
 def test_all_system():
     """
-    ./ReadRinex.py -q tests/demo3.10o -o r3all.nc
+    python -m georinex.read -q tests/demo3.10o -o r3all.nc
     """
     pytest.importorskip("netCDF4")
 
@@ -185,7 +185,7 @@ def test_all_system():
 
 def tests_all_indicators():
     """
-    ./ReadRinex.py -q tests/demo3.10o -useindicators -o r3all_indicators.nc
+    python -m georinex.read -q tests/demo3.10o -useindicators -o r3all_indicators.nc
     """
     pytest.importorskip("netCDF4")