Programme informatique de construction
des cadrans solaires bifilaires et à polos
2è partie (2/4)
par
Jean Pakhomoff
gnomoniste à Marseille
1510 ZO$="ENTREE EN CAPRICORNE LE 22/12 (lo=270°:HIVER)":LPRINT:RETURN
1520 ZO$="ENTREE EN VERSEAU LE 20/1 (LO=300°) ET EN SAGITTAIRE LE 22/11 (LO=240
°)":RETURN
1530 ZO$="ENTREE EN POISSONS LE 19/2 (LO=330°) ET EN SCORPION LE 23/10 (LO=210°)
":RETURN
1540 ZO$="ENTREE EN BELIER LE 20/3 (LO=0° PRINTEMPS) ET EN BALANCE LE 23/9 (LO=1
80° AUTOMNE)":RETURN
1550 ZO$="ENTREE EN TAUREAU LE 20/4 (LO=30°) ET EN VIERGE LE 23/8 (LO=150°)":RET
URN
1560 ZO$="ENTREE EN GEMEAUX LE 21/5 (LO=60°) ET EN LION LE 23/7 (LO=120°)":RETUR
N
1570 ZO$="ENTREE EN CANCER LE 21/6 (LO=90° ETE)":RETURN
1580 REM position de A et C dans l'espace.
1590 ABP=AB*COS(K):ACP=AC*COS(K):ACPP=AC*SIN(K)
1600 PRINT #1,"Distance BB'du mur au fil passant par B :";FIX(ABP*100#)/100#;" c
m."
1610 PRINT #1,"Distance CC'du mur au fil passant par C :";FIX(ACP*100#)/100#;" c
m."
1620 PRINT #1,"Ce fil méridien passe dans le plan parallèle à la méridienne et s
itué à"
1630 PRINT #1," AC' = ";FIX(ACPP*100#)/100#;" cm. : vers l'Est si DW ou INW, ver
s l'Ouest si DE ou INE"
1640 RETURN
1650 IF DD=B(4) THEN GOSUB 1510:GOTO 1730
1660 IF DD=B(3) THEN GOSUB 1520:GOTO 1730
1670 IF DD=B(2) THEN GOSUB 1530:GOTO 1730
1680 IF DD=B(1) THEN GOSUB 1540:GOTO 1730
1690 IF DD=DEC(2) THEN GOSUB 1550:GOTO 1730
1700 IF DD=DEC(3) THEN GOSUB 1560:GOTO 1730
1710 IF DD=DEC(4) THEN GOSUB 1570:GOTO 1730
1720 PRINT #1,"Arc diurne de déclinaison ";" ";DD;"°":GOTO 1740
1730 PRINT #1,ZO$
1740 ON N GOTO 2940,3430,4110
1750 REM LH DELE, DWLW sur f
1760 IF F<PI/2#-.0000000001# THEN 1780
1770 RETURN
1780 MM=ATN(1#/(COS(DG)/(COS(F)*TAN(R))-SIN(DG)*TAN(F))):IF MM<0 THEN MM=MM+PI:R
EM
1790 H=MM*180#/PI:RETURN
1800 REM LH INELW, INWLE sur f
1810 IF F<PI/2#-.0000000001# THEN 1830
1820 RETURN
1830 MMM=ATN(-1#/(COS(DG)/(COS(F)*TAN(R))-SIN(DG)*TAN(F))):IF MMM<0 THEN MMM=MMM
+PI
1840 H=MMM*180#/PI:RETURN
1850 REM LH DELW, DWLE sur f
1860 IF F<PI/2#-.0000000001# THEN 1880
1870 RETURN
1880 HH=ATN(1#/(COS(DG)/(COS(F)*TAN(R))+SIN(DG)*TAN(F))):IF HH <0 THEN HH=HH+PI
1890 H=HH*180#/PI:RETURN
1900 REM LH INELE, INWLW sur f
1910 IF F<PI/2#-.0000000001# THEN 1930
1920 RETURN
1930 HHH=ATN(-1#/(COS(DG)/(COS(F)*TAN(R))+SIN(DG)*TAN(F))):IF HHH<0 THEN HHH=HHH
+PI
1940 H=HHH*180#/PI:RETURN
1950 REM Azimut. Si f=pi/2 ====>azimut aa = angle horaire rrr
1960 IF F<.00000000001# THEN IF ABS(D)<.0000000001# THEN A=PI/2#:AA=90#:RETURN
1970 A=ATN(SIN(R)/(COS(R)*SIN(F)-TAN(D)*COS(F))):IF A<0 THEN A=A+PI
1980 AA=A*180#/PI:RETURN
1990 RETURN
2000 REM Hauteur
2010 BB=COS(D)*COS(R)*COS(F)+SIN(D)*SIN(F):BB=ATN(BB/SQR(1#-BB^2#))
2020 HT=BB*180#/PI:RETURN
2030 REM LH HORIZONTAL
2040 IF F>.000000000001# THEN 2060
2050 RETURN
2060 S=ATN(SIN(F)*TAN(R)):IF S<0 THEN S=S+PI
2070 H=S*180#/PI:RETURN
2080 REM X et Y horizontaux
2090 IF BB<.0000001# THEN X=.001#:Y=.001#:HV=RRR:RETURN
2100 X=AC*SIN(A)/TAN(BB)
2110 Y=AB*COS(A)/TAN(BB):YO=Y+LSTC*COS(F)
2115 IF ELISABETHPAKHOMOFF$="p" THEN IF F<.000000001# THEN YO=Y
2120 RETURN
2130 REM X et Y verticaux déclinants et inclinants
2140 IF ABS(A+K-PI/2#)<.0000000001# THEN X=5.540624E+34:GOTO 2160:REM x = infini
car a+k=pi/2 (aw=pi/2+k avec K>0 ou <0 selon le cas : voir programme. Cas des l
evers ou couchers des déclinaisons magnifiques où a=aw).
2150 X=AC*SIN(A)/COS(A+K):XXX=X*AB/AC:REM xxx=AB''
2160 IF ABS(BB)<.000000001# THEN IF ABS(A+K-PI/2#)<.000000001# THEN Y=0:PRINT #1
,"Calcul de y avec impossibilité de forme 0/0 car bb=0 et a+k=90°";" k=";K:RETUR
N
2170 REM avec k >0 ou <0. (voir ligne 2140)
2180 Y=AB*COS(K)*TAN(BB)/COS(A+K)
2190 RETURN
2200 PRINT"Pour l'hémisphère Sud, entrer les déclinaisons X -1 par rapport aux s
ignes zodiacaux. Les LW sont à droite du cadran en regardant le Sud, les LE sont
à gauche (déclinant Sud).":PRINT
2210 PRINT" Inverse en regardant le Nord (inclinant au pôle Sud). Les calculs et
résultats sont identiques.":PRINT
2220 INPUT "Votre latitude f en ° ?",FF:F=FF*PI/180#:LPRINT:IF F>0 THEN IF F<PI/
2#-.0000000001# THEN INPUT "Nord= n, Sud= s ?",GH$