Changeset 9322 in ntrip for branches/BNC_2.12/src/bncephuser.cpp

Timestamp:

Dec 17, 2020, 4:11:09 PM (4 years ago)

Author:

stuerze

Message:

some lines are added to consider corrected GPS BRDC data sets which come 16 seconds before nominal TOC

File:

• branches/BNC_2.12/src/bncephuser.cpp (modified) (6 diffs)

branches/BNC_2.12/src/bncephuser.cpp

@@ -158 +158 @@
   }
 
-  if (ephOld == 0 || newEph->isNewerThan(ephOld)) {
+  if ( (ephOld == 0) ||
+       newEph->isNewerThan(ephOld) ||
+       newCorrectedGpsEphSet(newEph)  ) {
     checkEphemeris(newEph, realTime);
     eph->setCheckState(newEph->checkState());
@@ -170 +172 @@
       newEph->checkState() != t_eph::outdated) {
     deque<t_eph*>& qq = _eph[prn];
-    if (qq.empty() ||
-        newEph->isNewerThan(qq.back())){
-      qq.push_back(newEph);
-      if (qq.size() > _maxQueueSize) {
-        delete qq.front();
-        qq.pop_front();
-      }
-      ephBufferChanged();
-      return success;
-    }
-    else {
-      delete newEph;
-      return failure;
-    }
-
+    qq.push_back(newEph);
+    if (qq.size() > _maxQueueSize) {
+      delete qq.front();
+      qq.pop_front();
+    }
+    ephBufferChanged();
+    return success;
   }
   else {
@@ -280 +274 @@
     ColumnVector vvL(3);
     ephL->getCrd(eph->TOC(), xcL, vvL, false);
-    double dt = eph->TOC() - ephL->TOC();
-    if (dt < 0.0) {
-      dt += 604800.0;
-    }
+
     double diff  = (xc.Rows(1,3) - xcL.Rows(1,3)).norm_Frobenius();
     double diffC = fabs(xc(4) - xcL(4)) * t_CST::c;
 
     if (diff < MAXDIFF && diffC < MAXDIFF) {
-      if (diff < MINDIFF && diffC < MINDIFF && ephL->checkState() == t_eph::ok) {
+      if (diff < MINDIFF && diffC < MINDIFF &&
+          ephL->checkState() == t_eph::ok) {
         // to prevent the same data sets with different TOC values
         eph->setCheckState(t_eph::bad);
@@ -302 +294 @@
       }
     }
+
+    double dt = eph->TOC() - ephL->TOC();
     // some lines to allow update of ephemeris data sets after an outage
+    // -----------------------------------------------------------------
     // update interval: 2h,
+    // sometimes corrected BRDC sets are sent with TOC values 16 seconds
+    // before the nominal TOC (integer time)
     if      (eph->type() == t_eph::GPS     && dt > (2*3600.0 + 60.0)) {
       ephL->setCheckState(t_eph::outdated);
@@ -314 +311 @@
     }
     // updated every 30 minutes
-    else if (eph->type() == t_eph::GLONASS && dt > (1800.0 + 60.0)) {
+    else if (eph->type() == t_eph::GLONASS && dt > (  1800.0 + 60.0)) {
       ephL->setCheckState(t_eph::outdated);
       return;
@@ -340 +337 @@
   }
 }
+
+//
+////////////////////////////////////////////////////////////////////////////
+bool bncEphUser::newCorrectedGpsEphSet(const t_eph* eph) {
+
+  bool correctedGpsEphSet = false;
+  bool newCorrectedGpsEphSet = false;
+
+  if (eph->prn().system() != 'G') {
+    return newCorrectedGpsEphSet;
+  }
+
+  if ((fmod(eph->TOC().daysec()+16.0, 7200) == 0.0)) {
+    // sometimes corrected BRDC sets are sent with TOC values 16 seconds
+    // before the nominal TOC (integer time)
+    correctedGpsEphSet = true;
+  }
+
+  QString      prn     = QString(eph->prn().toInternalString().c_str());
+  t_eph*       ephL    = ephLast(prn);
+  if (!ephL) { // first set for GPS Sat
+    return correctedGpsEphSet;
+  }
+
+
+  double dt = eph->TOC() - ephL->TOC();
+  if (correctedGpsEphSet){
+    if (dt == 0.0) {
+      newCorrectedGpsEphSet = false;
+    }
+    else {
+      newCorrectedGpsEphSet = true;
+    }
+  }
+  return  newCorrectedGpsEphSet;
+}