Changeset 2581 in ntrip

Timestamp:

Aug 29, 2010, 3:06:26 PM (14 years ago)

Author:

mervart

Message:

 

File:

• trunk/BNC/bnctides.cpp (modified) (2 diffs)

trunk/BNC/bnctides.cpp

@@ -174 +174 @@
 
   static double       lastMjd = 0.0;
-  static ColumnVector xSun(3);
-  static ColumnVector xMoon(3);
+  static ColumnVector xSun;
+  static ColumnVector xMoon;
   static double       rSun;
   static double       rMoon;
@@ -216 +216 @@
   double facSun  = gms / gmWGS * 
                    (rRec * rRec * rRec * rRec) / (rSun * rSun * rSun);
+
   double facMoon = gmm / gmWGS * 
                    (rRec * rRec * rRec * rRec) / (rMoon * rMoon * rMoon);
 
-  double theta = GMST(Mjd);
-
-  double Ell[3]; xyz2ell(xyz.data(), Ell);
-
   ColumnVector dX = facSun  * (x2Sun  * xSun  + p2Sun  * xyzUnit) + 
-                    facMoon * (x2Moon * xMoon + p2Moon * xyzUnit) -
-              0.025 * sin(Ell[0]) * cos(Ell[0]) * sin(theta+Ell[1]) * xyzUnit;
+                    facMoon * (x2Moon * xMoon + p2Moon * xyzUnit);
 
   xyz += dX;