Changeset 9133 in ntrip for branches/BNC_2.12/src

Timestamp:

Sep 23, 2020, 2:56:18 PM (4 years ago)

Author:

stuerze

Message:

cleanup of the relativistic effects

File:

• branches/BNC_2.12/src/ephemeris.cpp (modified) (5 diffs)

branches/BNC_2.12/src/ephemeris.cpp

@@ -404 +404 @@
   // Relativistic Correction
   // -----------------------
-  // correspondent to IGS convention and GPS ICD (and SSR standard)
-  // xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;
   xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
-  // correspondent to IGS-SSR
+
   xc[4] = _clock_drift + _clock_driftrate*tc;
   xc[5] = _clock_driftrate;
@@ -429 +427 @@
   QString fmt = version < 3.0 ? "   %1%2%3%4\n" : "    %1%2%3%4\n";
 
-   out << QString(fmt)
+  out << QString(fmt)
     .arg(_IODE,    19, 'e', 12)
     .arg(_Crs,     19, 'e', 12)
@@ -644 +642 @@
   acc[1] = _y_acceleration * 1.e3;
   acc[2] = _z_acceleration * 1.e3;
+
   for (int ii = 1; ii <= nSteps; ii++) {
     _xv = rungeKutta4(_tt.gpssec(), _xv, step, acc, glo_deriv);
@@ -1000 +999 @@
   // Relativistic Correction
   // -----------------------
-  // correspondent to Galileo ICD and to SSR standard
   xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
-  // correspondent to IGS convention
-  //xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;
 
   xc[4] = _clock_drift + _clock_driftrate*tc;
@@ -1662 +1658 @@
   xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc;
 
-  // dotC  = _clock_drift + _clock_driftrate*tc
-  //       - 4.442807633e-10*_e*sqrt(a0)*cos(E) * dEdM * n;
+//  dotC  = _clock_drift + _clock_driftrate*tc
+//        - 4.442807633e-10*_e * sqrt(a0) * cos(E) * dEdM * n;
 
   // Relativistic Correction
   // -----------------------
-  // correspondent to BDS ICD and to SSR standard
   xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
-  // correspondent to IGS convention
-  // xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;
 
   xc[4] = _clock_drift + _clock_driftrate*tc;