Changeset 8961 in ntrip for trunk/BNC/src/PPP/pppSatObs.cpp

Timestamp:

Jun 29, 2020, 10:09:16 AM (4 years ago)

Author:

stuerze

Message:

PPP update: pseudo obs tropo added

File:

• trunk/BNC/src/PPP/pppSatObs.cpp (modified) (7 diffs)

trunk/BNC/src/PPP/pppSatObs.cpp

@@ -92 +92 @@
   for (unsigned ii = 0; ii < OPT->LCs(_prn.system()).size(); ii++) {
     t_lc::type tLC = OPT->LCs(_prn.system())[ii];
-    if (tLC == t_lc::GIM) {continue;}
+    if (tLC == t_lc::GIM || tLC == t_lc::Tz0) {continue;}
     if (!isValid(tLC)) {
       _valid = false;
@@ -190 +190 @@
     return;
   case t_lc::GIM:
+  case t_lc::Tz0:
   case t_lc::dummy:
   case t_lc::maxLc:
@@ -218 +219 @@
     }
     else {
-      return _stecRefSat - _stecSat;
-    }
+      return _stecRefSat;
+    }
+  }
+
+  if (tLC == t_lc::Tz0) {
+    return _model._tropo0;
   }
 
@@ -292 +297 @@
   if (tLC == t_lc::GIM) {
     retVal = OPT->_sigmaGIMdiff * OPT->_sigmaGIMdiff;
+  }
+
+  if (tLC == t_lc::Tz0) {
+    retVal = OPT->_sigmaTz0 * OPT->_sigmaTz0;
   }
 
@@ -407 +416 @@
   // Tropospheric Delay
   // ------------------
-  _model._tropo = t_tropo::delay_saast(rRec, _model._eleSat);
+  _model._tropo  = t_tropo::delay_saast(rRec, _model._eleSat);
+  _model._tropo0 = t_tropo::delay_saast(rRec, M_PI/2.0);
 
   // Code Biases
@@ -570 +580 @@
   }
   else if (tLC == t_lc::GIM) {
-    cmpValue = 0.0;
+    cmpValue =  _stecSat;
+  }
+  else if (tLC == t_lc::Tz0) {
+    cmpValue = _model._tropo0;
   }
   else {
@@ -633 +646 @@
   _stecRefSat = stecRefSat;
 }
+
+
+//
+////////////////////////////////////////////////////////////////////////////
+void  t_pppSatObs::setPseudoObsTropo() {qDebug() << "setPseudoObsTropo";
+  _tropo0 = _model._tropo0;
+}