programmes gwbasic donnant la valeur de l'équation du temps et la
transformation de l'heure lue au cadran
en heure civile (et réciproquement).
jean pakhomoff
20 rem EQTPER 30 3 2020
30 REM VALEUR MOYENNE DE L'EQUATION DU TEMPS somme de l'équation du centre
et de la réduction à l'équateur. 460" de temps=1°55'=1.916666667°
40 REM Le moyen mouvement=n*t avec n=360/365.242199 et t=temps écoulé en jours depuis
le passage au périgé le Zè jour.
45 REM L=(281.2187#+A*61.717#/3600)#*PI/180# + nt + (1.916666667#*PI/180#)*sin nt
50 REM 148' d'angle = 592" de temps. La formule de l'équation du
temps est EQT=C+R=460*sin(n*t)-592*sin(2*L)
60 DEFDBL A-Z:KEY OFF:CLS:PI=4#*ATN(1#)
70 OPEN "EQcal.doc" FOR OUTPUT AS #1
80 INPUT"Numéro de l'année à partir de 1900: 1920, 2050...?",A
90 A=A-1900
92 Z=(61.717#*A)/3600#: REM déplacement du périhélie
94 ZZ=Z/.9856473321#: REM nombre de jours correspondants
100 INPUT "Date: jour et mois: 1 à 31 et 1 à 12: a,b ?",J,N
110 IF N=1 THEN M=0:GOTO 230
120 IF N=2 THEN M=31:GOTO 230
130 IF N=3 THEN M=59:GOTO 230
140 IF N=4 THEN M=90:GOTO 230
150 IF N=5 THEN M=120:GOTO 230
160 IF N=6 THEN M=151:GOTO 230
170 IF N=7 THEN M=181:GOTO 230
180 IF N=8 THEN M=212:GOTO 230
190 IF N=9 THEN M=243:GOTO 230
200 IF N=10 THEN M=273:GOTO 230
210 IF N=11 THEN M=304:GOTO 230
220 IF N=12 THEN M=334
230 M=M+J
232 Q=M-ZZ: IF Q<0 THEN Q=365-ZZ+M
240 MM=.9856473321#*Q*PI/180#
245 LA=(281.2187#+A*61.717#/3600#)*PI/180#
247 IF LA>2#*PI THEN LA = LA-2#*PI
250 EQT=460#*SIN(MM)-592#*SIN(2#*(LA+MM+(1.916666667#*PI/180#)*SIN(MM)))
255 REM eqt en " transformé en' en divisant par 60.
260 EQT=FIX(EQT*1000#/6#)/10000#
270 PRINT:PRINT"Le ";J;"/";N;" la valeur de l'équation du temps est de ";EQT;"
280 PRINT # 1, J;"/";N;" LA VALEUR DE L'équation du temps est de ";EQT;" '"
282 INPUT "Recommencer, Arrêt : r ou a", Z$
284 IF Z$="r" THEN 290
286 CLOSE:END
290 GOTO 100
*****************************************************************************
programme transformant lheure lue au cadran en heure civile et inversement:
80 REM HSOLCIVF = HSOLCIV adapté à tous les fuseaux 10 12 2007
90 REM HSOLCIV = HCADHCIV AVEC LA RECIPROQUE 22 3 2002
100 REM HCADHCIV donne l'heure civile correspondante à l'heure solaire lue sur le cadran
110 KEY OFF:CLS:DEFDBL A-Z
111 INPUT "heure cadran => heure civile = 1; heure civile => heure cadran = 2 :1 ou 2 ?",N
120 INPUT "Longitude : - pour l'est, + pour l'ouest en ° décimaux ?",L
122 INPUT "LONGITUDE DU FUSEAU par rapport à Greenwich: 0,+15,+30... à L'W; -15,-30...à l'Est", LF
123 DL=L-LF
125 DL=4#*DL: REM transformation en minutes de temps
130 INPUT "Equation du temps en ' et '': +a,+b en Novembre ; -a,-b en Février par exemple",A,B
140 EQ=(60#*A+B)/60#
150 INPUT "heure d'hiver = h; été = e; pas de correction=s",A$
160 IF A$="h" THEN 190
165 IF A$="s" THEN Q=0:GOTO 191
170 Q=120#
180 GOTO 191
190 Q=60#
191 IF N=1 THEN 200
192 G$="heure du cadran":INPUT "heure de la montre en heures et minutes : h,m",H,M
193 H=60#*H+M
194 H=H-Q-DL+EQ
195 GOTO 230
200 G$="heure de la montre":INPUT "Heure solaire du cadran en heures et minutes
210 H=60#*H+M
220 H=H+DL-EQ+Q
230 H=H/60#
235 REM la ligne 240 correspond au programme dechexa
240 DD=H-FIX(H):R=DD*60#:MM=R-FIX(R):S=MM*60#
250 PRINT G$;" = ";FIX(H);"H ";FIX(R);"M ";FIX(S);"S"
**********************************************************************************************