public void GetCoordinates()
{
double jd = 2448724.5;
// geocentrical coordinates
CrdsEcliptical ecl = LunarMotion.GetCoordinates(jd);
Assert.AreEqual(133.162655, ecl.Lambda, 1e-6);
Assert.AreEqual(-3.229126, ecl.Beta, 1e-6);
Assert.AreEqual(368409.7, ecl.Distance, 1e-1);
// get nutation elements
var nutation = Nutation.NutationElements(jd);
// apparent geocentrical ecliptical coordinates
ecl += Nutation.NutationEffect(nutation.deltaPsi);
// true obliquity of the Earth orbit
double epsilon = Date.TrueObliquity(jd, nutation.deltaEpsilon);
// equatorial geocentrical coordinates
CrdsEquatorial eq = ecl.ToEquatorial(epsilon);
// Max error in Right Ascention is 0.1" of time
double errAlpha = new HMS("0h 0m 00.1s").ToDecimalAngle();
// Max error in Declination is 1" of arc
double errDelta = new DMS("0* 0' 01''").ToDecimalAngle();
Assert.AreEqual(new HMS("8h 58m 45.2s").ToDecimalAngle(), eq.Alpha, errAlpha);
Assert.AreEqual(new DMS("+13* 46' 06''").ToDecimalAngle(), eq.Delta, errDelta);
}
/// <summary>
/// Gets apparent geocentrical ecliptical coordinates of the Moon
/// </summary>
private CrdsEcliptical Ecliptical0(SkyContext c)
{
// geocentrical coordinates of the Moon
CrdsEcliptical ecliptical0 = LunarMotion.GetCoordinates(c.JulianDay);
// apparent geocentrical ecliptical coordinates
ecliptical0 += Nutation.NutationEffect(c.NutationElements.deltaPsi);
return(ecliptical0);
}