/// <summary>
/// Calculate approximate position of the sun for a local date and time.
/// </summary>
/// <param name="lctHours">Local civil time, in hours.</param>
/// <param name="lctMinutes">Local civil time, in minutes.</param>
/// <param name="lctSeconds">Local civil time, in seconds.</param>
/// <param name="localDay">Local day, day part.</param>
/// <param name="localMonth">Local day, month part.</param>
/// <param name="localYear">Local day, year part.</param>
/// <param name="isDaylightSaving">Is daylight savings in effect?</param>
/// <param name="zoneCorrection">Time zone correction, in hours.</param>
/// <returns>
/// <para>sunRAHour -- Right Ascension of Sun, hour part</para>
/// <para>sunRAMin -- Right Ascension of Sun, minutes part</para>
/// <para>sunRASec -- Right Ascension of Sun, seconds part</para>
/// <para>sunDecDeg -- Declination of Sun, degrees part</para>
/// <para>sunDecMin -- Declination of Sun, minutes part</para>
/// <para>sunDecSec -- Declination of Sun, seconds part</para>
/// </returns>
public (double sunRAHour, double sunRAMin, double sunRASec, double sunDecDeg, double sunDecMin, double sunDecSec) ApproximatePositionOfSun(double lctHours, double lctMinutes, double lctSeconds, double localDay, int localMonth, int localYear, bool isDaylightSaving, int zoneCorrection)
{
var daylightSaving = (isDaylightSaving == true) ? 1 : 0;
var greenwichDateDay = PAMacros.LocalCivilTimeGreenwichDay(lctHours, lctMinutes, lctSeconds, daylightSaving, zoneCorrection, localDay, localMonth, localYear);
var greenwichDateMonth = PAMacros.LocalCivilTimeGreenwichMonth(lctHours, lctMinutes, lctSeconds, daylightSaving, zoneCorrection, localDay, localMonth, localYear);
var greenwichDateYear = PAMacros.LocalCivilTimeGreenwichYear(lctHours, lctMinutes, lctSeconds, daylightSaving, zoneCorrection, localDay, localMonth, localYear);
var utHours = PAMacros.LocalCivilTimeToUniversalTime(lctHours, lctMinutes, lctSeconds, daylightSaving, zoneCorrection, localDay, localMonth, localYear);
var utDays = utHours / 24;
var jdDays = PAMacros.CivilDateToJulianDate(greenwichDateDay, greenwichDateMonth, greenwichDateYear) + utDays;
var dDays = jdDays - PAMacros.CivilDateToJulianDate(0, 1, 2010);
var nDeg = 360 * dDays / 365.242191;
var mDeg1 = nDeg + PAMacros.SunELong(0, 1, 2010) - PAMacros.SunPeri(0, 1, 2010);
var mDeg2 = mDeg1 - 360 * (mDeg1 / 360).Floor();
var eCDeg = 360 * PAMacros.SunEcc(0, 1, 2010) * mDeg2.ToRadians().Sine() / Math.PI;
var lSDeg1 = nDeg + eCDeg + PAMacros.SunELong(0, 1, 2010);
var lSDeg2 = lSDeg1 - 360 * (lSDeg1 / 360).Floor();
var raDeg = PAMacros.EcRA(lSDeg2, 0, 0, 0, 0, 0, greenwichDateDay, greenwichDateMonth, greenwichDateYear);
var raHours = PAMacros.DecimalDegreesToDegreeHours(raDeg);
var decDeg = PAMacros.EcDec(lSDeg2, 0, 0, 0, 0, 0, greenwichDateDay, greenwichDateMonth, greenwichDateYear);
var sunRAHour = PAMacros.DecimalHoursHour(raHours);
var sunRAMin = PAMacros.DecimalHoursMinute(raHours);
var sunRASec = PAMacros.DecimalHoursSecond(raHours);
var sunDecDeg = PAMacros.DecimalDegreesDegrees(decDeg);
var sunDecMin = PAMacros.DecimalDegreesMinutes(decDeg);
var sunDecSec = PAMacros.DecimalDegreesSeconds(decDeg);
return(sunRAHour, sunRAMin, sunRASec, sunDecDeg, sunDecMin, sunDecSec);
}
/// <summary>
/// Calculate approximate position of a planet.
/// </summary>
public (double planetRAHour, double planetRAMin, double planetRASec, double planetDecDeg, double planetDecMin, double planetDecSec) ApproximatePositionOfPlanet(double lctHour, double lctMin, double lctSec, bool isDaylightSaving, int zoneCorrectionHours, double localDateDay, int localDateMonth, int localDateYear, string planetName)
{
var daylightSaving = (isDaylightSaving) ? 1 : 0;
var planetInfo = PlanetInfo.GetPlanetInfo(planetName);
var gdateDay = PAMacros.LocalCivilTimeGreenwichDay(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var gdateMonth = PAMacros.LocalCivilTimeGreenwichMonth(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var gdateYear = PAMacros.LocalCivilTimeGreenwichYear(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var utHours = PAMacros.LocalCivilTimeToUniversalTime(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var dDays = PAMacros.CivilDateToJulianDate(gdateDay + (utHours / 24), gdateMonth, gdateYear) - PAMacros.CivilDateToJulianDate(0, 1, 2010);
var npDeg1 = 360 * dDays / (365.242191 * planetInfo.tp_PeriodOrbit);
var npDeg2 = npDeg1 - 360 * (npDeg1 / 360).Floor();
var mpDeg = npDeg2 + planetInfo.long_LongitudeEpoch - planetInfo.peri_LongitudePerihelion;
var lpDeg1 = npDeg2 + (360 * planetInfo.ecc_EccentricityOrbit * mpDeg.ToRadians().Sine() / Math.PI) + planetInfo.long_LongitudeEpoch;
var lpDeg2 = lpDeg1 - 360 * (lpDeg1 / 360).Floor();
var planetTrueAnomalyDeg = lpDeg2 - planetInfo.peri_LongitudePerihelion;
var rAU = planetInfo.axis_AxisOrbit * (1 - Math.Pow(planetInfo.ecc_EccentricityOrbit, 2)) / (1 + planetInfo.ecc_EccentricityOrbit * planetTrueAnomalyDeg.ToRadians().Cosine());
var earthInfo = PlanetInfo.GetPlanetInfo("Earth");
var neDeg1 = 360 * dDays / (365.242191 * earthInfo.tp_PeriodOrbit);
var neDeg2 = neDeg1 - 360 * (neDeg1 / 360).Floor();
var meDeg = neDeg2 + earthInfo.long_LongitudeEpoch - earthInfo.peri_LongitudePerihelion;
var leDeg1 = neDeg2 + earthInfo.long_LongitudeEpoch + 360 * earthInfo.ecc_EccentricityOrbit * meDeg.ToRadians().Sine() / Math.PI;
var leDeg2 = leDeg1 - 360 * (leDeg1 / 360).Floor();
var earthTrueAnomalyDeg = leDeg2 - earthInfo.peri_LongitudePerihelion;
var rAU2 = earthInfo.axis_AxisOrbit * (1 - Math.Pow(earthInfo.ecc_EccentricityOrbit, 2)) / (1 + earthInfo.ecc_EccentricityOrbit * earthTrueAnomalyDeg.ToRadians().Cosine());
var lpNodeRad = (lpDeg2 - planetInfo.node_LongitudeAscendingNode).ToRadians();
var psiRad = ((lpNodeRad).Sine() * planetInfo.incl_OrbitalInclination.ToRadians().Sine()).ASine();
var y = lpNodeRad.Sine() * planetInfo.incl_OrbitalInclination.ToRadians().Cosine();
var x = lpNodeRad.Cosine();
var ldDeg = PAMacros.Degrees(y.AngleTangent2(x)) + planetInfo.node_LongitudeAscendingNode;
var rdAU = rAU * psiRad.Cosine();
var leLdRad = (leDeg2 - ldDeg).ToRadians();
var atan2Type1 = (rdAU * leLdRad.Sine()).AngleTangent2(rAU2 - rdAU * leLdRad.Cosine());
var atan2Type2 = (rAU2 * (-leLdRad).Sine()).AngleTangent2(rdAU - rAU2 * leLdRad.Cosine());
var aRad = (rdAU < 1) ? atan2Type1 : atan2Type2;
var lamdaDeg1 = (rdAU < 1) ? 180 + leDeg2 + PAMacros.Degrees(aRad) : PAMacros.Degrees(aRad) + ldDeg;
var lamdaDeg2 = lamdaDeg1 - 360 * (lamdaDeg1 / 360).Floor();
var betaDeg = PAMacros.Degrees((rdAU * psiRad.Tangent() * ((lamdaDeg2 - ldDeg).ToRadians()).Sine() / (rAU2 * (-leLdRad).Sine())).AngleTangent());
var raHours = PAMacros.DecimalDegreesToDegreeHours(PAMacros.EcRA(lamdaDeg2, 0, 0, betaDeg, 0, 0, gdateDay, gdateMonth, gdateYear));
var decDeg = PAMacros.EcDec(lamdaDeg2, 0, 0, betaDeg, 0, 0, gdateDay, gdateMonth, gdateYear);
var planetRAHour = PAMacros.DecimalHoursHour(raHours);
var planetRAMin = PAMacros.DecimalHoursMinute(raHours);
var planetRASec = PAMacros.DecimalHoursSecond(raHours);
var planetDecDeg = PAMacros.DecimalDegreesDegrees(decDeg);
var planetDecMin = PAMacros.DecimalDegreesMinutes(decDeg);
var planetDecSec = PAMacros.DecimalDegreesSeconds(decDeg);
return(planetRAHour, planetRAMin, planetRASec, planetDecDeg, planetDecMin, planetDecSec);
}
/// <summary>
/// Calculate approximate position of the Moon.
/// </summary>
/// <returns>
/// <para>moon_ra_hour -- Right ascension of Moon (hour part)</para>
/// <para>moon_ra_min -- Right ascension of Moon (minutes part)</para>
/// <para>moon_ra_sec -- Right ascension of Moon (seconds part)</para>
/// <para>moon_dec_deg -- Declination of Moon (degrees part)</para>
/// <para>moon_dec_min -- Declination of Moon (minutes part)</para>
/// <para>moon_dec_sec -- Declination of Moon (seconds part)</para>
/// </returns>
public (double moonRAHour, double moonRAMin, double moonRASec, double moonDecDeg, double moonDecMin, double moonDecSec) ApproximatePositionOfMoon(double lctHour, double lctMin, double lctSec, bool isDaylightSaving, int zoneCorrectionHours, double localDateDay, int localDateMonth, int localDateYear)
{
var daylightSaving = (isDaylightSaving) ? 1 : 0;
var l0 = 91.9293359879052;
var p0 = 130.143076320618;
var n0 = 291.682546643194;
var i = 5.145396;
var gdateDay = PAMacros.LocalCivilTimeGreenwichDay(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var gdateMonth = PAMacros.LocalCivilTimeGreenwichMonth(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var gdateYear = PAMacros.LocalCivilTimeGreenwichYear(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var utHours = PAMacros.LocalCivilTimeToUniversalTime(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var dDays = PAMacros.CivilDateToJulianDate(gdateDay, gdateMonth, gdateYear) - PAMacros.CivilDateToJulianDate(0.0, 1, 2010) + utHours / 24;
var sunLongDeg = PAMacros.SunLong(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var sunMeanAnomalyRad = PAMacros.SunMeanAnomaly(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var lmDeg = PAMacros.UnwindDeg(13.1763966 * dDays + l0);
var mmDeg = PAMacros.UnwindDeg(lmDeg - 0.1114041 * dDays - p0);
var nDeg = PAMacros.UnwindDeg(n0 - (0.0529539 * dDays));
var evDeg = 1.2739 * ((2.0 * (lmDeg - sunLongDeg) - mmDeg).ToRadians()).Sine();
var aeDeg = 0.1858 * (sunMeanAnomalyRad).Sine();
var a3Deg = 0.37 * (sunMeanAnomalyRad).Sine();
var mmdDeg = mmDeg + evDeg - aeDeg - a3Deg;
var ecDeg = 6.2886 * mmdDeg.ToRadians().Sine();
var a4Deg = 0.214 * (2.0 * (mmdDeg).ToRadians()).Sine();
var ldDeg = lmDeg + evDeg + ecDeg - aeDeg + a4Deg;
var vDeg = 0.6583 * (2.0 * (ldDeg - sunLongDeg).ToRadians()).Sine();
var lddDeg = ldDeg + vDeg;
var ndDeg = nDeg - 0.16 * (sunMeanAnomalyRad).Sine();
var y = ((lddDeg - ndDeg).ToRadians()).Sine() * i.ToRadians().Cosine();
var x = (lddDeg - ndDeg).ToRadians().Cosine();
var moonLongDeg = PAMacros.UnwindDeg(PAMacros.Degrees(y.AngleTangent2(x)) + ndDeg);
var moonLatDeg = PAMacros.Degrees(((lddDeg - ndDeg).ToRadians().Sine() * i.ToRadians().Sine()).ASine());
var moonRAHours1 = PAMacros.DecimalDegreesToDegreeHours(PAMacros.EcRA(moonLongDeg, 0, 0, moonLatDeg, 0, 0, gdateDay, gdateMonth, gdateYear));
var moonDecDeg1 = PAMacros.EcDec(moonLongDeg, 0, 0, moonLatDeg, 0, 0, gdateDay, gdateMonth, gdateYear);
var moonRAHour = PAMacros.DecimalHoursHour(moonRAHours1);
var moonRAMin = PAMacros.DecimalHoursMinute(moonRAHours1);
var moonRASec = PAMacros.DecimalHoursSecond(moonRAHours1);
var moonDecDeg = PAMacros.DecimalDegreesDegrees(moonDecDeg1);
var moonDecMin = PAMacros.DecimalDegreesMinutes(moonDecDeg1);
var moonDecSec = PAMacros.DecimalDegreesSeconds(moonDecDeg1);
return(moonRAHour, moonRAMin, moonRASec, moonDecDeg, moonDecMin, moonDecSec);
}
