/// <summary>
/// Determine if a lunar eclipse is likely to occur.
/// </summary>
/// <returns>
/// <para>status -- One of "Lunar eclipse certain", "Lunar eclipse possible", or "No lunar eclipse".</para>
/// <para>eventDateDay -- Date of eclipse event (day).</para>
/// <para>eventDateMonth -- Date of eclipse event (month).</para>
/// <para>eventDateYear -- Date of eclipse event (year).</para>
/// </returns>
public (string status, double eventDateDay, int eventDateMonth, int eventDateYear) LunarEclipseOccurrence(double localDateDay, int localDateMonth, int localDateYear, bool isDaylightSaving, int zoneCorrectionHours)
{
var daylightSaving = (isDaylightSaving) ? 1 : 0;
var julianDateOfFullMoon = PAMacros.FullMoon(daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var gDateOfFullMoonDay = PAMacros.JulianDateDay(julianDateOfFullMoon);
var integerDay = (gDateOfFullMoonDay).Floor();
var gDateOfFullMoonMonth = PAMacros.JulianDateMonth(julianDateOfFullMoon);
var gDateOfFullMoonYear = PAMacros.JulianDateYear(julianDateOfFullMoon);
var utOfFullMoonHours = gDateOfFullMoonDay - integerDay;
var localCivilDateDay = PAMacros.UniversalTime_LocalCivilDay(utOfFullMoonHours, 0.0, 0.0, daylightSaving, zoneCorrectionHours, integerDay, gDateOfFullMoonMonth, gDateOfFullMoonYear);
var localCivilDateMonth = PAMacros.UniversalTime_LocalCivilMonth(utOfFullMoonHours, 0.0, 0.0, daylightSaving, zoneCorrectionHours, integerDay, gDateOfFullMoonMonth, gDateOfFullMoonYear);
var localCivilDateYear = PAMacros.UniversalTime_LocalCivilYear(utOfFullMoonHours, 0.0, 0.0, daylightSaving, zoneCorrectionHours, integerDay, gDateOfFullMoonMonth, gDateOfFullMoonYear);
var eclipseOccurrence = PAMacros.LunarEclipseOccurrence(daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var status = eclipseOccurrence;
var eventDateDay = localCivilDateDay;
var eventDateMonth = localCivilDateMonth;
var eventDateYear = localCivilDateYear;
return(status, eventDateDay, eventDateMonth, eventDateYear);
}
/// <summary>
/// Calculate new moon and full moon instances.
/// </summary>
/// <returns>
/// <para>nmLocalTimeHour -- new Moon instant - local time (hour)</para>
/// <para>nmLocalTimeMin -- new Moon instant - local time (minutes)</para>
/// <para>nmLocalDateDay -- new Moon instance - local date (day)</para>
/// <para>nmLocalDateMonth -- new Moon instance - local date (month)</para>
/// <para>nmLocalDateYear -- new Moon instance - local date (year)</para>
/// <para>fmLocalTimeHour -- full Moon instant - local time (hour)</para>
/// <para>fmLocalTimeMin -- full Moon instant - local time (minutes)</para>
/// <para>fmLocalDateDay -- full Moon instance - local date (day)</para>
/// <para>fmLocalDateMonth -- full Moon instance - local date (month)</para>
/// <para>fmLocalDateYear -- full Moon instance - local date (year)</para>
/// </returns>
public (double nmLocalTimeHour, double nmLocalTimeMin, double nmLocalDateDay, int nmLocalDateMonth, int nmLocalDateYear, double fmLocalTimeHour, double fmLocalTimeMin, double fmLocalDateDay, int fmLocalDateMonth, int fmLocalDateYear) TimesOfNewMoonAndFullMoon(bool isDaylightSaving, int zoneCorrectionHours, double localDateDay, int localDateMonth, int localDateYear)
{
var daylightSaving = (isDaylightSaving) ? 1 : 0;
var jdOfNewMoonDays = PAMacros.NewMoon(daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var jdOfFullMoonDays = PAMacros.FullMoon(3, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var gDateOfNewMoonDay = PAMacros.JulianDateDay(jdOfNewMoonDays);
var integerDay1 = gDateOfNewMoonDay.Floor();
var gDateOfNewMoonMonth = PAMacros.JulianDateMonth(jdOfNewMoonDays);
var gDateOfNewMoonYear = PAMacros.JulianDateYear(jdOfNewMoonDays);
var gDateOfFullMoonDay = PAMacros.JulianDateDay(jdOfFullMoonDays);
var integerDay2 = gDateOfFullMoonDay.Floor();
var gDateOfFullMoonMonth = PAMacros.JulianDateMonth(jdOfFullMoonDays);
var gDateOfFullMoonYear = PAMacros.JulianDateYear(jdOfFullMoonDays);
var utOfNewMoonHours = 24.0 * (gDateOfNewMoonDay - integerDay1);
var utOfFullMoonHours = 24.0 * (gDateOfFullMoonDay - integerDay2);
var lctOfNewMoonHours = PAMacros.UniversalTimeToLocalCivilTime(utOfNewMoonHours + 0.008333, 0, 0, daylightSaving, zoneCorrectionHours, integerDay1, gDateOfNewMoonMonth, gDateOfNewMoonYear);
var lctOfFullMoonHours = PAMacros.UniversalTimeToLocalCivilTime(utOfFullMoonHours + 0.008333, 0, 0, daylightSaving, zoneCorrectionHours, integerDay2, gDateOfFullMoonMonth, gDateOfFullMoonYear);
var nmLocalTimeHour = PAMacros.DecimalHoursHour(lctOfNewMoonHours);
var nmLocalTimeMin = PAMacros.DecimalHoursMinute(lctOfNewMoonHours);
var nmLocalDateDay = PAMacros.UniversalTime_LocalCivilDay(utOfNewMoonHours, 0, 0, daylightSaving, zoneCorrectionHours, integerDay1, gDateOfNewMoonMonth, gDateOfNewMoonYear);
var nmLocalDateMonth = PAMacros.UniversalTime_LocalCivilMonth(utOfNewMoonHours, 0, 0, daylightSaving, zoneCorrectionHours, integerDay1, gDateOfNewMoonMonth, gDateOfNewMoonYear);
var nmLocalDateYear = PAMacros.UniversalTime_LocalCivilYear(utOfNewMoonHours, 0, 0, daylightSaving, zoneCorrectionHours, integerDay1, gDateOfNewMoonMonth, gDateOfNewMoonYear);
var fmLocalTimeHour = PAMacros.DecimalHoursHour(lctOfFullMoonHours);
var fmLocalTimeMin = PAMacros.DecimalHoursMinute(lctOfFullMoonHours);
var fmLocalDateDay = PAMacros.UniversalTime_LocalCivilDay(utOfFullMoonHours, 0, 0, daylightSaving, zoneCorrectionHours, integerDay2, gDateOfFullMoonMonth, gDateOfFullMoonYear);
var fmLocalDateMonth = PAMacros.UniversalTime_LocalCivilMonth(utOfFullMoonHours, 0, 0, daylightSaving, zoneCorrectionHours, integerDay2, gDateOfFullMoonMonth, gDateOfFullMoonYear);
var fmLocalDateYear = PAMacros.UniversalTime_LocalCivilYear(utOfFullMoonHours, 0, 0, daylightSaving, zoneCorrectionHours, integerDay2, gDateOfFullMoonMonth, gDateOfFullMoonYear);
return(nmLocalTimeHour, nmLocalTimeMin, nmLocalDateDay, nmLocalDateMonth, nmLocalDateYear, fmLocalTimeHour, fmLocalTimeMin, fmLocalDateDay, fmLocalDateMonth, fmLocalDateYear);
}
/// <summary>
/// Calculate the circumstances of a lunar eclipse.
/// </summary>
/// <returns>
/// <para>lunarEclipseCertainDateDay -- Lunar eclipse date (day)</para>
/// <para>lunarEclipseCertainDateMonth -- Lunar eclipse date (month)</para>
/// <para>lunarEclipseCertainDateYear -- Lunar eclipse date (year)</para>
/// <para>utstartPenPhaseHour -- Start of penumbral phase (hour)</para>
/// <para>utStartPenPhaseMinutes -- Start of penumbral phase (minutes)</para>
/// <para>utStartUmbralPhaseHour -- Start of umbral phase (hour)</para>
/// <para>utStartUmbralPhaseMinutes -- Start of umbral phase (minutes)</para>
/// <para>utStartTotalPhaseHour -- Start of total phase (hour)</para>
/// <para>utStartTotalPhaseMinutes -- Start of total phase (minutes)</para>
/// <para>utMidEclipseHour -- Mid-eclipse (hour)</para>
/// <para>utMidEclipseMinutes -- Mid-eclipse (minutes)</para>
/// <para>utEndTotalPhaseHour -- End of total phase (hour)</para>
/// <para>utEndTotalPhaseMinutes -- End of total phase (minutes)</para>
/// <para>utEndUmbralPhaseHour -- End of umbral phase (hour)</para>
/// <para>utEndUmbralPhaseMinutes -- End of umbral phase (minutes)</para>
/// <para>utEndPenPhaseHour -- End of penumbral phase (hour)</para>
/// <para>utEndPenPhaseMinutes -- End of penumbral phase (minutes)</para>
/// <para>eclipseMagnitude -- Eclipse magnitude</para>
/// </returns>
public (double lunarEclipseCertainDateDay, double lunarEclipseCertainDateMonth, double lunarEclipseCertainDateYear, double utStartPenPhaseHour, double utStartPenPhaseMinutes, double utStartUmbralPhaseHour, double utStartUmbralPhaseMinutes, double utStartTotalPhaseHour, double utStartTotalPhaseMinutes, double utMidEclipseHour, double utMidEclipseMinutes, double utEndTotalPhaseHour, double utEndTotalPhaseMinutes, double utEndUmbralPhaseHour, double utEndUmbralPhaseMinutes, double utEndPenPhaseHour, double utEndPenPhaseMinutes, double eclipseMagnitude) LunarEclipseCircumstances(double localDateDay, int localDateMonth, int localDateYear, bool isDaylightSaving, int zoneCorrectionHours)
{
var daylightSaving = (isDaylightSaving) ? 1 : 0;
var julianDateOfFullMoon = PAMacros.FullMoon(daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var gDateOfFullMoonDay = PAMacros.JulianDateDay(julianDateOfFullMoon);
var integerDay = gDateOfFullMoonDay.Floor();
var gDateOfFullMoonMonth = PAMacros.JulianDateMonth(julianDateOfFullMoon);
var gDateOfFullMoonYear = PAMacros.JulianDateYear(julianDateOfFullMoon);
var utOfFullMoonHours = gDateOfFullMoonDay - integerDay;
var localCivilDateDay = PAMacros.UniversalTime_LocalCivilDay(utOfFullMoonHours, 0.0, 0.0, daylightSaving, zoneCorrectionHours, integerDay, gDateOfFullMoonMonth, gDateOfFullMoonYear);
var localCivilDateMonth = PAMacros.UniversalTime_LocalCivilMonth(utOfFullMoonHours, 0.0, 0.0, daylightSaving, zoneCorrectionHours, integerDay, gDateOfFullMoonMonth, gDateOfFullMoonYear);
var localCivilDateYear = PAMacros.UniversalTime_LocalCivilYear(utOfFullMoonHours, 0.0, 0.0, daylightSaving, zoneCorrectionHours, integerDay, gDateOfFullMoonMonth, gDateOfFullMoonYear);
var utMaxEclipse = PAMacros.UTMaxLunarEclipse(localDateDay, localDateMonth, localDateYear, daylightSaving, zoneCorrectionHours);
var utFirstContact = PAMacros.UTFirstContactLunarEclipse(localDateDay, localDateMonth, localDateYear, daylightSaving, zoneCorrectionHours);
var utLastContact = PAMacros.UTLastContactLunarEclipse(localDateDay, localDateMonth, localDateYear, daylightSaving, zoneCorrectionHours);
var utStartUmbralPhase = PAMacros.UTStartUmbraLunarEclipse(localDateDay, localDateMonth, localDateYear, daylightSaving, zoneCorrectionHours);
var utEndUmbralPhase = PAMacros.UTEndUmbraLunarEclipse(localDateDay, localDateMonth, localDateYear, daylightSaving, zoneCorrectionHours);
var utStartTotalPhase = PAMacros.UTStartTotalLunarEclipse(localDateDay, localDateMonth, localDateYear, daylightSaving, zoneCorrectionHours);
var utEndTotalPhase = PAMacros.UTEndTotalLunarEclipse(localDateDay, localDateMonth, localDateYear, daylightSaving, zoneCorrectionHours);
var eclipseMagnitude1 = PAMacros.MagLunarEclipse(localDateDay, localDateMonth, localDateYear, daylightSaving, zoneCorrectionHours);
var lunarEclipseCertainDateDay = localCivilDateDay;
var lunarEclipseCertainDateMonth = localCivilDateMonth;
var lunarEclipseCertainDateYear = localCivilDateYear;
var utStartPenPhaseHour = (utFirstContact == -99.0) ? -99.0 : PAMacros.DecimalHoursHour(utFirstContact + 0.008333);
var utStartPenPhaseMinutes = (utFirstContact == -99.0) ? -99.0 : PAMacros.DecimalHoursMinute(utFirstContact + 0.008333);
var utStartUmbralPhaseHour = (utStartUmbralPhase == -99.0) ? -99.0 : PAMacros.DecimalHoursHour(utStartUmbralPhase + 0.008333);
var utStartUmbralPhaseMinutes = (utStartUmbralPhase == -99.0) ? -99.0 : PAMacros.DecimalHoursMinute(utStartUmbralPhase + 0.008333);
var utStartTotalPhaseHour = (utStartTotalPhase == -99.0) ? -99.0 : PAMacros.DecimalHoursHour(utStartTotalPhase + 0.008333);
var utStartTotalPhaseMinutes = (utStartTotalPhase == -99.0) ? -99.0 : PAMacros.DecimalHoursMinute(utStartTotalPhase + 0.008333);
var utMidEclipseHour = (utMaxEclipse == -99.0) ? -99.0 : PAMacros.DecimalHoursHour(utMaxEclipse + 0.008333);
var utMidEclipseMinutes = (utMaxEclipse == -99.0) ? -99.0 : PAMacros.DecimalHoursMinute(utMaxEclipse + 0.008333);
var utEndTotalPhaseHour = (utEndTotalPhase == -99.0) ? -99.0 : PAMacros.DecimalHoursHour(utEndTotalPhase + 0.008333);
var utEndTotalPhaseMinutes = (utEndTotalPhase == -99.0) ? -99.0 : PAMacros.DecimalHoursMinute(utEndTotalPhase + 0.008333);
var utEndUmbralPhaseHour = (utEndUmbralPhase == -99.0) ? -99.0 : PAMacros.DecimalHoursHour(utEndUmbralPhase + 0.008333);
var utEndUmbralPhaseMinutes = (utEndUmbralPhase == -99.0) ? -99.0 : PAMacros.DecimalHoursMinute(utEndUmbralPhase + 0.008333);
var utEndPenPhaseHour = (utLastContact == -99.0) ? -99.0 : PAMacros.DecimalHoursHour(utLastContact + 0.008333);
var utEndPenPhaseMinutes = (utLastContact == -99.0) ? -99.0 : PAMacros.DecimalHoursMinute(utLastContact + 0.008333);
var eclipseMagnitude = (eclipseMagnitude1 == -99.0) ? -99.0 : Math.Round(eclipseMagnitude1, 2);
return(lunarEclipseCertainDateDay, lunarEclipseCertainDateMonth, lunarEclipseCertainDateYear, utStartPenPhaseHour, utStartPenPhaseMinutes, utStartUmbralPhaseHour, utStartUmbralPhaseMinutes, utStartTotalPhaseHour, utStartTotalPhaseMinutes, utMidEclipseHour, utMidEclipseMinutes, utEndTotalPhaseHour, utEndTotalPhaseMinutes, utEndUmbralPhaseHour, utEndUmbralPhaseMinutes, utEndPenPhaseHour, utEndPenPhaseMinutes, eclipseMagnitude);
}