double magnitude;
double lonEcl, latEcl;
- double fgCalcEccAnom(double M, double e);
- double fgCalcActTime(double mjd);
+ double sgCalcEccAnom(double M, double e);
+ double sgCalcActTime(double mjd);
void updateOrbElements(double mjd);
public:
eFirst = ef; eSec = es;
MFirst = Mf; MSec = Ms;
updateOrbElements(mjd);
-};
+}
inline CelestialBody::CelestialBody(double Nf, double Ns,
double If, double Is,
aFirst = af; aSec = as;
eFirst = ef; eSec = es;
MFirst = Mf; MSec = Ms;
-};
+}
/****************************************************************************
* inline void CelestialBody::updateOrbElements(double mjd)
***************************************************************************/
inline void CelestialBody::updateOrbElements(double mjd)
{
- double actTime = fgCalcActTime(mjd);
+ double actTime = sgCalcActTime(mjd);
M = DEG_TO_RAD * (MFirst + (MSec * actTime));
w = DEG_TO_RAD * (wFirst + (wSec * actTime));
N = DEG_TO_RAD * (NFirst + (NSec * actTime));
a = aFirst + (aSec * actTime);
}
/*****************************************************************************
- * inline double CelestialBody::fgCalcActTime(double mjd)
+ * inline double CelestialBody::sgCalcActTime(double mjd)
* this private member function returns the offset in days from the epoch for
* wich the orbital elements are calculated (Jan, 1st, 2000).
*
*
* return value: the (fractional) number of days until Jan 1, 2000.
****************************************************************************/
-inline double CelestialBody::fgCalcActTime(double mjd)
+inline double CelestialBody::sgCalcActTime(double mjd)
{
return (mjd - 36523.5);
}