add more orbit sources to nmea2snr

CHANGELOG.md

@@ -7,8 +7,11 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/).
 Fixed bugs in highrate RINEX 3 downloads/translations that conflicted with parallel processing
 (GA and BKG archives). This was mostly related to how I stored interim Rinex 2 files.
 
+added ultra rapid orbits to nmea code
+
+
 ## 3.1.1
-?force update of user manual
+nmea2snr checks rapid, final, and ultra GNSS orbits from GFZ.
 
 ## 3.1.0
 

README.md

@@ -1,8 +1,9 @@
-# gnssrefl v3.1.0 
+# gnssrefl v3.1.1 
 
 [![PyPI Version](https://img.shields.io/pypi/v/gnssrefl.svg)](https://pypi.python.org/pypi/gnssrefl) 
 
-[![DOI](https://zenodo.org/badge/doi/10.5281/zenodo.10644225.svg)](http://dx.doi.org/10.5281/zenodo.10644225) 
+[![DOI](https://zenodo.org/badge/doi/10.5281/zenodo.5601494.svg)](http://dx.doi.org/10.5281/zenodo.5601494) 
+
 
 gnssrefl is an open source software package for GNSS interferometric reflectometry (GNSS-IR). 
 

gnssrefl/nmea2snr.py

@@ -101,14 +101,20 @@ def nmea_translate(locdir, fname, snrfile, csnr, dec, year, doy, sp3, recv, gzip
     if sp3:
         # first try to find precise because they have beidou
         xf,orbdir,foundit=g.gbm_orbits_direct(year,month,day)
+        # try to find it elsewhere
         if not foundit: 
             print('Could not find the precise GNSS orbits from GFZ. ')
             xf,orbdir,foundit=g.rapid_gfz_orbits(year,month,day)
             if not foundit: 
-                print('Could not find the rapid orbits from GFZ. Exiting')
-            return
-        else:
-            orbfile = orbdir + '/' + xf # hopefully
+                print('Could not find the rapid orbits from GFZ. ')
+                xf,orbdir,foundit = g.ultra_gfz_orbits(year,month,day,0)
+                if not foundit:
+                    print('Could not find the ultrarapid orbits from GFZ. Exiting')
+                    return
+
+        # define orbit file name
+        orbfile = orbdir + '/' + xf # hopefully
+        print('Found a sp3 file: ', orbfile)
 
     # this is to help a colleague 
     if station == 'argt':

gnssrefl/sd_libs.py

@@ -996,8 +996,13 @@ def RH_ortho_plot2( station, H0, year,  txtdir, fs, time_rh, rh, gap_min_val,th,
 
     # looks like I identified the gaps in day of year units - 
     # but then did the implementation in mjd and then datetime ...
-    splinefileout =  txtdir + '/' + station + '_' + str(year) + '_spline_out.txt'
-    splinefileout =  txtdir + '/' + station +  '_spline_out.txt'
+    #splinefileout =  txtdir + '/' + station + '_' + str(year) + '_spline_out.txt'
+    if csvfile:
+        splinefileout =  txtdir + '/' + station +  '_spline_out.csv'
+    else:
+        splinefileout =  txtdir + '/' + station +  '_spline_out.txt'
+
+
     print('Writing evenly sampled file to: ', splinefileout)
     fout = open(splinefileout,'w+')
     vn = station + ' gnssrefl v' + str(g.version('gnssrefl'))

setup.py

@@ -42,7 +42,7 @@
 ]
 setup(
     name="gnssrefl",
-    version="3.1.0",
+    version="3.1.1",
     author="Kristine Larson",
     author_email="[email protected]",
     description="A GNSS reflectometry software package ",