Index: trunk/BNC/bnctides.cpp
===================================================================
--- trunk/BNC/bnctides.cpp	(revision 2580)
+++ trunk/BNC/bnctides.cpp	(revision 2581)
@@ -174,6 +174,6 @@
 
   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,14 +216,10 @@
   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;