Page copy protected against web site content infringement by Copyscape

N° de dépôt 00037236

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 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 équation du temps sur un cadran vertical à style oiseau

multi indicateur prêt à être posé.

Retour à la page

Retour au menu