Sonnenbahnberechnung mit BASCOM

[mat-sche]

Hi Dolivo,

Dein Code sieht ganz chick aus! Nun würde ich diesen gern bei mir in meine Haussteuerung einsetzen. Leider nur kann ich nicht einen Bezugspunkt im Code finden, bei dem ich meine Standortkoordinaten einfügen kann. Könntest Du mir bitte mal kurz unter die Arme greifen und helfen .
Danke und Grüße MAT

[dolivo]

Hallo,
ich habe die im lexikon.astronomie aufgeführten Formeln einfach nur "aufgetröselt", damit sie in Bascom verarbeitet werden können. Wenn Du Dir das Zahlenbeispiel von Berlin ansiehst, wird auch manches klarer. Ich hänge aber noch einmal mein Programm an, in dem ich ein paar Erläuterungen mehr eingefügt habe.

[mat-sche]

Hi,

dank für Deine Antwort!
habe mal versucht den Code im Simulator laufen zu lassen, natürlich mit meinen Koordinaten.
Was ich bei Deinem Code nicht so recht verstehe ist:
'Sonnenaufgang
Sa = 12 - Zeitdiff
Sa = Sa - Woz_moz
Sa = Sa - 0.9067 '13°37'/15°=0,9067
...wie kommst Du auf 0,9067?

Und dann ist da noch... welche Zeit wird genau ausgerechnet? Wenn ich die Ergebnisse mit http://www.suncalc.net/#/51.7236,12....12.10.01/18:37 vergleiche, gibt es mehrere Arten der Sonnen auf/untergänge...
hmmmmm.... werde wohl noch weiter testen.

[dolivo]

Hallo, mat-sche,

ja, da habe ich falsch gerechnet. Richtig ist 0,9078.
Ich habe meine Werte bisher nur mit den Zahlen verglichen, die im Kalender stehen und die Zeitdifferenz zwischen meinem Ort und dem angegebenen (hier Kassel, ermittelt mit der Längengraddifferenz) abgezogen (da östlich). Das stimmt recht gut überein innerhalb des genannten Fehlerbereichs.

Nun noch zum Beitrag von U-1110 v. 28.8.: Ich berücksichtige die Sommer- bzw. Normalzeit beim Sonnenauf- bzw. -untergang. An der Stelle interessiert mich ja nur das Zwischenergebnis.

dolivo

[mat-sche]

Hi Dolivo,

nun endlich habe ich Deinen Code in meine Haussteuerung integriert und nun läuft die Berechnung dank dcf77 um die 1min genau.
Danke für die Codezeilen!
Grüße MAT

[Accenter]

Bin gerade auf diesen Code hier gestoßen und frage mich wie ich das ganze für Wien umschreiben kann?

Meine Werte sind:
Breitengrad : 48.2667
Längengrad : 16.4


Würde mich über Hilfe freuen, ich verstehe nur Teile von diesem Code.

[dolivo]

Ich habe die Erklärung etwas kurz gelassen, weil die Berechnung im Astrolexikon gut beschrieben ist. Ich hatte das ja nur auf Bascom umgesetzt und die Werte, die für den Ort konstant sind, gleich eingesetzt um die Rechenzeit zu kürzen.
Für Wien muss man bei der Berechnung der Zeitdifferenz so rechnen:
B=pi*48,267/180 = 0,843, sin B= 0,747 (natürlich mit Bogenmaß rechnen);
Sonnenaufgang: 16,4/15 = 1,093
Viel Erfolg
dolivo

[Accenter]

Danke! Teste ich gleich mal

[Accenter]

Ich habe jetzt mal den Code nach deinen Werten die du mir errechnet hast umgebaut. Die Sonnenuntergangszeit stimmt auf die Minute genau, jedoch liegt die Sonnenaufgangszeit um 10 Minuten daneben.
Was habe ich falsch gemacht?

Hier mal mein Code:

 Zone = 1
 Term1 = 0.0337 * Tagimjahr
 Term1 = Term1 + 0.465
 Term1 = Sin(term1)
 Term1 = -0.171 * Term1
 Term2 = 0.01787 * Tagimjahr
 Term2 = Term2 - 0.168
 Term2 = Sin(term2)
 Term2 = 0.1299 * Term2
 Woz_moz = Term1 - Term2                                    'Zeitgleichung /in Stunden
 Diff2 = Woz_moz * 60                                       'auf Minuten gerechnet

 'Berechnung Deklination
 Term1 = Tagimjahr - 80.086
 Term1 = Term1 * 0.016906
 Term1 = Sin(term1)
 Deklination = Term1 * 0.4095

 'Berechnung Zeitdifferenz
 Term1 = Sin(deklination)
 Term1 = 0.747 * Term1                                      'B=pi*48,267/180 = 0,843, sin B= 0,747
 Term1 = -0.0145 - Term1                                    'sin(h)=-0,0145
 Term2 = Cos(deklination)
 Term2 = 0.6024 * Term2                                     'cos(B)=0,6294
 Zeitdiff = Term1 / Term2
 Zeitdiff = Acos(zeitdiff)
 Zeitdiff = 3.8197 * Zeitdiff                               '12/pi=3,8197

 'Sonnenaufgang
 Sa = 12 - Zeitdiff
 Sa = Sa - Woz_moz
 Sa = Sa - 1.093                                            '16.4 / 15 = 1.093
 Sa = Sa + Zone                                             'MEZ, bei MESZ +2
 Samin = Frac(sa)
 Sahour = Sa - Samin
 Samin1 = Samin * 60
 Sastrmin = Str(samin1)
 Sastrmin = Format(sastrmin , "00")

 'Sonnenuntergang
 Su = 12 + Zeitdiff
 Su = Su - Woz_moz
 Su = Su - 0.9067
 Su = Su + Zone
 Sumin = Frac(su)
 Suhour = Su - Sumin
 Sumin1 = Sumin * 60
 Sustrmin = Str(sumin1)
 Sustrmin = Format(sustrmin , "00")

[dolivo]

Natürlich muss auch bei der Sonnenuntergangsberechnung Dein Wert von Wien eingesetzt werden. Also auch 1,093 statt o,9067!

'Sonnenuntergang Su = 12 + Zeitdiff
Su = Su - Woz_moz
Su = Su - 0.9067
...