- Timestamp:
- Aug 29, 2010, 3:06:26 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified trunk/BNC/bnctides.cpp ¶
r2579 r2581 174 174 175 175 static double lastMjd = 0.0; 176 static ColumnVector xSun (3);177 static ColumnVector xMoon (3);176 static ColumnVector xSun; 177 static ColumnVector xMoon; 178 178 static double rSun; 179 179 static double rMoon; … … 216 216 double facSun = gms / gmWGS * 217 217 (rRec * rRec * rRec * rRec) / (rSun * rSun * rSun); 218 218 219 double facMoon = gmm / gmWGS * 219 220 (rRec * rRec * rRec * rRec) / (rMoon * rMoon * rMoon); 220 221 221 double theta = GMST(Mjd);222 223 double Ell[3]; xyz2ell(xyz.data(), Ell);224 225 222 ColumnVector dX = facSun * (x2Sun * xSun + p2Sun * xyzUnit) + 226 facMoon * (x2Moon * xMoon + p2Moon * xyzUnit) - 227 0.025 * sin(Ell[0]) * cos(Ell[0]) * sin(theta+Ell[1]) * xyzUnit; 223 facMoon * (x2Moon * xMoon + p2Moon * xyzUnit); 228 224 229 225 xyz += dX;
Note:
See TracChangeset
for help on using the changeset viewer.