public void NutationInLongitudeTest()
{
var date = new AASDate(1987, 4, 10, 0, 0, 0, true);
double nutationInLongitude = AASNutation.NutationInLongitude(date.Julian);
Assert.Equal(-3.7879310766037446, nutationInLongitude);
}
protected void CorrectForFK5SystemReduction()
{
double appGeoLon = GetGeocentricLongitude(vectorToEarthCorrected.x, vectorToEarthCorrected.y);
double appGeoLat = GetGeocentricLatitude(vectorToEarthCorrected.x, vectorToEarthCorrected.y, vectorToEarthCorrected.z);
double deltaLambdaFK5 = AASFK5.CorrectionInLongitude(appGeoLon, appGeoLat, jdeCorrected);
double deltaBetaFK5 = AASFK5.CorrectionInLatitude(appGeoLon, jdeCorrected);
lambdaFK5 = lambdaAberration + deltaLambdaFK5;
betaFK5 = betaAberration + deltaBetaFK5;
//correct for nutation
double deltaPhi = AASNutation.NutationInLongitude(jdeCorrected) / 3600;
double deltaEpsilon = AASNutation.NutationInObliquity(jdeCorrected) / 3600;
epsilon = AASNutation.TrueObliquityOfEcliptic(jdeCorrected);
double lambdaNutation = lambdaFK5 + deltaPhi;
double betaNutation = betaFK5 + deltaEpsilon;
//x=RA in decimal hours, y=dec in degrees
AAS2DCoordinate equatorial = AASCoordinateTransformation.Ecliptic2Equatorial(lambdaNutation, betaNutation, epsilon);
equatorialCoords.RA = new HourAngle(equatorial.X);
equatorialCoords.Declination = new DegreesAngle(equatorial.Y);
}