#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