Observation callsign changes for GPSRO

src/opsinputs/GnssroStationIDMetOffice.cc

@@ -37,6 +37,7 @@ static ObsFunctionMaker<GnssroStationIDMetOffice>
 GnssroStationIDMetOffice::GnssroStationIDMetOffice(const eckit::LocalConfiguration & conf)
   : invars_() {
   invars_ += Variable("MetaData/satelliteIdentifier");
+  invars_ += Variable("MetaData/qualityFlags");
 }
 
 // -----------------------------------------------------------------------------
@@ -55,6 +56,20 @@ void GnssroStationIDMetOffice::compute(const ObsFilterData & in,
   std::vector<int> satid;
   in.get(Variable("MetaData/satelliteIdentifier"), satid);
 
+  const int missing = util::missingValue<int>();
+  std::vector<int> qualityFlags;
+  in.get(Variable("MetaData/qualityFlags"), qualityFlags);
+  std::vector<int> isRising(nlocs, 2);  // setting = 0, rising = 1, unknown = 2
+  for (size_t i = 0; i < nlocs; ++i) {
+    if (qualityFlags[i] == missing) {
+      isRising[i] = 2;
+    } else if ((qualityFlags[i] & (1 << 13)) != 0) {
+      isRising[i] = 1;
+    } else {
+      isRising[i] = 0;
+    }
+  }
+
   // Get the record number of each profile
   std::vector<size_t> recordNumbers = out.space().recidx_all_recnums();
   oops::Log::debug() << "Unique record numbers" << std::endl;
@@ -69,7 +84,8 @@ void GnssroStationIDMetOffice::compute(const ObsFilterData & in,
     for (size_t jobs : obsNumbers) {
       out[0][jobs] = stringFormat(satid[jobs], 4) +
                      stringFormat(iProfile, 5) +
-                     stringFormat(jobs, 7);
+                     stringFormat(jobs, 6) +
+                     stringFormat(isRising[jobs], 1);
     }
   }
 }

test/generate_unittest_netcdfs.py

@@ -852,7 +852,8 @@ def copy_var_to_var(Group, invarname, outvarname, filename):
                                   'ObsValue/bendingAngle',
                                   'ObsError/bendingAngle'],
                                  [],
-                                 ['MetaData/satelliteIdentifier'],
+                                 ['MetaData/satelliteIdentifier',
+                                  'MetaData/qualityFlags'],
                                  'testinput/varobs_globalnamelist_gnssro.nc4')
 
     # Sonde