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;"

'":PRINT

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 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"

**********************************************************************************************

Retour à la page

Retour au menu