programmes donnant la valeur de l'équation du temps et la
transformation de l'heure lue au cadran
en heure civile (et réciproquement).
30 REM VALEUR MOYENNE DE L'EQUATION DU TEMPS somme de l'équation du centre
et de la réduction à l'équateur. 115' d'angle=460" de temps=1°55'=1.916666667°
40 REM Le moyen mouvement=n*t n=360/365.242199 et t=temps écoulé en jours depuis
le passage au périgé le 2 Janvier en moyenne.
50 REM 148'=592" de temps. lo=la+n*t+115'*sin(n*t) (voir valeur de "la" à la ligne 103).la formule de l'équation du temps est EQT=C+R=460*sin(n*t)-592*sin(2*lo)
60 DEFDBL A-Z:KEY OFF:CLS:PI=4#*ATN(1#)
70 OPEN "EQcal.doc" FOR OUTPUT AS #30
100 INPUT "Date: jour et mois: 1 à 31 et 1 à 12: a,b ?",J,N
102 INPUT " Année ?”, an
103 la=282.93306# + (62#/3600#)*(an-2000#)
104 rem la longitude du périhélie était de 282.93306° le 1/1/2000 et
s’accroît de 62’’ d’angle chaque année. j'appelle donc "la" cette longitude actualisée à l'année an.
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
225 REM m=-1<=>1er Janvier
230 M=M+J-2:IF M=-1 THEN M=364
240 MM=.9856473321#*M*PI/180#
250 EQT=460#*SIN(MM)-592#*SIN(2#*(la*PI/180#+MM+(1.916666667#*PI/180#)*SI
N(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 # 30, 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 l’heure 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"
**********************************************************************************************
une éqution du temps sur un cadran vertical à style oiseau
multi indicateur prêt à être posé.