/// <summary>
/// Calculate precise position of the Moon.
/// </summary>
/// <returns>
/// <para>moonRAHour -- Right ascension of Moon (hour part)</para>
/// <para>moonRAMin -- Right ascension of Moon (minutes part)</para>
/// <para>moonRASec -- Right ascension of Moon (seconds part)</para>
/// <para>moonDecDeg -- Declination of Moon (degrees part)</para>
/// <para>moonDecMin -- Declination of Moon (minutes part)</para>
/// <para>moonDecSec -- Declination of Moon (seconds part)</para>
/// <para>earthMoonDistKM -- Distance from Earth to Moon (km)</para>
/// <para>moonHorParallaxDeg -- Horizontal parallax of Moon (degrees)</para>
/// </returns>
public (double moonRAHour, double moonRAMin, double moonRASec, double moonDecDeg, double moonDecMin, double moonDecSec, double earthMoonDistKM, double moonHorParallaxDeg) PrecisePositionOfMoon(double lctHour, double lctMin, double lctSec, bool isDaylightSaving, int zoneCorrectionHours, double localDateDay, int localDateMonth, int localDateYear)
{
var daylightSaving = (isDaylightSaving) ? 1 : 0;
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 moonResult = PAMacros.MoonLongLatHP(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var nutationInLongitudeDeg = PAMacros.NutatLong(gdateDay, gdateMonth, gdateYear);
var correctedLongDeg = moonResult.moonLongDeg + nutationInLongitudeDeg;
var earthMoonDistanceKM = 6378.14 / moonResult.moonHorPara.ToRadians().Sine();
var moonRAHours1 = PAMacros.DecimalDegreesToDegreeHours(PAMacros.EcRA(correctedLongDeg, 0, 0, moonResult.moonLatDeg, 0, 0, gdateDay, gdateMonth, gdateYear));
var moonDecDeg1 = PAMacros.EcDec(correctedLongDeg, 0, 0, moonResult.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);
var earthMoonDistKM = Math.Round(earthMoonDistanceKM, 0);
var moonHorParallaxDeg = Math.Round(moonResult.moonHorPara, 6);
return(moonRAHour, moonRAMin, moonRASec, moonDecDeg, moonDecMin, moonDecSec, earthMoonDistKM, moonHorParallaxDeg);
}
/// <summary>
/// Convert Equatorial Coordinates to Ecliptic Coordinates
/// </summary>
/// <returns>Tuple (outEclLongDeg, outEclLongMin, outEclLongSec, outEclLatDeg, outEclLatMin, outEclLatSec)</returns>
public (double outEclLongDeg, double outEclLongMin, double outEclLongSec, double outEclLatDeg, double outEclLatMin, double outEclLatSec) EquatorialCoordinateToEclipticCoordinate(double raHours, double raMinutes, double raSeconds, double decDegrees, double decMinutes, double decSeconds, double gwDay, int gwMonth, int gwYear)
{
var raDeg = PAMacros.DegreeHoursToDecimalDegrees(PAMacros.HMStoDH(raHours, raMinutes, raSeconds));
var decDeg = PAMacros.DegreesMinutesSecondsToDecimalDegrees(decDegrees, decMinutes, decSeconds);
var raRad = raDeg.ToRadians();
var decRad = decDeg.ToRadians();
var obliqDeg = PAMacros.Obliq(gwDay, gwMonth, gwYear);
var obliqRad = obliqDeg.ToRadians();
var sinEclLat = decRad.Sine() * obliqRad.Cosine() - decRad.Cosine() * obliqRad.Sine() * raRad.Sine();
var eclLatRad = sinEclLat.ASine();
var eclLatDeg = PAMacros.Degrees(eclLatRad);
var y = raRad.Sine() * obliqRad.Cosine() + decRad.Tangent() * obliqRad.Sine();
var x = raRad.Cosine();
var eclLongRad = y.AngleTangent2(x);
var eclLongDeg1 = PAMacros.Degrees(eclLongRad);
var eclLongDeg2 = eclLongDeg1 - 360 * (eclLongDeg1 / 360).Floor();
var outEclLongDeg = PAMacros.DecimalDegreesDegrees(eclLongDeg2);
var outEclLongMin = PAMacros.DecimalDegreesMinutes(eclLongDeg2);
var outEclLongSec = PAMacros.DecimalDegreesSeconds(eclLongDeg2);
var outEclLatDeg = PAMacros.DecimalDegreesDegrees(eclLatDeg);
var outEclLatMin = PAMacros.DecimalDegreesMinutes(eclLatDeg);
var outEclLatSec = PAMacros.DecimalDegreesSeconds(eclLatDeg);
return(outEclLongDeg, outEclLongMin, outEclLongSec, outEclLatDeg, outEclLatMin, outEclLatSec);
}
/// <summary>
/// Convert Ecliptic Coordinates to Equatorial Coordinates
/// </summary>
/// <returns>Tuple (outRAHours, outRAMinutes, outRASeconds, outDecDegrees, outDecMinutes, outDecSeconds)</returns>
public (double outRAHours, double outRAMinutes, double outRASeconds, double outDecDegrees, double outDecMinutes, double outDecSeconds) EclipticCoordinateToEquatorialCoordinate(double eclipticLongitudeDegrees, double eclipticLongitudeMinutes, double eclipticLongitudeSeconds, double eclipticLatitudeDegrees, double eclipticLatitudeMinutes, double eclipticLatitudeSeconds, double greenwichDay, int greenwichMonth, int greenwichYear)
{
var eclonDeg = PAMacros.DegreesMinutesSecondsToDecimalDegrees(eclipticLongitudeDegrees, eclipticLongitudeMinutes, eclipticLongitudeSeconds);
var eclatDeg = PAMacros.DegreesMinutesSecondsToDecimalDegrees(eclipticLatitudeDegrees, eclipticLatitudeMinutes, eclipticLatitudeSeconds);
var eclonRad = eclonDeg.ToRadians();
var eclatRad = eclatDeg.ToRadians();
var obliqDeg = PAMacros.Obliq(greenwichDay, greenwichMonth, greenwichYear);
var obliqRad = obliqDeg.ToRadians();
var sinDec = eclatRad.Sine() * obliqRad.Cosine() + eclatRad.Cosine() * obliqRad.Sine() * eclonRad.Sine();
var decRad = sinDec.ASine();
var decDeg = PAMacros.Degrees(decRad);
var y = eclonRad.Sine() * obliqRad.Cosine() - eclatRad.Tangent() * obliqRad.Sine();
var x = eclonRad.Cosine();
var raRad = y.AngleTangent2(x);
var raDeg1 = PAMacros.Degrees(raRad);
var raDeg2 = raDeg1 - 360 * (raDeg1 / 360).Floor();
var raHours = PAMacros.DecimalDegreesToDegreeHours(raDeg2);
var outRAHours = PAMacros.DecimalHoursHour(raHours);
var outRAMinutes = PAMacros.DecimalHoursMinute(raHours);
var outRASeconds = PAMacros.DecimalHoursSecond(raHours);
var outDecDegrees = PAMacros.DecimalDegreesDegrees(decDeg);
var outDecMinutes = PAMacros.DecimalDegreesMinutes(decDeg);
var outDecSeconds = PAMacros.DecimalDegreesSeconds(decDeg);
return(outRAHours, outRAMinutes, outRASeconds, outDecDegrees, outDecMinutes, outDecSeconds);
}
/// <summary>
/// Convert Galactic Coordinates to Equatorial Coordinates
/// </summary>
/// <returns>Tuple (raHours, raMinutes, raSeconds, decDegrees, decMinutes, decSeconds)</returns>
public (double raHours, double raMinutes, double raSeconds, double decDegrees, double decMinutes, double decSeconds) GalacticCoordinateToEquatorialCoordinate(double galLongDeg, double galLongMin, double galLongSec, double galLatDeg, double galLatMin, double galLatSec)
{
var glongDeg = PAMacros.DegreesMinutesSecondsToDecimalDegrees(galLongDeg, galLongMin, galLongSec);
var glatDeg = PAMacros.DegreesMinutesSecondsToDecimalDegrees(galLatDeg, galLatMin, galLatSec);
var glongRad = glongDeg.ToRadians();
var glatRad = glatDeg.ToRadians();
var sinDec = glatRad.Cosine() * (27.4).ToRadians().Cosine() * (glongRad - (33.0).ToRadians()).Sine() + glatRad.Sine() * (27.4).ToRadians().Sine();
var decRadians = sinDec.ASine();
var decDeg = PAMacros.Degrees(decRadians);
var y = glatRad.Cosine() * (glongRad - (33.0).ToRadians()).Cosine();
var x = glatRad.Sine() * ((27.4).ToRadians()).Cosine() - (glatRad).Cosine() * ((27.4).ToRadians()).Sine() * (glongRad - (33.0).ToRadians()).Sine();
var raDeg1 = PAMacros.Degrees(y.AngleTangent2(x)) + 192.25;
var raDeg2 = raDeg1 - 360 * (raDeg1 / 360).Floor();
var raHours1 = PAMacros.DecimalDegreesToDegreeHours(raDeg2);
var raHours = PAMacros.DecimalHoursHour(raHours1);
var raMinutes = PAMacros.DecimalHoursMinute(raHours1);
var raSeconds = PAMacros.DecimalHoursSecond(raHours1);
var decDegrees = PAMacros.DecimalDegreesDegrees(decDeg);
var decMinutes = PAMacros.DecimalDegreesMinutes(decDeg);
var decSeconds = PAMacros.DecimalDegreesSeconds(decDeg);
return(raHours, raMinutes, raSeconds, decDegrees, decMinutes, decSeconds);
}
/// <summary>
/// Calculate the angle between two celestial objects
/// </summary>
/// <returns>Tuple (angleDeg, angleMin, angleSec)</returns>
public (double angleDeg, double angleMin, double angleSec) AngleBetweenTwoObjects(double raLong1HourDeg, double raLong1Min, double raLong1Sec, double decLat1Deg, double decLat1Min, double decLat1Sec, double raLong2HourDeg, double raLong2Min, double raLong2Sec, double decLat2Deg, double decLat2Min, double decLat2Sec, PAAngleMeasure hourOrDegree)
{
var raLong1Decimal = (hourOrDegree == PAAngleMeasure.Hours) ? PAMacros.HMStoDH(raLong1HourDeg, raLong1Min, raLong1Sec) : PAMacros.DegreesMinutesSecondsToDecimalDegrees(raLong1HourDeg, raLong1Min, raLong1Sec);
var raLong1Deg = (hourOrDegree == PAAngleMeasure.Hours) ? PAMacros.DegreeHoursToDecimalDegrees(raLong1Decimal) : raLong1Decimal;
var raLong1Rad = raLong1Deg.ToRadians();
var decLat1Deg1 = PAMacros.DegreesMinutesSecondsToDecimalDegrees(decLat1Deg, decLat1Min, decLat1Sec);
var decLat1Rad = decLat1Deg1.ToRadians();
var raLong2Decimal = (hourOrDegree == PAAngleMeasure.Hours) ? PAMacros.HMStoDH(raLong2HourDeg, raLong2Min, raLong2Sec) : PAMacros.DegreesMinutesSecondsToDecimalDegrees(raLong2HourDeg, raLong2Min, raLong2Sec);
var raLong2Deg = (hourOrDegree == PAAngleMeasure.Hours) ? PAMacros.DegreeHoursToDecimalDegrees(raLong2Decimal) : raLong2Decimal;
var raLong2Rad = raLong2Deg.ToRadians();
var decLat2Deg1 = PAMacros.DegreesMinutesSecondsToDecimalDegrees(decLat2Deg, decLat2Min, decLat2Sec);
var decLat2Rad = decLat2Deg1.ToRadians();
var cosD = decLat1Rad.Sine() * decLat2Rad.Sine() + decLat1Rad.Cosine() * decLat2Rad.Cosine() * (raLong1Rad - raLong2Rad).Cosine();
var dRad = cosD.ACosine();
var dDeg = PAMacros.Degrees(dRad);
var angleDeg = PAMacros.DecimalDegreesDegrees(dDeg);
var angleMin = PAMacros.DecimalDegreesMinutes(dDeg);
var angleSec = PAMacros.DecimalDegreesSeconds(dDeg);
return(angleDeg, angleMin, angleSec);
}
/// <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 position of a parabolic comet.
/// </summary>
/// <returns>
/// cometRAHour -- Right ascension of comet (hour part)
/// cometRAMin -- Right ascension of comet (minutes part)
/// cometRASec -- Right ascension of comet (seconds part)
/// cometDecDeg -- Declination of comet (degrees part)
/// cometDecMin -- Declination of comet (minutes part)
/// cometDecSec -- Declination of comet (seconds part)
/// cometDistEarth -- Comet's distance from Earth (AU)
/// </returns>
/// <returns></returns>
public (double cometRAHour, double cometRAMin, double cometRASec, double cometDecDeg, double cometDecMin, double cometDecSec, double cometDistEarth) PositionOfParabolicComet(double lctHour, double lctMin, double lctSec, bool isDaylightSaving, int zoneCorrectionHours, double localDateDay, int localDateMonth, int localDateYear, string cometName)
{
var daylightSaving = (isDaylightSaving) ? 1 : 0;
var greenwichDateDay = PAMacros.LocalCivilTimeGreenwichDay(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var greenwichDateMonth = PAMacros.LocalCivilTimeGreenwichMonth(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var greenwichDateYear = PAMacros.LocalCivilTimeGreenwichYear(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var cometInfo = CometInfoParabolic.GetCometParabolicInfo(cometName);
var perihelionEpochDay = cometInfo.EpochPeriDay;
var perihelionEpochMonth = cometInfo.EpochPeriMonth;
var perihelionEpochYear = cometInfo.EpochPeriYear;
var qAU = cometInfo.PeriDist;
var inclinationDeg = cometInfo.Incl;
var perihelionDeg = cometInfo.ArgPeri;
var nodeDeg = cometInfo.Node;
var cometLongLatDist = PAMacros.PCometLongLatDist(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear, perihelionEpochDay, perihelionEpochMonth, perihelionEpochYear, qAU, inclinationDeg, perihelionDeg, nodeDeg);
var cometRAHours = PAMacros.DecimalDegreesToDegreeHours(PAMacros.EcRA(cometLongLatDist.cometLongDeg, 0, 0, cometLongLatDist.cometLatDeg, 0, 0, greenwichDateDay, greenwichDateMonth, greenwichDateYear));
var cometDecDeg1 = PAMacros.EcDec(cometLongLatDist.cometLongDeg, 0, 0, cometLongLatDist.cometLatDeg, 0, 0, greenwichDateDay, greenwichDateMonth, greenwichDateYear);
var cometRAHour = PAMacros.DecimalHoursHour(cometRAHours);
var cometRAMin = PAMacros.DecimalHoursMinute(cometRAHours);
var cometRASec = PAMacros.DecimalHoursSecond(cometRAHours);
var cometDecDeg = PAMacros.DecimalDegreesDegrees(cometDecDeg1);
var cometDecMin = PAMacros.DecimalDegreesMinutes(cometDecDeg1);
var cometDecSec = PAMacros.DecimalDegreesSeconds(cometDecDeg1);
var cometDistEarth = Math.Round(cometLongLatDist.cometDistAU, 2);
return(cometRAHour, cometRAMin, cometRASec, cometDecDeg, cometDecMin, cometDecSec, cometDistEarth);
}
/// <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>
/// Convert Horizon Coordinates to Equatorial Coordinates
/// </summary>
/// <returns>Tuple (hourAngleHours, hourAngleMinutes, hourAngleSeconds, declinationDegrees, declinationMinutes, declinationSeconds)</returns>
public (double hour_angle_hours, double hour_angle_minutes, double hour_angle_seconds, double declination_degrees, double declination_minutes, double declinationseconds) HorizonCoordinatesToEquatorialCoordinates(double azimuthDegrees, double azimuthMinutes, double azimuthSeconds, double altitudeDegrees, double altitudeMinutes, double altitudeSeconds, double geographicalLatitude)
{
var hourAngleInDecimalDegrees = PAMacros.HorizonCoordinatesToHourAngle(azimuthDegrees, azimuthMinutes, azimuthSeconds, altitudeDegrees, altitudeMinutes, altitudeSeconds, geographicalLatitude);
var declinationInDecimalDegrees = PAMacros.HorizonCoordinatesToDeclination(azimuthDegrees, azimuthMinutes, azimuthSeconds, altitudeDegrees, altitudeMinutes, altitudeSeconds, geographicalLatitude);
var hourAngleHours = PAMacros.DecimalHoursHour(hourAngleInDecimalDegrees);
var hourAngleMinutes = PAMacros.DecimalHoursMinute(hourAngleInDecimalDegrees);
var hourAngleSeconds = PAMacros.DecimalHoursSecond(hourAngleInDecimalDegrees);
var declinationDegrees = PAMacros.DecimalDegreesDegrees(declinationInDecimalDegrees);
var declinationMinutes = PAMacros.DecimalDegreesMinutes(declinationInDecimalDegrees);
var declinationSeconds = PAMacros.DecimalDegreesSeconds(declinationInDecimalDegrees);
return(hourAngleHours, hourAngleMinutes, hourAngleSeconds, declinationDegrees, declinationMinutes, declinationSeconds);
}
/// <summary>
/// Convert Equatorial Coordinates to Horizon Coordinates
/// </summary>
/// <returns>Tuple (azimuthDegrees, azimuthMinutes, azimuthSeconds, altitudeDegrees, altitudeMinutes, altitudeSeconds)</returns>
public (double azimuthDegrees, double azimuthMinutes, double azimuthSeconds, double altitudeDegrees, double altitudeMinutes, double altitudeSeconds) EquatorialCoordinatesToHorizonCoordinates(double hourAngleHours, double hourAngleMinutes, double hourAngleSeconds, double declinationDegrees, double declinationMinutes, double declinationSeconds, double geographicalLatitude)
{
var azimuthInDecimalDegrees = PAMacros.EquatorialCoordinatesToAzimuth(hourAngleHours, hourAngleMinutes, hourAngleSeconds, declinationDegrees, declinationMinutes, declinationSeconds, geographicalLatitude);
var altitudeInDecimalDegrees = PAMacros.EquatorialCoordinatesToAltitude(hourAngleHours, hourAngleMinutes, hourAngleSeconds, declinationDegrees, declinationMinutes, declinationSeconds, geographicalLatitude);
var azimuthDegrees = PAMacros.DecimalDegreesDegrees(azimuthInDecimalDegrees);
var azimuthMinutes = PAMacros.DecimalDegreesMinutes(azimuthInDecimalDegrees);
var azimuthSeconds = PAMacros.DecimalDegreesSeconds(azimuthInDecimalDegrees);
var altitudeDegrees = PAMacros.DecimalDegreesDegrees(altitudeInDecimalDegrees);
var altitudeMinutes = PAMacros.DecimalDegreesMinutes(altitudeInDecimalDegrees);
var altitudeSeconds = PAMacros.DecimalDegreesSeconds(altitudeInDecimalDegrees);
return(azimuthDegrees, azimuthMinutes, azimuthSeconds, altitudeDegrees, altitudeMinutes, altitudeSeconds);
}
/// <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);
}
/// <summary>
/// Calculate position of an elliptical comet.
/// </summary>
/// <returns>
/// cometRAHour -- Right ascension of comet (hour part)
/// cometRAMin -- Right ascension of comet (minutes part)
/// cometDecDeg -- Declination of comet (degrees part)
/// cometDecMin -- Declination of comet (minutes part)
/// cometDistEarth -- Comet's distance from Earth (AU)
/// </returns>
public (double cometRAHour, double cometRAMin, double cometDecDeg, double cometDecMin, double cometDistEarth) PositionOfEllipticalComet(double lctHour, double lctMin, double lctSec, bool isDaylightSaving, int zoneCorrectionHours, double localDateDay, int localDateMonth, int localDateYear, string cometName)
{
var daylightSaving = (isDaylightSaving) ? 1 : 0;
var greenwichDateDay = PAMacros.LocalCivilTimeGreenwichDay(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var greenwichDateMonth = PAMacros.LocalCivilTimeGreenwichMonth(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var greenwichDateYear = PAMacros.LocalCivilTimeGreenwichYear(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var cometInfo = CometInfoElliptical.GetCometEllipticalInfo(cometName);
var timeSinceEpochYears = (PAMacros.CivilDateToJulianDate(greenwichDateDay, greenwichDateMonth, greenwichDateYear) - PAMacros.CivilDateToJulianDate(0.0, 1, greenwichDateYear)) / 365.242191 + greenwichDateYear - cometInfo.epoch_EpochOfPerihelion;
var mcDeg = 360 * timeSinceEpochYears / cometInfo.period_PeriodOfOrbit;
var mcRad = (mcDeg - 360 * (mcDeg / 360).Floor()).ToRadians();
var eccentricity = cometInfo.ecc_EccentricityOfOrbit;
var trueAnomalyDeg = PAMacros.Degrees(PAMacros.TrueAnomaly(mcRad, eccentricity));
var lcDeg = trueAnomalyDeg + cometInfo.peri_LongitudeOfPerihelion;
var rAU = cometInfo.axis_SemiMajorAxisOfOrbit * (1 - eccentricity * eccentricity) / (1 + eccentricity * ((trueAnomalyDeg).ToRadians()).Cosine());
var lcNodeRad = (lcDeg - cometInfo.node_LongitudeOfAscendingNode).ToRadians();
var psiRad = ((lcNodeRad).Sine() * ((cometInfo.incl_InclinationOfOrbit).ToRadians()).Sine()).ASine();
var y = (lcNodeRad).Sine() * ((cometInfo.incl_InclinationOfOrbit).ToRadians()).Cosine();
var x = (lcNodeRad).Cosine();
var ldDeg = PAMacros.Degrees(y.AngleTangent2(x)) + cometInfo.node_LongitudeOfAscendingNode;
var rdAU = rAU * (psiRad).Cosine();
var earthLongitudeLeDeg = PAMacros.SunLong(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear) + 180.0;
var earthRadiusVectorAU = PAMacros.SunDist(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var leLdRad = (earthLongitudeLeDeg - ldDeg).ToRadians();
var aRad = (rdAU < earthRadiusVectorAU) ? (rdAU * (leLdRad).Sine()).AngleTangent2(earthRadiusVectorAU - rdAU * (leLdRad).Cosine()) : (earthRadiusVectorAU * (-leLdRad).Sine()).AngleTangent2(rdAU - earthRadiusVectorAU * (leLdRad).Cosine());
var cometLongDeg1 = (rdAU < earthRadiusVectorAU) ? 180.0 + earthLongitudeLeDeg + PAMacros.Degrees(aRad) : PAMacros.Degrees(aRad) + ldDeg;
var cometLongDeg = cometLongDeg1 - 360 * (cometLongDeg1 / 360).Floor();
var cometLatDeg = PAMacros.Degrees((rdAU * (psiRad).Tangent() * ((cometLongDeg1 - ldDeg).ToRadians()).Sine() / (earthRadiusVectorAU * (-leLdRad).Sine())).AngleTangent());
var cometRAHours1 = PAMacros.DecimalDegreesToDegreeHours(PAMacros.EcRA(cometLongDeg, 0, 0, cometLatDeg, 0, 0, greenwichDateDay, greenwichDateMonth, greenwichDateYear));
var cometDecDeg1 = PAMacros.EcDec(cometLongDeg, 0, 0, cometLatDeg, 0, 0, greenwichDateDay, greenwichDateMonth, greenwichDateYear);
var cometDistanceAU = (Math.Pow(earthRadiusVectorAU, 2) + Math.Pow(rAU, 2) - 2.0 * earthRadiusVectorAU * rAU * ((lcDeg - earthLongitudeLeDeg).ToRadians()).Cosine() * (psiRad).Cosine()).SquareRoot();
var cometRAHour = PAMacros.DecimalHoursHour(cometRAHours1 + 0.008333);
var cometRAMin = PAMacros.DecimalHoursMinute(cometRAHours1 + 0.008333);
var cometDecDeg = PAMacros.DecimalDegreesDegrees(cometDecDeg1 + 0.008333);
var cometDecMin = PAMacros.DecimalDegreesMinutes(cometDecDeg1 + 0.008333);
var cometDistEarth = Math.Round(cometDistanceAU, 2);
return(cometRAHour, cometRAMin, cometDecDeg, cometDecMin, cometDistEarth);
}
/// <summary>
/// Calculate Moon's distance, angular diameter, and horizontal parallax.
/// </summary>
/// <returns>
/// <para>earth_moon_dist -- Earth-Moon distance (km)</para>
/// <para>ang_diameter_deg -- Angular diameter (degrees part)</para>
/// <para>ang_diameter_min -- Angular diameter (minutes part)</para>
/// <para>hor_parallax_deg -- Horizontal parallax (degrees part)</para>
/// <para>hor_parallax_min -- Horizontal parallax (minutes part)</para>
/// <para>hor_parallax_sec -- Horizontal parallax (seconds part)</para>
/// </returns>
public (double earthMoonDist, double angDiameterDeg, double angDiameterMin, double horParallaxDeg, double horParallaxMin, double horParallaxSec) MoonDistAngDiamHorParallax(double lctHour, double lctMin, double lctSec, bool isDaylightSaving, int zoneCorrectionHours, double localDateDay, int localDateMonth, int localDateYear)
{
var daylightSaving = (isDaylightSaving) ? 1 : 0;
var moonDistance = PAMacros.MoonDist(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var moonAngularDiameter = PAMacros.MoonSize(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var moonHorizontalParallax = PAMacros.MoonHP(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear);
var earthMoonDist = Math.Round(moonDistance, 0);
var angDiameterDeg = PAMacros.DecimalDegreesDegrees(moonAngularDiameter + 0.008333);
var angDiameterMin = PAMacros.DecimalDegreesMinutes(moonAngularDiameter + 0.008333);
var horParallaxDeg = PAMacros.DecimalDegreesDegrees(moonHorizontalParallax);
var horParallaxMin = PAMacros.DecimalDegreesMinutes(moonHorizontalParallax);
var horParallaxSec = PAMacros.DecimalDegreesSeconds(moonHorizontalParallax);
return(earthMoonDist, angDiameterDeg, angDiameterMin, horParallaxDeg, horParallaxMin, horParallaxSec);
}
/// <summary>
/// Calculate precise position of a planet.
/// </summary>
public (double planetRAHour, double planetRAMin, double planetRASec, double planetDecDeg, double planetDecMin, double planetDecSec) PrecisePositionOfPlanet(double lctHour, double lctMin, double lctSec, bool isDaylightSaving, int zoneCorrectionHours, double localDateDay, int localDateMonth, int localDateYear, string planetName)
{
var daylightSaving = (isDaylightSaving) ? 1 : 0;
var coordinateResults = PAMacros.PlanetCoordinates(lctHour, lctMin, lctSec, daylightSaving, zoneCorrectionHours, localDateDay, localDateMonth, localDateYear, planetName);
var planetRAHours = PAMacros.DecimalDegreesToDegreeHours(PAMacros.EcRA(coordinateResults.planetLongitude, 0, 0, coordinateResults.planetLatitude, 0, 0, localDateDay, localDateMonth, localDateYear));
var planetDecDeg1 = PAMacros.EcDec(coordinateResults.planetLongitude, 0, 0, coordinateResults.planetLatitude, 0, 0, localDateDay, localDateMonth, localDateYear);
var planetRAHour = PAMacros.DecimalHoursHour(planetRAHours);
var planetRAMin = PAMacros.DecimalHoursMinute(planetRAHours);
var planetRASec = PAMacros.DecimalHoursSecond(planetRAHours);
var planetDecDeg = PAMacros.DecimalDegreesDegrees(planetDecDeg1);
var planetDecMin = PAMacros.DecimalDegreesMinutes(planetDecDeg1);
var planetDecSec = PAMacros.DecimalDegreesSeconds(planetDecDeg1);
return(planetRAHour, planetRAMin, planetRASec, planetDecDeg, planetDecMin, planetDecSec);
}
/// <summary>
/// Calculate corrected RA/Dec, accounting for geocentric parallax.
/// </summary>
/// <returns>corrected RA hours,minutes,seconds and corrected Declination degrees,minutes,seconds</returns>
public (double correctedRAHour, double correctedRAMin, double correctedRASec, double correctedDecDeg, double correctedDecMin, double correctedDecSec) CorrectionsForGeocentricParallax(double raHour, double raMin, double raSec, double decDeg, double decMin, double decSec, PACoordinateType coordinateType, double equatorialHorParallaxDeg, double geogLongDeg, double geogLatDeg, double heightM, int daylightSaving, int timezoneHours, double lcdDay, int lcdMonth, int lcdYear, double lctHour, double lctMin, double lctSec)
{
var haHours = PAMacros.RightAscensionToHourAngle(raHour, raMin, raSec, lctHour, lctMin, lctSec, daylightSaving, timezoneHours, lcdDay, lcdMonth, lcdYear, geogLongDeg);
var correctedHAHours = PAMacros.ParallaxHA(haHours, 0, 0, decDeg, decMin, decSec, coordinateType, geogLatDeg, heightM, equatorialHorParallaxDeg);
var correctedRAHours = PAMacros.HourAngleToRightAscension(correctedHAHours, 0, 0, lctHour, lctMin, lctSec, daylightSaving, timezoneHours, lcdDay, lcdMonth, lcdYear, geogLongDeg);
var correctedDecDeg1 = PAMacros.ParallaxDec(haHours, 0, 0, decDeg, decMin, decSec, coordinateType, geogLatDeg, heightM, equatorialHorParallaxDeg);
var correctedRAHour = PAMacros.DecimalHoursHour(correctedRAHours);
var correctedRAMin = PAMacros.DecimalHoursMinute(correctedRAHours);
var correctedRASec = PAMacros.DecimalHoursSecond(correctedRAHours);
var correctedDecDeg = PAMacros.DecimalDegreesDegrees(correctedDecDeg1);
var correctedDecMin = PAMacros.DecimalDegreesMinutes(correctedDecDeg1);
var correctedDecSec = PAMacros.DecimalDegreesSeconds(correctedDecDeg1);
return(correctedRAHour, correctedRAMin, correctedRASec, correctedDecDeg, correctedDecMin, correctedDecSec);
}
/// <summary>
/// Correct ecliptic coordinates for the effects of aberration.
/// </summary>
/// <returns>
/// apparent ecliptic longitude (degrees, minutes, seconds),
/// apparent ecliptic latitude (degrees, minutes, seconds)
/// </returns>
public (double apparentEclLongDeg, double apparentEclLongMin, double apparentEclLongSec, double apparentEclLatDeg, double apparentEclLatMin, double apparentEclLatSec) CorrectForAberration(double utHour, double utMinutes, double utSeconds, double gwDay, int gwMonth, int gwYear, double trueEclLongDeg, double trueEclLongMin, double trueEclLongSec, double trueEclLatDeg, double trueEclLatMin, double trueEclLatSec)
{
var trueLongDeg = PAMacros.DegreesMinutesSecondsToDecimalDegrees(trueEclLongDeg, trueEclLongMin, trueEclLongSec);
var trueLatDeg = PAMacros.DegreesMinutesSecondsToDecimalDegrees(trueEclLatDeg, trueEclLatMin, trueEclLatSec);
var sunTrueLongDeg = PAMacros.SunLong(utHour, utMinutes, utSeconds, 0, 0, gwDay, gwMonth, gwYear);
var dlongArcsec = -20.5 * ((sunTrueLongDeg - trueLongDeg).ToRadians()).Cosine() / ((trueLatDeg).ToRadians()).Cosine();
var dlatArcsec = -20.5 * ((sunTrueLongDeg - trueLongDeg).ToRadians()).Sine() * ((trueLatDeg).ToRadians()).Sine();
var apparentLongDeg = trueLongDeg + (dlongArcsec / 3600);
var apparentLatDeg = trueLatDeg + (dlatArcsec / 3600);
var apparentEclLongDeg = PAMacros.DecimalDegreesDegrees(apparentLongDeg);
var apparentEclLongMin = PAMacros.DecimalDegreesMinutes(apparentLongDeg);
var apparentEclLongSec = PAMacros.DecimalDegreesSeconds(apparentLongDeg);
var apparentEclLatDeg = PAMacros.DecimalDegreesDegrees(apparentLatDeg);
var apparentEclLatMin = PAMacros.DecimalDegreesMinutes(apparentLatDeg);
var apparentEclLatSec = PAMacros.DecimalDegreesSeconds(apparentLatDeg);
return(apparentEclLongDeg, apparentEclLongMin, apparentEclLongSec, apparentEclLatDeg, apparentEclLatMin, apparentEclLatSec);
}
/// <summary>
/// Calculate corrected RA/Dec, accounting for atmospheric refraction.
/// </summary>
/// <remarks>
/// NOTE: Valid values for coordinate_type are "TRUE" and "APPARENT".
/// </remarks>
/// <returns>
/// <para>corrected RA hours,minutes,seconds</para>
/// <para>corrected Declination degrees,minutes,seconds</para>
/// </returns>
public (double correctedRAHour, double correctedRAMin, double correctedRASec, double correctedDecDeg, double correctedDecMin, double correctedDecSec) AtmosphericRefraction(double trueRAHour, double trueRAMin, double trueRASec, double trueDecDeg, double trueDecMin, double trueDecSec, PACoordinateType coordinateType, double geogLongDeg, double geogLatDeg, int daylightSavingHours, int timezoneHours, double lcdDay, int lcdMonth, int lcdYear, double lctHour, double lctMin, double lctSec, double atmosphericPressureMbar, double atmosphericTemperatureCelsius)
{
var haHour = PAMacros.RightAscensionToHourAngle(trueRAHour, trueRAMin, trueRASec, lctHour, lctMin, lctSec, daylightSavingHours, timezoneHours, lcdDay, lcdMonth, lcdYear, geogLongDeg);
var azimuthDeg = PAMacros.EquatorialCoordinatesToAzimuth(haHour, 0, 0, trueDecDeg, trueDecMin, trueDecSec, geogLatDeg);
var altitudeDeg = PAMacros.EquatorialCoordinatesToAltitude(haHour, 0, 0, trueDecDeg, trueDecMin, trueDecSec, geogLatDeg);
var correctedAltitudeDeg = PAMacros.Refract(altitudeDeg, coordinateType, atmosphericPressureMbar, atmosphericTemperatureCelsius);
var correctedHAHour = PAMacros.HorizonCoordinatesToHourAngle(azimuthDeg, 0, 0, correctedAltitudeDeg, 0, 0, geogLatDeg);
var correctedRAHour1 = PAMacros.HourAngleToRightAscension(correctedHAHour, 0, 0, lctHour, lctMin, lctSec, daylightSavingHours, timezoneHours, lcdDay, lcdMonth, lcdYear, geogLongDeg);
var correctedDecDeg1 = PAMacros.HorizonCoordinatesToDeclination(azimuthDeg, 0, 0, correctedAltitudeDeg, 0, 0, geogLatDeg);
var correctedRAHour = PAMacros.DecimalHoursHour(correctedRAHour1);
var correctedRAMin = PAMacros.DecimalHoursMinute(correctedRAHour1);
var correctedRASec = PAMacros.DecimalHoursSecond(correctedRAHour1);
var correctedDecDeg = PAMacros.DecimalDegreesDegrees(correctedDecDeg1);
var correctedDecMin = PAMacros.DecimalDegreesMinutes(correctedDecDeg1);
var correctedDecSec = PAMacros.DecimalDegreesSeconds(correctedDecDeg1);
return(correctedRAHour, correctedRAMin, correctedRASec, correctedDecDeg, correctedDecMin, correctedDecSec);
}
/// <summary>
/// Calculate precise position of the sun for a local date and time.
/// </summary>
public (double sunRAHour, double sunRAMin, double sunRASec, double sunDecDeg, double sunDecMin, double sunDecSec) PrecisePositionOfSun(double lctHours, double lctMinutes, double lctSeconds, double localDay, int localMonth, int localYear, bool isDaylightSaving, int zoneCorrection)
{
var daylightSaving = (isDaylightSaving == true) ? 1 : 0;
var gDay = PAMacros.LocalCivilTimeGreenwichDay(lctHours, lctMinutes, lctSeconds, daylightSaving, zoneCorrection, localDay, localMonth, localYear);
var gMonth = PAMacros.LocalCivilTimeGreenwichMonth(lctHours, lctMinutes, lctSeconds, daylightSaving, zoneCorrection, localDay, localMonth, localYear);
var gYear = PAMacros.LocalCivilTimeGreenwichYear(lctHours, lctMinutes, lctSeconds, daylightSaving, zoneCorrection, localDay, localMonth, localYear);
var sunEclipticLongitudeDeg = PAMacros.SunLong(lctHours, lctMinutes, lctSeconds, daylightSaving, zoneCorrection, localDay, localMonth, localYear);
var raDeg = PAMacros.EcRA(sunEclipticLongitudeDeg, 0, 0, 0, 0, 0, gDay, gMonth, gYear);
var raHours = PAMacros.DecimalDegreesToDegreeHours(raDeg);
var decDeg = PAMacros.EcDec(sunEclipticLongitudeDeg, 0, 0, 0, 0, 0, gDay, gMonth, gYear);
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 distance to the Sun (in km), and angular size.
/// </summary>
/// <returns>
/// <para>sunDistKm -- Sun's distance, in kilometers</para>
/// <para>sunAngSizeDeg -- Sun's angular size (degrees part)</para>
/// <para>sunAngSizeMin -- Sun's angular size (minutes part)</para>
/// <para>sunAngSizeSec -- Sun's angular size (seconds part)</para>
/// </returns>
public (double sunDistKm, double sunAngSizeDeg, double sunAngSizeMin, double sunAngSizeSec) SunDistanceAndAngularSize(double lctHours, double lctMinutes, double lctSeconds, double localDay, int localMonth, int localYear, bool isDaylightSaving, int zoneCorrection)
{
var daylightSaving = (isDaylightSaving) ? 1 : 0;
var gDay = PAMacros.LocalCivilTimeGreenwichDay(lctHours, lctMinutes, lctSeconds, daylightSaving, zoneCorrection, localDay, localMonth, localYear);
var gMonth = PAMacros.LocalCivilTimeGreenwichMonth(lctHours, lctMinutes, lctSeconds, daylightSaving, zoneCorrection, localDay, localMonth, localYear);
var gYear = PAMacros.LocalCivilTimeGreenwichYear(lctHours, lctMinutes, lctSeconds, daylightSaving, zoneCorrection, localDay, localMonth, localYear);
var trueAnomalyDeg = PAMacros.SunTrueAnomaly(lctHours, lctMinutes, lctSeconds, daylightSaving, zoneCorrection, localDay, localMonth, localYear);
var trueAnomalyRad = trueAnomalyDeg.ToRadians();
var eccentricity = PAMacros.SunEcc(gDay, gMonth, gYear);
var f = (1 + eccentricity * trueAnomalyRad.Cosine()) / (1 - eccentricity * eccentricity);
var rKm = 149598500 / f;
var thetaDeg = f * 0.533128;
var sunDistKm = Math.Round(rKm, 0);
var sunAngSizeDeg = PAMacros.DecimalDegreesDegrees(thetaDeg);
var sunAngSizeMin = PAMacros.DecimalDegreesMinutes(thetaDeg);
var sunAngSizeSec = PAMacros.DecimalDegreesSeconds(thetaDeg);
return(sunDistKm, sunAngSizeDeg, sunAngSizeMin, sunAngSizeSec);
}
/// <summary>
/// Calculate precession (corrected coordinates between two epochs)
/// </summary>
/// <returns>Tuple (correctedRAHour, correctedRAMinutes, correctedRASeconds, correctedDecDeg, correctedDecMinutes, correctedDecSeconds)</returns>
public (double correctedRAHour, double correctedRAMinutes, double correctedRASeconds, double correctedDecDeg, double correctedDecMinutes, double correctedDecSeconds) CorrectForPrecession(double raHour, double raMinutes, double raSeconds, double decDeg, double decMinutes, double decSeconds, double epoch1Day, int epoch1Month, int epoch1Year, double epoch2Day, int epoch2Month, int epoch2Year)
{
var ra1Rad = (PAMacros.DegreeHoursToDecimalDegrees(PAMacros.HMStoDH(raHour, raMinutes, raSeconds))).ToRadians();
var dec1Rad = (PAMacros.DegreesMinutesSecondsToDecimalDegrees(decDeg, decMinutes, decSeconds)).ToRadians();
var tCenturies = (PAMacros.CivilDateToJulianDate(epoch1Day, epoch1Month, epoch1Year) - 2415020) / 36525;
var mSec = 3.07234 + (0.00186 * tCenturies);
var nArcsec = 20.0468 - (0.0085 * tCenturies);
var nYears = (PAMacros.CivilDateToJulianDate(epoch2Day, epoch2Month, epoch2Year) - PAMacros.CivilDateToJulianDate(epoch1Day, epoch1Month, epoch1Year)) / 365.25;
var s1Hours = ((mSec + (nArcsec * (ra1Rad).Sine() * (dec1Rad).Tangent() / 15)) * nYears) / 3600;
var ra2Hours = PAMacros.HMStoDH(raHour, raMinutes, raSeconds) + s1Hours;
var s2Deg = (nArcsec * (ra1Rad).Cosine() * nYears) / 3600;
var dec2Deg = PAMacros.DegreesMinutesSecondsToDecimalDegrees(decDeg, decMinutes, decSeconds) + s2Deg;
var correctedRAHour = PAMacros.DecimalHoursHour(ra2Hours);
var correctedRAMinutes = PAMacros.DecimalHoursMinute(ra2Hours);
var correctedRASeconds = PAMacros.DecimalHoursSecond(ra2Hours);
var correctedDecDeg = PAMacros.DecimalDegreesDegrees(dec2Deg);
var correctedDecMinutes = PAMacros.DecimalDegreesMinutes(dec2Deg);
var correctedDecSeconds = PAMacros.DecimalDegreesSeconds(dec2Deg);
return(correctedRAHour, correctedRAMinutes, correctedRASeconds, correctedDecDeg, correctedDecMinutes, correctedDecSeconds);
}
/// <summary>
/// Convert Equatorial Coordinates to Galactic Coordinates
/// </summary>
/// <returns>Tuple (galLongDeg, galLongMin, galLongSec, galLatDeg, galLatMin, galLatSec)</returns>
public (double galLongDeg, double galLongMin, double galLongSec, double galLatDeg, double galLatMin, double galLatSec) EquatorialCoordinateToGalacticCoordinate(double raHours, double raMinutes, double raSeconds, double decDegrees, double decMinutes, double decSeconds)
{
var raDeg = PAMacros.DegreeHoursToDecimalDegrees(PAMacros.HMStoDH(raHours, raMinutes, raSeconds));
var decDeg = PAMacros.DegreesMinutesSecondsToDecimalDegrees(decDegrees, decMinutes, decSeconds);
var raRad = raDeg.ToRadians();
var decRad = decDeg.ToRadians();
var sinB = decRad.Cosine() * (27.4).ToRadians().Cosine() * (raRad - (192.25).ToRadians()).Cosine() + decRad.Sine() * (27.4).ToRadians().Sine();
var bRadians = sinB.ASine();
var bDeg = PAMacros.Degrees(bRadians);
var y = decRad.Sine() - sinB * (27.4).ToRadians().Sine();
var x = decRad.Cosine() * (raRad - (192.25).ToRadians()).Sine() * (27.4).ToRadians().Cosine();
var longDeg1 = PAMacros.Degrees(y.AngleTangent2(x)) + 33;
var longDeg2 = longDeg1 - 360 * (longDeg1 / 360).Floor();
var galLongDeg = PAMacros.DecimalDegreesDegrees(longDeg2);
var galLongMin = PAMacros.DecimalDegreesMinutes(longDeg2);
var galLongSec = PAMacros.DecimalDegreesSeconds(longDeg2);
var galLatDeg = PAMacros.DecimalDegreesDegrees(bDeg);
var galLatMin = PAMacros.DecimalDegreesMinutes(bDeg);
var galLatSec = PAMacros.DecimalDegreesSeconds(bDeg);
return(galLongDeg, galLongMin, galLongSec, galLatDeg, galLatMin, galLatSec);
}