czas temperaturę imieniny i inne wydarzenia i pomyślałem że fajnie by było aby pokazywał
też godzinę wschodu i zachodu słońca oraz długość dnia.
Po kilku dniach wertowania internetu i studiowaniu rożnych wzorów powstał taki oto kawałek kodu
[syntax=basic4gl]$regfile = "m328pdef.dat"
$crystal = 16000000
$hwstack = 100
$swstack = 100
$framesize = 40
$baud = 57600
Config Serialin = Buffered , Size = 250
Config Serialout = Buffered , Size = 250
Config Pind.2 = Output
Config Pind.3 = Output
Reset Portd.2
Reset Portd.3
Enable Interrupts
'w programie wykorzystałem część wzorów podanych na stronie http://cybermoon.pl/wiedza/algorithms/wschody_slonca.html
'inne zmodyfikowałem lub przekształciłem
' wyniki są z dokładnością do jednej dwóch minut
'czas wykonywania samych obliczeń przy kwarcu 16 MHz to około 3,5milisekundy
'wynik jest podany w czasie UTC dla Polski należy doliczyć dla czasu zimowego 1 godzinę a dla letniego 2 godziny
'*******DANE WEJŚCIOWE ROK MIESIĄC DZIEŃ SZEROKOŚĆ I DŁUGOŚĆ GEOGRAFICZNA
'Zmienna Req w zależności jaki wschód i zachód liczymy (połozenie słońca pod horyzontem)
'=-0.833 oficjalny
'=-6 cywilny
'=-12 żeglarski
'=-18 astronomiczny
Dim Req As Single : Req = -0.833
Dim Rok As Word : Rok = 2019
Dim Miesiac As Byte : Miesiac = 1
Dim Dzien As Byte : Dzien = 2
Dim Szerokosc As Single : Szerokosc = 50.069495
Dim Dlugosc As Single : Dlugosc = 20.480322
'****ZMIENNE POMOCNICZE
Dim J As Single : Dim J1 As Single : Dim J2 As Single
Dim Cent As Single
Dim L As Single : Dim L1 As Single
Dim G As Single : Dim G1 As Single
Dim O As Single
Dim F As Single : Dim F1 As Single : Dim F2 As Single
Dim E As Single : Dim E1 As Single : Dim E2 As Single
Dim A As Single : Dim A1 As Single : Dim A2 As Single
Dim C As Single : Dim C1 As Single : Dim C2 As Single
Dim C3 As Single : Dim C4 As Single : Dim C5 As Single
Dim C6 As Single : Dim C7 As Single : Dim C8 As Single
Dim Wsch As Single : Dim Wsch1 As Single : Dim Wsch2 As Single
Dim Godzina As Single : Dim Minuta As Single
Dim Zach As Single : Dim Zach1 As Single : Dim Zach2 As Single
Dim Godzinazach As Single : Dim Minutazach As Single
Dim Gzach As Byte 'GODZINA ZACHODU
Dim Mzach As Byte 'MINUTA ZACHODU
Dim Gwsch As Byte 'GODZINA WSCHODU
Dim Mwsch As Byte 'MINUTA WSCHODU
Do
J1 = Miesiac + 9 : J1 = J1 / 12 : J1 = Int(j1) : J1 = J1 + Rok
J1 = J1 / 4 : J1 = J1 * 7 : J1 = Int(j1)
J2 = 275 * Miesiac : J2 = J2 / 9 : J2 = Int(j2)
J = 367 * Rok : J = J - J1 : J = J + J2
J = J + Dzien : J = J - 730531.5
Cent = J / 36525
L = 628.331969753 * Cent : L = L + 4.894950420 : L1 = L / 6.283185307
L1 = Int(l1) : L1 = L1 * 6.28318530 : L = L - L1
G = 628.30 * Cent : G = G + 6.2400408 : G1 = G / 6.283185307
G1 = Int(g1) : G1 = G1 * 6.283185307 : G = G - G1
O = 0.0002269 * Cent : O = 0.409093 - O
F1 = 2 * G : F1 = Sin(f1) : F1 = F1 * 0.00034907
F2 = Sin(g) : F2 = F2 * 0.033423 : F = F1 + F2
E2 = L + F : E2 = E2 * 2 : E2 = Sin(e2) : E2 = E2 * 0.0430398
E1 = L + F : E1 = E1 * 4 : E1 = Sin(e1) : E1 = E1 * 0.00092502
E = E2 - E1 : E = E - F
A1 = L + F : A1 = Sin(a1) : A = Sin(o)
A = A * A1 : A = Asin(a)
C1 = 0.017453293 * Req : C1 = Sin(c1)
C2 = 0.017453293 * Szerokosc : C2 = Sin(c2) : C3 = Sin(a)
C4 = 0.0174553293 * Szerokosc : C4 = Cos(c4) : C5 = Cos(a)
C6 = C2 * C3 : C7 = C1 - C6 : C8 = C4 * C5 : C = C7 / C8
Wsch1 = Acos(c) : Wsch2 = Dlugosc * 0.0174453293 : Wsch2 = Wsch2 + E
Wsch2 = Wsch2 + Wsch1 : Wsch = 3.141592653 - Wsch2
Wsch = Wsch * 57.2957795 : Wsch = Wsch / 15
Godzina = Int(wsch) : Gwsch = Godzina : Minuta = Wsch - Gwsch
Minuta = Minuta * 100 : Minuta = Int(minuta) : Minuta = Minuta * 6
Minuta = Minuta / 10 : Mwsch = Minuta
Zach1 = Acos(c) : Zach1 = -1 * Zach1 : Zach2 = Dlugosc * 0.0174453293
Zach2 = Zach2 + E : Zach2 = Zach2 + Zach1 : Zach = 3.141592653 - Zach2
Zach = Zach * 57.2957795 : Zach = Zach / 15
Godzinazach = Int(zach) : Gzach = Godzinazach : Minutazach = Zach - Gzach
Minutazach = Minutazach * 100 : Minutazach = Int(minutazach)
Minutazach = Minutazach * 6 : Minutazach = Minutazach / 10 : Mzach = Minutazach
Print "Rok--------- " ; Rok
Print "Miesiac----- " ; Miesiac
Print "Dzien------- " ; Dzien
Print "dlugosc----- " ; Dlugosc
Print "szerokosc--- " ; Szerokosc
Print "Wsch-------- " ; Gwsch ; ":" ; Mwsch
Print "zach-------- " ; Gzach ; ":" ; Mzach
Print " "
Wait 1
Loop[/syntax]
Może komuś się przyda w jakiejś stacji meteo, może ktoś to bardziej zoptymalizuje.
Oczywiście obliczenia nie przewidują nocy i dnia polarnego (należałoby sprawdzać najpierw szerokość geograficzną i wtedy po wystąpieniu błędu sprawdzić czy to jest noc czy dzień polarny)
Dla szerokości geograficznej południowej i długości zachodniej dane wprowadzamy w postaci liczb ujemnych.
Ogólnie wyniki generowane dla miejscowości w Polsce są przyzwoite, przynajmniej dla mnie.
Zachęcam do sprawdzania a wszystkie ewentualne poprawki mile widziane.Statystyki: Napisane przez mrek33 — 2 sty 2019, o 15:35
]]>