Changeset 8965 in ntrip for trunk/BNC/src/PPP

Timestamp:

Jun 30, 2020, 3:47:40 PM (4 years ago)

Author:

stuerze

Message:

minor changes

Location:

trunk/BNC/src/PPP

Files:

• pppFilter.cpp (modified) (8 diffs)
• pppSatObs.cpp (modified) (1 diff)

trunk/BNC/src/PPP/pppFilter.cpp

@@ -118 +118 @@
   // ----------------------------------------
   if ((OPT->_obsModelType == OPT->DCMcodeBias   ||
-       OPT->_obsModelType == OPT->DCMphaseBias) && (_numEpoProcessing == 1)) {
+       OPT->_obsModelType == OPT->DCMphaseBias) &&(_numEpoProcessing == 1))  {
     _numAllUsedLCs = 0;
     for (unsigned iSys = 0; iSys < OPT->systems().size(); iSys++) {
@@ -125 +125 @@
       if (OPT->_pseudoObsIono && epoch->pseudoObsIono() == false) {
         _numAllUsedLCs -= 1;  // GIM not used
-      }
-    }
-    int modify = 0;
-    if (OPT->_pseudoObsTropo) {
-      modify  = -1;
-    }
-    // max Obs
-    int maxObs = allObs.size() * (_numAllUsedLCs + modify);
-    if (OPT->_pseudoObsIono && epoch->pseudoObsIono() == true) { // stecDiff w.r.t refSat
-      maxObs -= 1;
-    }
-    if (OPT->_pseudoObsTropo) {
-      maxObs += 1;  // once per station
-    }
-    _datumTrafo->initAA(maxObs, _parlist->nPar());
+        if (OPT->_pseudoObsTropo) { // if no VTEC values /pseudo obs iono then no pseudo obs tropo
+          _numAllUsedLCs -= 1;
+        }
+      }
+      int modifyLCs = 0;
+      if (OPT->_pseudoObsTropo &&
+          OPT->_pseudoObsIono && epoch->pseudoObsIono() == true) {
+        modifyLCs  = -1;
+      }
+      // max Obs
+      int maxObs = allObs.size() * (_numAllUsedLCs + modifyLCs);
+      if (OPT->_pseudoObsIono && epoch->pseudoObsIono() == true) {
+        maxObs -= 1; // pseudo obs iono with respect to refSat
+      }
+      if (OPT->_pseudoObsIono && epoch->pseudoObsIono() == true &&
+          OPT->_pseudoObsTropo) {
+        maxObs += 1;   // pseudo obs tropo once per station
+      }
+      _datumTrafo->initAA(maxObs, _parlist->nPar());
+    }
   }
 
@@ -225 +230 @@
   // Detect Cycle Slips
   // ------------------
-
   if (detectCycleSlips(LCs, obsVector, refPrn, preProcessing) != success) {
     return failure;
@@ -231 +235 @@
 
   unsigned usedLCs = LCs.size();
-  if (OPT->_pseudoObsIono && !pseudoObsIonoAvailable) {
-    usedLCs -= 1;  // GIM not used
-  }
 
   ColumnVector               xSav       = _xFlt;
@@ -239 +240 @@
   string                     epoTimeStr = string(_epoTime);
   const vector<t_pppParam*>& params     = _parlist->params();
-  int modify = 0;
-  if (OPT->_pseudoObsTropo) {
-    modify  = -1;
+
+  if (OPT->_pseudoObsIono && !pseudoObsIonoAvailable) {
+    usedLCs -= 1;  // GIM not used
+    if (OPT->_pseudoObsTropo) { // if no VTEC values /pseudo obs iono then no pseudo obs tropo
+      usedLCs -= 1;
+    }
+  }
+
+  int modifyLCs = 0;
+  if (OPT->_pseudoObsTropo &&
+      OPT->_pseudoObsIono && pseudoObsIonoAvailable) {
+    modifyLCs  = -1;
   }
   // max Obs
-  unsigned maxObs = obsVector.size() * (usedLCs + modify);
-  if (OPT->_pseudoObsIono && pseudoObsIonoAvailable) { // stecDiff w.r.t refSat
-    maxObs -= 1;
-  }
-  if (OPT->_pseudoObsTropo) {
-    maxObs +=1;   // tropo pseudo obs once per station
+  unsigned maxObs = obsVector.size() * (usedLCs + modifyLCs);
+  if (OPT->_pseudoObsIono && pseudoObsIonoAvailable) {
+    maxObs -= 1; // pseudo obs iono with respect to refSat
+  }
+  if (OPT->_pseudoObsIono && pseudoObsIonoAvailable &&
+      OPT->_pseudoObsTropo) {
+    maxObs += 1;   // pseudo obs tropo once per station
   }
 
@@ -272 +283 @@
     for (unsigned ii = 0; ii < obsVector.size(); ii++) {
       t_pppSatObs* obs = obsVector[ii];
+      unsigned hlp = ii+1;
       if (!obs->outlier()) {
         for (unsigned jj = 0; jj < usedLCs; jj++) {
           const t_lc::type tLC = LCs[jj];
-          if (tLC == t_lc::GIM && obs->isReference()) {continue;}
-          if (tLC == t_lc::Tz0 && ii+1 != obsVector.size()) {continue;}
+          if (tLC == t_lc::GIM &&  obs->isReference()) {continue;}
+          if (tLC == t_lc::Tz0 && hlp != obsVector.size()) {continue;}
           ++iObs;
           usedObs.push_back(obs);
@@ -465 +477 @@
           double ll = obs->obsValue(tLC) - obs->cmpValue(tLC) - DotProduct(_x0, AA);
           double vv = DotProduct(AA, _xFlt) - ll;
+
           if (fabs(vv) > SLIP) {
             if (preProcessing) {
@@ -486 +499 @@
 t_irc t_pppFilter::resetAmb(t_prn prn, const vector<t_pppSatObs*>& obsVector,
                             SymmetricMatrix* QSav, ColumnVector* xSav) {
-                            
+
   t_irc irc = failure;
   vector<t_pppParam*>& params = _parlist->params();

trunk/BNC/src/PPP/pppSatObs.cpp

@@ -650 +650 @@
 //
 ////////////////////////////////////////////////////////////////////////////
-void  t_pppSatObs::setPseudoObsTropo() {qDebug() << "setPseudoObsTropo";
+void  t_pppSatObs::setPseudoObsTropo() {
   _tropo0 = _model._tropo0;
 }