Das ist mir zu viel Arbeit.
Hm, also bei mir waren das vielleicht 10 Minuten. Den Text hier zu tippen dauert länger. Ich hab mir von CalSky die Daten geschnappt und den Winkel chi - q, also den Zenitwinkel des beleuchteten Mondteils berechnet.
Mond, 14.1.2019, 23:00 MEZ = 22:00 UT
Länge: lambda = 10°E
Breite: phi = 50°N
(= in der Nähe von Würzburg)
Mond topozentrisch:
RA = 2h 04m 0.6s = 31.0025°
DEC = 6°30'19" = 6.505°
Sternzeit in GW:
theta0 = 7h 32m 41s = 7.54472h = 113.171°
Stundenwinkel des Mondes:
H = theta0 + lambda - RA = 113.171° + 10° - 31° = 92.171°
Damit folgt für den parallaktischen Winkel mit der Formel für
tan(q) = 0.99928/1.1884 = 0.841
=> q = arctan(0.841) = 40.06°
Für Winkel phi benötigt man die geozentrischen Sonnen- und Mondkoordinaten.
CalSky sagt
Sonne geozentrisch:
RA0 = 19h 45m 24.9s = 296.354°
DEC0 = -21° 13m 39S = -21.228°
Mond geozentrisch:
RA = 2h 06m 5.2s = 31.522°
DEC = 7°09'16" = 7.154°
Damit folgt für tan(chi) = -0.928358/-0.34881 (3. Quadrant) = 2.6615
chi = arctan(2.6615) + 180° = 249.41°
Und damit
chi - q = 249.41° - 40.06° = 209.35°
Diese 209° werden vom Zenitpunkt gerechnet, und wenn ich mir die Grafik von GUIDE ansehe, dann kommt das auch hin, also die Mondsichel liegt ca. 30° zum Horizont, von links unten nach rechts oben. Das sind von der positiven mathematischen x-Achse aus dann 209.35° + 90° = 299.35° ≈ 300°
Mit der obigen Angabe ergibt Peters Formel Folgendes:
(CalSky berechnet den Azimut von Nord, also für Südrichtung 180° abziehen! Ist aber egal, da ohnehin die Differenz der Azimutwerte benötigt wird)
Sonne:
Az = x1 = 320°27'15" = 320.454° - 180° = 140.454°
h = y1 = -56°17'09" = -56.286°
Mond:
Az = x2 = 252°56'17" = 252.938° - 180° = 72.938°
h = y2 = 22°05'15" = 22.0875°
tan(alpha) = tan(y2 - y1) / sin(x2 - x1)
tan(alpha) = tan(22.0875° - (-56.286°)) / sin(72.938° - 140.454°)
tan(alpha) = 4.860 / -0.924 (4. Quadrant)
alpha = -79.236° = 280.764°
Die Mondsichel müsste nach dieser Formel also ca. 11° gegen den Horizont geneigt sein.
cs,
harald
--