Changeset 8961 in ntrip for trunk/BNC/src/PPP/pppParlist.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/pppParlist.cpp (modified) (6 diffs)

trunk/BNC/src/PPP/pppParlist.cpp

@@ -141 +141 @@
   switch (_type) {
   case crdX:
-    if (tLC == t_lc::GIM) {return 0.0;}
+    if (tLC == t_lc::GIM || tLC == t_lc::Tz0) {return 0.0;}
     return (sta->xyzApr()[0] - obs->xc()[0]) / rhoV.NormFrobenius();
   case crdY:
-    if (tLC == t_lc::GIM) {return 0.0;}
+    if (tLC == t_lc::GIM || tLC == t_lc::Tz0) {return 0.0;}
     return (sta->xyzApr()[1] - obs->xc()[1]) / rhoV.NormFrobenius();
   case crdZ:
-    if (tLC == t_lc::GIM) {return 0.0;}
+    if (tLC == t_lc::GIM || tLC == t_lc::Tz0) {return 0.0;}
     return (sta->xyzApr()[2] - obs->xc()[2]) / rhoV.NormFrobenius();
   case clkR:
-    if (tLC == t_lc::GIM) {return 0.0;}
+    if (tLC == t_lc::GIM || tLC == t_lc::Tz0) {return 0.0;}
     return 1.0;
   case offGG:
-    if (tLC == t_lc::GIM) {return 0.0;}
+    if (tLC == t_lc::GIM || tLC == t_lc::Tz0) {return 0.0;}
     return (obs->prn().system() == 'R') ? 1.0 : 0.0;
   case amb:
-    if      (tLC == t_lc::GIM) {return 0.0;}
+    if      (tLC == t_lc::GIM || tLC == t_lc::Tz0) {return 0.0;}
     else if ((OPT->_obsModelType == OPT->IF)     ||
              (OPT->_obsModelType == OPT->PPPRTK) ||
@@ -181 +181 @@
     break;
   case trp:
-    if (tLC == t_lc::GIM) {return 0.0;}
-    return  1.0 / sin(obs->eleSat());
+    if      (tLC == t_lc::GIM) {
+      return 0.0;
+    }
+    else if (tLC == t_lc::Tz0) {
+      return 1.0;
+    }
+    else {
+      return 1.0 / sin(obs->eleSat());
+    }
   case ion:
-
     if (obs->prn() == _prn) {
       if      (tLC == t_lc::c1) {
@@ -203 +209 @@
     }
     if (tLC == t_lc::GIM && _prn == refPrn) {
-      return  1.0;
+      return 1.0;
     }
     break;
   case cBias1:
     if  (tLC == t_lc::c1) {
-      return  1.0;
+      return 1.0;
     }
     else {
@@ -216 +222 @@
   case cBias2:
      if (tLC == t_lc::c2) {
-      return  1.0;
+      return 1.0;
     }
     else {
@@ -224 +230 @@
   case pBias1:
     if  (tLC == t_lc::l1) {
-      return  1.0;
+      return 1.0;
     }
     else {
@@ -232 +238 @@
   case pBias2:
     if  (tLC == t_lc::l2) {
-      return  1.0;
+      return 1.0;
     }
     else {