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$

SUITE DU PROGRAMME

RETOUR A MES TRAVAUX

RETOUR AU MENU

RETOUR PAGE PRECEDENTE