* September 1998. This code is based upon algorithms and data kindly
* provided by Mr. Paul Schlyter. (pausch@saaf.se).
*
- * This program is free software; you can redistribute it and/or
- * modify it under the terms of the GNU General Public License as
- * published by the Free Software Foundation; either version 2 of the
- * License, or (at your option) any later version.
+ * This library is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU Library General Public
+ * License as published by the Free Software Foundation; either
+ * version 2 of the License, or (at your option) any later version.
*
- * This program is distributed in the hope that it will be useful, but
- * WITHOUT ANY WARRANTY; without even the implied warranty of
+ * This library is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
- * General Public License for more details.
+ * Library General Public License for more details.
*
- * You should have received a copy of the GNU General Public License
- * along with this program; if not, write to the Free Software
- * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
+ * You should have received a copy of the GNU Library General Public
+ * License along with this library; if not, write to the
+ * Free Software Foundation, Inc., 59 Temple Place - Suite 330,
+ * Boston, MA 02111-1307, USA.
*
* $Id$
**************************************************************************/
#endif
-#include <Time/fg_time.hxx>
-#include <Include/fg_constants.h>
+#include <simgear/constants.h>
class Star;
class CelestialBody
{
-protected: // make the data protected, in order to give the inherited
- // classes direct access to the data
+protected: // make the data protected, in order to give the
+ // inherited classes direct access to the data
double NFirst; /* longitude of the ascending node first part */
double NSec; /* longitude of the ascending node second part */
double iFirst; /* inclination to the ecliptic first part */
double magnitude;
double lonEcl, latEcl;
- double fgCalcEccAnom(double M, double e);
- double fgCalcActTime(FGTime *t);
- void updateOrbElements(FGTime *t);
+ double sgCalcEccAnom(double M, double e);
+ double sgCalcActTime(double mjd);
+ void updateOrbElements(double mjd);
public:
CelestialBody(double Nf, double Ns,
double wf, double ws,
double af, double as,
double ef, double es,
- double Mf, double Ms, FGTime *t);
+ double Mf, double Ms, double mjd);
+ CelestialBody(double Nf, double Ns,
+ double If, double Is,
+ double wf, double ws,
+ double af, double as,
+ double ef, double es,
+ double Mf, double Ms);
void getPos(double *ra, double *dec);
void getPos(double *ra, double *dec, double *magnitude);
+ double getRightAscension();
+ double getDeclination();
+ double getMagnitude();
double getLon();
double getLat();
- void updatePosition(FGTime *t, Star *ourSun);
+ void updatePosition(double mjd, Star *ourSun);
};
/*****************************************************************************
double wf, double ws,
double af, double as,
double ef, double es,
- double Mf, double Ms, FGTime *t)
+ double Mf, double Ms, double mjd)
{
NFirst = Nf; NSec = Ns;
iFirst = If; iSec = Is;
aFirst = af; aSec = as;
eFirst = ef; eSec = es;
MFirst = Mf; MSec = Ms;
- updateOrbElements(t);
-};
+ updateOrbElements(mjd);
+}
+
+inline CelestialBody::CelestialBody(double Nf, double Ns,
+ double If, double Is,
+ double wf, double ws,
+ double af, double as,
+ double ef, double es,
+ double Mf, double Ms)
+{
+ NFirst = Nf; NSec = Ns;
+ iFirst = If; iSec = Is;
+ wFirst = wf; wSec = ws;
+ aFirst = af; aSec = as;
+ eFirst = ef; eSec = es;
+ MFirst = Mf; MSec = Ms;
+}
/****************************************************************************
- * inline void CelestialBody::updateOrbElements(FGTime *t)
+ * inline void CelestialBody::updateOrbElements(double mjd)
* given the current time, this private member calculates the actual
* orbital elements
*
- * Arguments: FGTime *t: the current time:
+ * Arguments: double mjd: the current modified julian date:
*
* return value: none
***************************************************************************/
-inline void CelestialBody::updateOrbElements(FGTime *t)
+inline void CelestialBody::updateOrbElements(double mjd)
{
- double actTime = fgCalcActTime(t);
+ 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(FGTime *t)
+ * 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(FGTime *t)
+inline double CelestialBody::sgCalcActTime(double mjd)
{
- return (t->getMjd() - 36523.5);
+ return (mjd - 36523.5);
}
/*****************************************************************************
*magn = magnitude;
}
+inline double CelestialBody::getRightAscension() { return rightAscension; }
+inline double CelestialBody::getDeclination() { return declination; }
+inline double CelestialBody::getMagnitude() { return magnitude; }
+
inline double CelestialBody::getLon()
{
return lonEcl;