From f21109960d9d5bd135a98e789bb612abf9961421 Mon Sep 17 00:00:00 2001 From: Daniel Warner Date: Sat, 9 Apr 2011 19:56:35 +0100 Subject: [PATCH] Comments --- Eci.cpp | 16 +++++++++++++++- 1 file changed, 15 insertions(+), 1 deletion(-) diff --git a/Eci.cpp b/Eci.cpp index 266a3f5..6ddffb2 100644 --- a/Eci.cpp +++ b/Eci.cpp @@ -13,16 +13,30 @@ Eci::Eci(const Julian &date, const CoordGeodetic &geo) */ const double theta = date_.ToLocalMeanSiderealTime(longitude); + /* + * take into account earth flattening + */ const double c = 1.0 / sqrt(1.0 + Globals::F() * (Globals::F() - 2.0) * pow(sin(latitude), 2.0)); const double s = pow(1.0 - Globals::F(), 2.0) * c; const double achcp = (Globals::XKMPER() * c + altitude) * cos(latitude); + /* + * X position in km + * Y position in km + * Z position in km + * W magnitude in km + */ position_.SetX(achcp * cos(theta)); position_.SetY(achcp * sin(theta)); - position_.SetZ((Globals::XKMPER() * s + altitude) * sin(latitude)); position_.SetW(position_.GetMagnitude()); + /* + * X velocity in km/s + * Y velocity in km/s + * Z velocity in km/s + * W magnitude in km/s + */ velocity_.SetX(-mfactor * position_.GetY()); velocity_.SetY(mfactor * position_.GetX()); velocity_.SetZ(0.0);