Changeset 5752 in ntrip for trunk/BNC/src/PPP/parlist.cpp

Timestamp:

Aug 3, 2014, 11:52:30 AM (10 years ago)

Author:

mervart

Message:

 

File:

• trunk/BNC/src/PPP/parlist.cpp (modified) (12 diffs)

trunk/BNC/src/PPP/parlist.cpp

@@ -7 +7 @@
 #include "parlist.h"
 #include "satobs.h"
+
 #include "station.h"
-#include "tropo.h"
-#include "iono.h"
-#include "utils.h"
-#include "genconst.h"
+#include "bncutils.h"
+#include "bncconst.h"
+#include "pppClient.h"
 
 using namespace BNC;
@@ -33 +33 @@
    case crdX:
      _epoSpec = true;
-     _sigma0  = OPT->sigmaCrd();
+     _sigma0  = OPT->_sigCrd[0];
      break;
    case crdY:
      _epoSpec = true;
-     _sigma0  = OPT->sigmaCrd();
+     _sigma0  = OPT->_sigCrd[1];
      break;
    case crdZ:
      _epoSpec = true;
-     _sigma0  = OPT->sigmaCrd();
+     _sigma0  = OPT->_sigCrd[2];
      break;
    case clkR:
@@ -57 +57 @@
              offGG = PPP_CLIENT->offGG();
            }
-           _x0 = nint((obs->obsValue(tLC) - offGG - obs->cmpValue(tLC)) / obs->lambda(tLC));
+           _x0 = floor((obs->obsValue(tLC) - offGG - obs->cmpValue(tLC)) / obs->lambda(tLC) + 0.5);
            break;
          }
@@ -72 +72 @@
    case trp:
      _epoSpec = false;
-     _sigma0  = OPT->sigmaTropo();
-     _noise   = OPT->noiseTropo();
-     break;
-   case iono:
-     if (OPT->noiseIono() > 0.0) {
-       _epoSpec = false;
-     }
-     else {
-       _epoSpec = true;
-     }
-     _sigma0  = OPT->sigmaIono();
-     _noise   = OPT->noiseIono();
+     _sigma0  = OPT->_sigTropo;
+     _noise   = OPT->_noiseTropo;
      break;
    case bias:
@@ -100 +90 @@
 // 
 ////////////////////////////////////////////////////////////////////////////
-double t_param::partial(const bncTime& epoTime, const t_satObs* obs, 
+double t_param::partial(const bncTime& /* epoTime */, const t_satObs* obs, 
                         const t_lc::type& tLC) const {
 
@@ -144 +134 @@
     return 0.0;
   case trp:
-    return t_tropo::mf(OPT->tropoMF(), epoTime, sta->ellApr()[0],
-                       sta->ellApr()[1], sta->ellApr()[2], obs->eleSat());
-  case iono:
-    if (obs->prn() == _prn) {
-      return t_iono::dLCdI(tLC, obs->prn().system(), obs->channel());
-    }
-    else {
-      return 0.0;
-    }
+    return 1.0 / sin(obs->eleSat()); 
   case bias:
     if (tLC == _tLC && obs->prn().system() == _prn.system()) {
@@ -191 +173 @@
     ss << "TRP        ";
     break;
-  case iono:
-    ss << "IONO " << _prn.toString();
-    break;
   case bias:
     ss << "BIAS " << _prn.system() << ' ' << left << setw(3) << t_lc::toString(_tLC);
@@ -237 +216 @@
     }
 
-    else if (par->type() == t_param::amb || par->type() == t_param::iono) {
+    else if (par->type() == t_param::amb) {
       if (par->lastObsTime().valid() && (epoTime - par->lastObsTime() > 3600.0)) {
         remove = true;
@@ -298 +277 @@
   // Receiver Biases
   // ---------------
-  for (unsigned ii = 0; ii < OPT->estBias().size(); ii++) {
-    const t_options::t_optBias& optBias = OPT->estBias()[ii];
+  for (unsigned ii = 0; ii < OPT->_estBias.size(); ii++) {
+    const t_options::t_optBias& optBias = OPT->_estBias[ii];
     required.push_back(new t_param(t_param::bias, t_prn(optBias._system, 1), optBias._tLC));
   }
@@ -316 +295 @@
       const t_satObs* satObs = obsVector[jj];
       required.push_back(new t_param(t_param::amb, satObs->prn(), tLC, &obsVector));
-    }
-  }
-
-  // Stochastic Ionosphere
-  //----------------------
-  if (OPT->estIono()) {
-    for (unsigned jj = 0; jj < obsVector.size(); jj++) {
-      const t_satObs* satObs = obsVector[jj];
-      required.push_back(new t_param(t_param::iono, satObs->prn(), t_lc::dummy));
     }
   }
@@ -365 +335 @@
   }
 
-  return t_irc::success;
+  return success;
 }
 
@@ -398 +368 @@
           << setw(8)  << setprecision(4) << sqrt(QQ[ind][ind]);
       if (par->type() == t_param::amb) {
-        LOG << " el = " << setw(6) << setprecision(2) << par->ambEleSat() * t_genConst::rho_deg 
+        LOG << " el = " << setw(6) << setprecision(2) << par->ambEleSat() * 180.0 / M_PI
             << " epo = " << setw(4) << par->ambNumEpo();
       }