#sonne.bas '-- berechnet Sonnenauf-/-untergang '-- minutengenau in sa$ und su$ '-- Format hhmm00 '-- Geogr. Koordinaten br=53.0961 la=8.4658 '- Bogenmaß brad=br/180*pi '- Höhe Sonne -50 Bogenmin. h0=-50/60/180*pi '- Deklination dek, Tag des Jahres d0 year$=left$(fn date$(0),4) d0=fn days(fn date$(0)) - fn days(year$+"0101")+1 dek=0.40954 * sin(0.0172 * (d0 - 79.35)) zh1=sin(h0) - sin(brad) * sin(dek) zh2=cos(brad) * cos(dek) zd=12*acos(zh1/zh2) / pi zgl=-0.1752 * sin(0.03343 * d0 + 0.5474) - 0.134 * sin(0.018234 * d0 - 0.1939) '-Sonnenuntergang tsu=12 + zd - zgl su=(tsu + (15 - la) /15) hh$=str$(int(su,0)) if len(hh$)<2 then hh$="0"+hh$ endif su = (su - val(hh$))*60 mm$=str$(int(su,0)) if len(mm$)<2 then mm$="0"+mm$ endif su$=hh$+mm$+"00" '- Sonnenaufgang tsa=12 - zd - zgl sa=(tsa + (15 - la) /15) hh$=str$(int(sa,0)) if len(hh$)<2 then hh$="0"+hh$ endif sa = (sa - val(hh$))*60 mm$=str$(int(sa,0)) if len(mm$)<2 then mm$="0"+mm$ endif sa$=hh$+mm$+"00" '---Ende