Suche: Skript für Position des Mondes

slart

Aktives Mitglied
Hallo Gemeinde,

Für ein bisschen kreatives Ausleben suche ich ein Skript, das die Position des Mondes auf der Erde berechnet. Da gibt es dieses schon etwas ältere Skript, das aber scheinbar einen Fehler hat?
Der Mond bewegt sich diesem Skript nach nur 1° in 2 Stunden? Gibt es ähnliche Skripte mit plausibleren Ergebnissen?
 
Hallo slart,

soweit ich sehe ist das eine Umsetzung von Jean Meeus' Algorithmus aus seinem Buch "Astronomical Algorithms".

Ich habe jetzt mal mit dem Testdatum 12.04.1992, 00:00 TD getestet (19920412.0000), das Meeus auch im Buch verwendet, es stimmen Rektaszension und Deklination überein, allerdings nicht die geozentrische Entfernung. Dieses Script liefert Δ = 368417 km, das korrekte Ergebnis laut Script wäre aber Δ = 368409.7 km, und das noch genauere (mit der ungekürzten ELP-2000/82) wäre Δ = 368405.6 km

Zitat von slart:
Der Mond bewegt sich diesem Skript nach nur 1° in 2 Stunden?
Klicke in dieses Feld, um es in vollständiger Größe anzuzeigen.

Woraus schließt du das?

Zitat von slart:
Gibt es ähnliche Skripte mit plausibleren Ergebnissen?
Klicke in dieses Feld, um es in vollständiger Größe anzuzeigen.

Sicher, siehe EphemPedia: Sonne, Mond, geozentrisch

cs,
harald

--
 
Stimmt, ich vergaß die alte Bauernregel "....ungefähr seinen scheinbaren Durchmesser pro Stunde".
Und ~ bedeutet "proportional", ≈ hingegen "ungefähr" ;) :coffee:

cs,
harald

--
 
> suche ich ein Skript, das die Position des Mondes auf der Erde berechnet.

...wenn du dieses erstaunliche Script gefunden hast, laß es uns wissen! :unsure: :sneaky:

*Durch Moderation aufgrund fehlender Copyrightangabe entfernt*

(moon landing on earth)
 
Zuletzt von einem Moderator bearbeitet:
Hallo h_c_greier und P_E_T_E_R,

vielen Dank für die Antworten und Hinweise.

Das gesuchte Skript soll das Ergebnis der Mondposition in equatorial liefern. Das von mir eingangs verlinkte Skript macht dies, auch wenn Longitude nur in Hours (was in Degree umgerechnet werden kann). Aber die Ergebnisse kommen mir seltsam vor.

Zur leichteren Erklärung im Folgenden "fliegt" der Mond (von der Erde aus gesehen) in etwa 24h50min einmal um die Erde. Tatsächlich (vom Mond aus gesehen) in 29.5 Tagen.

Eine "Umdrehung" des Mondes um die Erde von der Erde aus gesehen, also dass der Mond wieder an der selben Longitude ist, beträgt im Ø 24h50min.
Ergebnisse des Skripts:
20250928.1230
- Lon: (ec): 259.253 Deg
- Lon (eq): 17.1888 Hrs
20250928.1430
- Lon: (ec): 260.25 Deg
- Lon (eq): 17.2638 Hrs
20250928.1930
- Lon: (ec): 262.748 Deg
- Lon (eq): 17.4519 Hrs
20250929.0230
- Lon: (ec): 266.257 Deg
- Lon (eq): 17.7169 Hrs
20250929.1230
- Lon: (ec): 271.296 Deg
- Lon (eq): 18.0979 Hrs
Der Mond hat sich nach diesen Ergebnissen von der Erde aus gesehen in 24 Stunden nur 22 Deg bewegt, müsste aber einmal fast drumherum "geflogen" sein. Es sind also scheinbar andere Ergebnisse als die von mir gesuchten?

