Changeset 2429 in ntrip for trunk/BNS

Timestamp:

Apr 15, 2010, 1:42:35 PM (14 years ago)

Author:

mervart

Message:

* empty log message *

Location:

trunk/BNS

Files:

• bns.cpp (modified) (1 diff)
• bnseph.cpp (modified) (2 diffs)

trunk/BNS/bns.cpp

@@ -426 +426 @@
               xx(5) *= 1e-6;         // rel. corr.
                                      // xx(6), xx(7), xx(8) ... PhaseCent - CoM
-              xx(9)  /= 0.299792458; // xx(9)  ... P1-C1 DCB in ns
-              xx(10) /= 0.299792458; // xx(10) ... P1-P2 DCB in ns
+                                     // xx(9)  ... P1-C1 DCB in meters
+                                     // xx(10) ... P1-P2 DCB in meters
                                      // xx(11) ... dT
               xx(12) *= 1e3;         // x-crd at time + dT

trunk/BNS/bnseph.cpp

@@ -454 +454 @@
     tc += (GPSweek - _GPSweek) * secPerWeek;
   }
-  xc(4) = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc 
-          - 4.442807633e-10 * _e * sqrt(a0) *sin(E);
+  xc(4) = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc;
 
   // Velocity
@@ -480 +479 @@
 
   vv(3)  = sini    *doty  + yp*cosi      *doti;
+
+  // Relativistic Correction
+  // -----------------------
+  //  xc(4) -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
+  const static double c_vel = 299792458.0;
+  xc(4) -= 2.0 * DotProduct(xc.Rows(1,3),vv) / c_vel / c_vel;
 }