programmes 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
10 REM EQT2015S 1 7 2015_2
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=0.9856473348 et t=temps écoulé en jours
depuis le passage au périgée le 2 Janvier en moyenne.
50 REM 148'=592" de temps. L0=282°56'+n*t+115'*sin(n*t). la formule de l'équation
du temps est EQT=C+R=460*sin(n*t)-592*sin(2*L0) 282°56'=282,933° Lo le 1 1 2000
60 DEFDBL A-Z:KEY OFF:CLS:PI=4#*ATN(1#)
70 OPEN "EQcal.doc" FOR OUTPUT AS #1
100 INPUT "Année ?",AN
102 INPUT "Date: jour et mois: 1 à 31 et 1 à 12: a,b ?",J,N
103 LA=282.93306# + (61.717#/3600#)*(AN-2000#):LA=LA*pi/180
104 LAA=282.93306# + (61.717#/3600#)*(AN-1#-2000#):LAA=LAA*pi/180
105 REM LA : A LA s'ajoutent les valeurs journalières du moyen mouvement: exemple le 3 1 2015 on part ave la valeur de LA du 2 1 2015:
(282.93306# + (61.717#/3600#)*(2015-2000)) qui est la longitude du périgée au 1 1 2000 à 61.717#/3600#/365 degré près ou 0.0171416#/365 radian près
106 REM LAA correspond aux 1er et 2 janvier de la fin de l'AN-1 : en 2015 :
(282.93306# + (61.717#/3600#)*(2014-2000)) voir les lignes 110 à 240
107 REM la longitude du périhélie était de 282.93306° le 1/1/2000 et s'accroît
de 61.717'' d'angle chaque année. J'appelle donc "LA" cette longitude actualisée
110REM le moyen mouvement journalier correspond à 0.98564° d'écliptique et le mouvement annuel du périgée dans le même sens correspond à 61.717'' ou 0.01714°. Donc
112 REM pour parcourir ces 0.98564° (1 jour ) il faudra au périgée 0,98564/0.01714 = 57ans et demi.
On pourra donc pendant plusieurs années conserver le 2 janvier comme origine des longitudes.
115 REM 1er janvier 1900, 2 janvier 2000
118 IF N=1 THEN M=0:GOTO 230:REM janvier
120 IF N=2 THEN M=31:GOTO 230: REM février
130 IF N=3 THEN M=59:GOTO 230: REM mars...etc...
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 M=364 de l'AN-1
226 REM M=0 <=>2 janvier alors M=365 de l'AN-1
230 M=M+J-2:IF M=-1 THEN M=364#:LA=LAA:GOTO 240
235 IF M=0 THEN M=365#:LA=LAA
240 MM=.98564#*M*PI/180#
243 REM lo=la+nt+115' sin nt<=>lo =la+nt+1.916666667*pi/180 sin nt 115' d'angle =.033452rd
247 LO= LA+MM+.033452#*SIN(MM)
250 EQT=460#*SIN(MM)-592#*SIN(2#*LO)
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;" ' ":PRINT #1,
282 INPUT "Recommencer, Arrêt : r ou a", Z$
284 IF Z$="r" THEN 290
286 CLOSE:END
290 GOTO 102
***************************************************************************************
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"
**********************************************************************************************
une équation du temps sur un cadran vertical à style oiseau
multi indicateur prêt à être posé.