EDIT: jetzt habe ich das Skript verstanden. Es berechnet die Position des Mondes um die Erde herum vom Mond aus gesehen, also ohne Erdrotation mit einem Koordinatensystem einer nicht rotierenden Erde.

Gesucht wird ein Skript, dass die Position des Mondes "von der Erde aus" berechnet, mit einem Ergebnis in equatorial sowie Erdkoordinatensystem.
 
Zuletzt bearbeitet:
> ....und wenn diese erstaunliche Wortspende irgendwann mal Sinn ergibt

Du verstehst auch gar keinen Spaß, was?! :cautious:
 
Moment mal. Kann es sein, dass die Longitude-Werte des eingangs von mir verlinkten Skripts ein Äquivalent zur "equation of time" für die Sonne sind, das an die aktuelle Uhrzeit angefügt werden muss? Während die "equation of time" für die Sonne meist nur für je einen ganzen Tag berechnet wird, werden die Werte für den Mond pro Minute berechnet, weil der mehr schwankt. Aber an welche Zeit müssten die Werte addiert/subtrahiert werden?
 
Versuch einer Lösung: Erweiterung des Skripts Moon position Javascript: finds the apparent position of the Moon to within 10 arcsec allowing for nutation and aberration. Based on the truncated version of the ELP-2000 series by Jean Meeus. Includes Javascript code from Arkkana Peck's 'Hitchhiker's Guide to the Moon' page.


Die Funktion Aries() Neu hinzufügen:

// https://repositorio.unican.es/xmlui/bitstream/handle/10902/13373/Pascual Cabañero, María.pdf
function Aries(day) {
GHAAtrue = range(range(280.46061837 + 360.98564736629 * (day - 2451545) + 0.000387933 * (t*t) - (t*t*t) / 38710000) + delta_phi * Math.cos(eps));
return GHAAtrue;
}



Die Funktion Moon() erweitern:

function Moon(day) {

// ...

GHAAtrue = Aries(day);

// [gelöst] Kennt jemand Mond/Sonnen-Komponenten? - Deutsches Lazarusforum
var lambda = Lprime + sigma_l / 1000000;
var beta = sigma_b / 1000000;

// Moon position Javascript: finds the apparent position of the Moon to within 10 arcsec allowing for nutation and aberration. Based on the truncated version of the ELP-2000 series by Jean Meeus. Includes Javascript code from Arkkana Peck's 'Hitchhiker's Guide to the Moon' page.
Lprime_range = range(Lprime);
var lambda_range = range(Lprime_range + sigma_l / 1000000);

var lat = lambda_range + delta_phi;

// ecliptic to equatorial
// Calculating azimuth from equatorial coordinates
var RA = Math.atan2((Math.cos(eps) * Math.sin(lat) - Math.sin(eps) * Math.tan(beta)), Math.cos(lat));
var DEC = Math.asin(Math.cos(eps) * Math.sin(beta) + Math.sin(eps) * Math.cos(beta) * Math.sin(lat));

// https://repositorio.unican.es/xmlui/bitstream/handle/10902/13373/Pascual Cabañero, María.pdf
// Star()
RAmoon = range(RA / DEG2RAD);
DECmoon = (DEC / RAD2DEG); // Latitude
GHAmoon = range(GHAAtrue - RAmoon); // Longitude
}


Die Funktion nutation() erweitern:

var delta_phi = 0;
var eps = 0;
function nutation(day) {

// ...

// https://repositorio.unican.es/xmlui/bitstream/handle/10902/13373/Pascual Cabañero, María.pdf
delta_phi = dphi / 36000000;
delta_eps = deps / 36000000;
eps = ((84381.448 - 46.815 * t - 0.00059 * (t*t) + 0.001813 * (t*t*t)) / 3600) + delta_eps;
}


Keine Validierung. Bitte selbst testen.
 
Um antworten zu können musst du eingeloggt sein.
Zurück
Oben