Changeset 4908 in ntrip

Timestamp:

Feb 12, 2013, 11:29:38 AM (12 years ago)

Author:

weber

Message:

Clock correction due to heavy scale in GDA94 transformation

File:

• trunk/BNC/src/upload/bncrtnetuploadcaster.cpp (modified) (2 diffs)

trunk/BNC/src/upload/bncrtnetuploadcaster.cpp

@@ -574 +574 @@
                                     double& dc) {
 
-  // Clock Correction 
-  // ----------------
-  dc = 0.0; // TODO
-
   // Current epoch minus 2000.0 in years
   // ------------------------------------
@@ -595 +591 @@
 
   double sc = 1.0 + _sc * 1e-9 + dt * _scr * 1e-9;
+
+  // Clock correction due to heavy scale in GDA94 transformation
+  // -----------------------------------------------------------
+  if (_crdTrafo != "GDA94") {
+    dc = 0.0;
+  }
+  else {
+
+    // Spefify mean of Australia (here: GDA94 coordinates of site ALIC)
+    // ----------------------------------------------------------------
+    ColumnVector meanSta(3);
+    meanSta(1) = -4052051.767;
+    meanSta(2) =  4212836.197;
+    meanSta(3) = -2545106.022;
+
+    // Correction proportional to topocentric distance to satellites
+    // -------------------------------------------------------------
+    double rho = (xyz - meanSta).norm_Frobenius();
+    dc = rho * (sc - 1.0) / t_CST::c;
+  }
 
   Matrix rMat(3,3);