Changeset 7252 in ntrip

Timestamp:

Aug 24, 2015, 6:01:02 PM (9 years ago)

Author:

stuerze

Message:

ionospheric delay is added to the dispersive model part

File:

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

trunk/BNC/src/PPP/pppSatObs.cpp

@@ -396 +396 @@
   // -----------------
   const t_vTec* vTec = PPP_CLIENT->obsPool()->vTec();
-  if (vTec) {
+  if (vTec && (!isValid(t_lc::cIF) || !isValid(t_lc::lIF))) {
     double stec = station->stec(vTec, _signalPropagationTime, rSat);
     for (unsigned iFreq = 1; iFreq < t_frequency::max; iFreq++) {
@@ -411 +411 @@
   // ------------------
   _model._set = true;
+
+  //printModel();
 
   return success;
@@ -432 +434 @@
   for (unsigned iFreq = 1; iFreq < t_frequency::max; iFreq++) {
     if (_obs[iFreq]) {
-      LOG << "PCO:        " << t_frequency::toString(t_frequency::type(iFreq)) << setw(12) << setprecision(3) << _model._antPCO[iFreq]    << endl
-          << "BIAS CODE:  " << t_frequency::toString(t_frequency::type(iFreq)) << setw(12) << setprecision(3) << _model._codeBias[iFreq]  << endl
-          << "BIAS PHASE: " << t_frequency::toString(t_frequency::type(iFreq)) << setw(12) << setprecision(3) << _model._phaseBias[iFreq] << endl;
+      LOG << "PCO           : " << t_frequency::toString(t_frequency::type(iFreq)) << setw(12) << setprecision(3) << _model._antPCO[iFreq]    << endl
+          << "BIAS CODE     : " << t_frequency::toString(t_frequency::type(iFreq)) << setw(12) << setprecision(3) << _model._codeBias[iFreq]  << endl
+          << "BIAS PHASE    : " << t_frequency::toString(t_frequency::type(iFreq)) << setw(12) << setprecision(3) << _model._phaseBias[iFreq]  << endl
+          << "IONO CODEDELAY: " << t_frequency::toString(t_frequency::type(iFreq)) << setw(12) << setprecision(3) << _model._ionoCodeDelay[iFreq] << endl;
     }
   }
@@ -484 +487 @@
   for (it = codeCoeff.begin(); it != codeCoeff.end(); it++) {
     t_frequency::type tFreq = it->first;
-    dispPart += it->second * (_model._antPCO[tFreq] + _model._codeBias[tFreq]);
+    dispPart += it->second * (_model._antPCO[tFreq] + _model._codeBias[tFreq] +
+                              _model._ionoCodeDelay[tFreq]);
   }
   for (it = phaseCoeff.begin(); it != phaseCoeff.end(); it++) {
     t_frequency::type tFreq = it->first;
     dispPart += it->second * (_model._antPCO[tFreq] + _model._phaseBias[tFreq] +
-                              _model._windUp * t_CST::lambda(tFreq, _channel));
+                              _model._windUp * t_CST::lambda(tFreq, _channel)  -
+                              _model._ionoCodeDelay[tFreq]);
   }