spacecruft
/
CelestiaContent
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Accept either MeanAnomaly or MeanLongitude in orbital elements.

pull/3/head
Chris Laurel 2001-03-21 22:39:16 +00:00
parent 2c08c63e40
commit 5dbbc4424d
1 changed files with 19 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
src/solarsys.cpp
View File

@ -153,25 +153,36 @@ static Body* CreatePlanet(PlanetarySystem* system,
double ascendingNode = 0.0;
planetData->getNumber("AscendingNode", ascendingNode);
double argOfPeriapsis = 0.0;
planetData->getNumber("ArgOfPeriapsis", argOfPeriapsis);
double longOfPericenter = 0.0;
planetData->getNumber("ArgOfPeriapsis", longOfPericenter);
double trueAnomaly = 0.0;
planetData->getNumber("TrueAnomaly", trueAnomaly);
double epoch = astro::J2000;
planetData->getNumber("Epoch", epoch);
// Accept either the mean anomaly or mean longitude--use mean anomaly
// if both are specified.
double anomalyAtEpoch = 0.0;
if (!planetData->getNumber("MeanAnomaly", anomalyAtEpoch))
{
double longAtEpoch = 0.0;
if (planetData->getNumber("MeanLongitude", longAtEpoch))
anomalyAtEpoch = longAtEpoch - longOfPericenter;
}
if (usePlanetUnits)
{
semiMajorAxis = astro::AUtoKilometers(semiMajorAxis);
period = period * 265.25f;
period = period * 365.25f;
}
body->setOrbit(new EllipticalOrbit(semiMajorAxis,
eccentricity,
degToRad(inclination),
degToRad(ascendingNode),
degToRad(argOfPeriapsis),
degToRad(trueAnomaly),
period));
degToRad(longOfPericenter),
degToRad(anomalyAtEpoch),
period,
epoch));
double obliquity = 0.0;
planetData->getNumber("Obliquity", obliquity);