Selbst rechnen
----------------
In der Astrophysik-Ecke habe ich ein Programm zur
Berechnung des Zentrallinie vorgestellt, erst mal ohne Erklärungen. Es dient eigentlich zur Vorbereitung des Sonnenfinsternis 2026. Wie angedeutet, kommen die erläuterungen dazu in den Journalen des Jahres 2026.
Mit diesem Programm lassen sich historische Finsternisse nachrechnen, wenn man die Bessel-Elemente hat. Wichtig: Gleich am Anfang des Hauptprogramms muss die Zeitabweichung ΔT eingesetzt werden
dt=16953
Diese gilt für alle Finsternisse in dieser Zeit. Dann braucht man noch die Bessel-Elemente: Ich habe sie für alle Frühjahrsfinsternisse der fraglichen Zeit herausgesucht:
Code:
yr=-480
mnt=4
day=19
T0 = 9
x = [ -0.1491210, 0.5499079, 0.0000129, -0.0000089]
y = [ 0.2520590, 0.1620816, -0.0000729, -0.0000025]
d = [ 9.0285900, 0.0146180, -0.0000020, 0.0]
l1 = [ 0.5353800, 0.0000837, -0.0000124, 0.0]
l2 = [ -0.0107040, 0.0000833, -0.0000123, 0.0]
u = [ 315.5606100, 15.0049000, 0.0000000, 0.0]
tanf1 = 0.0046113
tanf2 = 0.0045884
yr=-479
mnt=4
day=9
T0 = 2
x = [ 0.1971560, 0.5565842, 0.0000009, -0.0000095]
y = [ -0.4111460, 0.1764906, -0.0000062, -0.0000029]
d = [ 5.1239600, 0.0153990, -0.0000010, 0.0]
l1 = [ 0.5314260, -0.0000184, -0.0000129, 0.0]
l2 = [ -0.0146390, -0.0000183, -0.0000128, 0.0]
u = [ 209.5904800, 15.0053800, 0.0000000, 0.0]
tanf1 = 0.0046211
tanf2 = 0.0045981
yr=-478
mnt=2
day=28
T0 = 8
x = [ -0.3093620, 0.5276379, 0.0000152, -0.0000075]
y = [ 1.4298691, 0.1516152, -0.0000019, -0.0000021]
d = [ -10.2092100, 0.0147030, 0.0000020, 0.0]
l1 = [ 0.5489280, -0.0001269, -0.0000115, 0.0]
l2 = [ 0.0027770, -0.0001263, -0.0000114, 0.0]
u = [ 296.0756800, 15.0036100, 0.0000000, 0.0]
tanf1 = 0.0046724
tanf2 = 0.0046491
yr=-478
mnt=3
day=29
T0 = 18
x = [ 0.3977080, 0.5376028, -0.0000053, -0.0000084]
y = [ -1.1365520, 0.1775236, 0.0000602, -0.0000027]
d = [ 1.0502000, 0.0158010, -0.0000010, 0.0]
l1 = [ 0.5394540, -0.0001037, -0.0000121, 0.0]
l2 = [ -0.0066500, -0.0001032, -0.0000120, 0.0]
u = [ 88.5397800, 15.0054300, 0.0000000, 0.0]
tanf1 = 0.0046330
tanf2 = 0.0046099
yr=-477
mnt=2
day=17
T0 = 15
x = [ 0.0025210, 0.5021353, -0.0000057, -0.0000060]
y = [ 0.8112590, 0.1221124, 0.0000536, -0.0000014]
d = [ -14.0009700, 0.0134030, 0.0000030, 0.0]
l1 = [ 0.5649260, -0.0001018, -0.0000103, 0.0]
l2 = [ 0.0186950, -0.0001013, -0.0000102, 0.0]
u = [ 40.6422900, 15.0021900, 0.0000000, 0.0]
tanf1 = 0.0046886
tanf2 = 0.0046652
yr=-476
mnt=2
day=6
T0 = 15
x = [ -0.0405820, 0.4952079, -0.0000183, -0.0000055]
y = [ 0.0758150, 0.0926106, 0.0001019, -0.0000010]
d = [ -17.4424710, 0.0114930, 0.0000040, 0.0]
l1 = [ 0.5723300, 0.0000021, -0.0000098, 0.0]
l2 = [ 0.0260610, 0.0000021, -0.0000098, 0.0]
u = [ 40.6103900, 15.0005400, 0.0000000, 0.0]
tanf1 = 0.0047046
tanf2 = 0.0046812
yr=-475
mnt=1
day=25
T0 = 16
x = [ 0.2204360, 0.5155033, -0.0000381, -0.0000063]
y = [ -0.5989650, 0.0626911, 0.0001538, -0.0000007]
d = [ -20.2609410, 0.0090140, 0.0000050, 0.0]
l1 = [ 0.5655600, 0.0000876, -0.0000104, 0.0]
l2 = [ 0.0193250, 0.0000871, -0.0000104, 0.0]
u = [ 56.0337300, 14.9989100, 0.0000000, 0.0]
tanf1 = 0.0047192
tanf2 = 0.0046957
Für die Finsternis des Jahres -477 taucht in der Ergebnistabelle die Nähe von Sardes (φ=38,48°N, λ=28,03°O) mit auf:
Code:
Zeit Lambda phi Dauer Hohe Pfad Monddurchm. Typ
°,<0:O °,+:N s ° (km)
---------------------------------------------------------------
-6 -28.1569 39.3621 355 36.5 314.2 0.945 r
Jetzt geht das Zeit-Hickhack wieder los:
Die Zeit ist gerechnet in Minuten von der Referenzstunde , es ergibt sich 15h-6m=14h54m TD. Hier muss man ΔT berücksichtigen, das sind 4h43m. UT=TD-ΔT, also fand die Finsternis 10h11m UT statt. Sardes liegt östlich von Greenwich, und zwar geht die Sonne dort 1h52m früher auf, der Tag ist weiter fortgeschritten. Die Finsternis war also dort 12h3m mittlere Sardes-Zeit. Die Zeitgleichung ist um diese Jahreszeit fast maximal negativ, ca. 14 Minuten. Damit erhalten wir 11h49m wahrer Ortszeit.
Zech gibt die zentrale Zeit nicht an, nur Anfang und Ende der ringförmigen Phase. Der Mittelwert ist 23h52m w. Sard. Z. (wie es im Text heißt). Hiervon müssen wir 12h subtrahieren, da Zech die Tage ab mittag rechnet. Es folgt 11h52m bei ihm, nur 3 Minuten Differenz zu unserer Rechnung.
Voilá!