Method Load

src/models/FGInertial.cpp:92–134

//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

bool FGInertial::Load(Element* el)
{
  if (!Upload(el, true)) return false;

  Name = el->GetAttributeValue("name");

  if (el->FindElement("semimajor_axis"))
    a = el->FindElementValueAsNumberConvertTo("semimajor_axis", "FT");
  else if (el->FindElement("equatorial_radius"))
    a = el->FindElementValueAsNumberConvertTo("equatorial_radius", "FT");
  if (el->FindElement("semiminor_axis"))
    b = el->FindElementValueAsNumberConvertTo("semiminor_axis", "FT");
  else if (el->FindElement("polar_radius"))
    b = el->FindElementValueAsNumberConvertTo("polar_radius", "FT");
  // Trigger GeographicLib exceptions if the equatorial or polar radii are
  // ill-defined.
  // This intercepts the exception before being thrown by a destructor.
  GeographicLib::Geodesic geod(a, 1.-b/a);

  if (el->FindElement("rotation_rate")) {
    double RotationRate = el->FindElementValueAsNumberConvertTo("rotation_rate", "RAD/SEC");
    vOmegaPlanet = {0., 0., RotationRate};
  }
  if (el->FindElement("GM"))
    GM = el->FindElementValueAsNumberConvertTo("GM", "FT3/SEC2");
  if (el->FindElement("J2"))
    J2 = el->FindElementValueAsNumber("J2"); // Dimensionless

  GroundCallback->SetEllipse(a, b);

  // Messages to warn the user about possible inconsistencies.
  if (debug_lvl > 0) {
    FGLogging log(LogLevel::WARN);
    if (a != b && J2 == 0.0)
      log << "Gravitational constant J2 is null for a non-spherical planet." << endl;
    if (a == b && J2 != 0.0)
      log << "Gravitational constant J2 is non-zero for a spherical planet." << endl;
  }

  Debug(2);

  return true;
}

//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%