private void calculateRiseSetMethodM(double D, GPLocationProvider ed)
{
GPLocation obs = ed.getLocation(D);
double a1, a2, a3;
double d1, d2, d3;
double siderealTime = GPAstroEngine.GetSiderealTime(D);
double h0 = -0.833333;
calculateCoordinatesMethodM(D - 1);
a1 = rightAscession;
d1 = declination;
calculateCoordinatesMethodM(D);
a2 = rightAscession;
d2 = declination;
calculateCoordinatesMethodM(D + 1);
a3 = rightAscession;
d3 = declination;
double longitude = -ed.GetLongitudeEastPositive();
double latitude = ed.GetLatitudeNorthPositive();
double cosH0 = (GPMath.sinDeg(h0) - GPMath.sinDeg(latitude) * GPMath.sinDeg(d2))
/ (GPMath.cosDeg(latitude) * GPMath.cosDeg(d2));
double H0 = GPMath.arccosDeg(cosH0);
H0 = GPMath.putIn180(H0);
double m0 = (a2 + longitude - siderealTime) / 360;
double m1 = m0 - H0 / 360;
double m2 = m0 + H0 / 360;
double deltaM = 0;
deltaM = getCorrection(D, a1, a2, a3, d1, d2, d3, siderealTime, h0, longitude, latitude, m0, true);
m0 += deltaM;
deltaM = getCorrection(D, a1, a2, a3, d1, d2, d3, siderealTime, h0, longitude, latitude, m1, false);
m1 += deltaM;
deltaM = getCorrection(D, a1, a2, a3, d1, d2, d3, siderealTime, h0, longitude, latitude, m2, false);
m2 += deltaM;
julianDayRise = julianDay + m1;
julianDayNoon = julianDay + m0;
julianDaySet = julianDay + m2;
sunrise_deg = GPMath.putIn360(m1 * 360);
noon_deg = GPMath.putIn360(m0 * 360);
sunset_deg = GPMath.putIn360(m2 * 360);
}
