FUNCTION TEST(X,Y)
C DIMENSION QA(2) 
  VECTOR VCAS(1)
C
C X=QA(1)
C Y=QA(2)
C WRITE(6,101) X,Y
101  FORMAT(' X,Y=',2F6.2)
C
C LES 2 FONCTIONS ``qtest.for'' ET ``test.for'' DOIVENT A TOUT 
C MOMENT ETRE QUASI-IDENTIQUES.
C
C POUR PASSER DE ``QTEST'' A ``TEST'' IL Y A 4 CHANGEMENTS:
C   1) QTEST -> TEST
C   2) QA -> X,Y
C   3) C DEVANT ``DIMENSION QA(2)''
C   4) C DEVANT ``X=QA(1)'' ET DEVANT ``Y=QA(2)''
C
C   ``TEST'' EST POUR TRACER UN HISTOGRAMME PAR: FUN2 10 TEST.FOR
C   ``QTEST'' EST POUR LES TESTS DE BALAYAGE DE L'HISTO. ON Y PASSE
C LES ARGUMENTS SOUS FORME VECTORIELLE PUISQU'ON NE
C PEUT PAS LES PASSER SOUS FORME SCALAIRE PAR [X],[Y] (message d'erreur)
C
C La fonction TEST(X,Y) est utilisee dans MAPO#3D (et ANIM#3D) pour
C visualisation 3d par HI/PLOT
C
C CAS=1: 2 GAUSSIENNES de centres (1,2) (-3,-2)-> -6<X<4 -6<Y<4
C CAS=2: TEST DE DIVERSES FONCTIONS ELEMENTAIRES
C CAS=3: ENCEINTE+ABBAYE: 
C             ENCEINTE: -8<X<8 -8<Y<8
C               ABBAYE: 
C CAS=4: EFFETS DE PERSPECTIVE, MAISONS DANS RUE
C CAS=5: CHAPELLE DES PENITENTS: -5<X<5 -10<Y<5
C
 CAS=VCAS(1)
C WRITE(6,201) CAS
201 FORMAT(' CAS=',F6.0)
C
C  1) DEUX GAUSSIENNES
 IF(CAS.NE.1) GOTO 1
 X1=1
 Y1=2
 P1=EXP(-(X-X1)**2-(Y-Y1)**2)
 X2=-3
 Y2=-2
 P2=EXP(-(X-X2)**2-(Y-Y2)**2)
 TEST=P1+3*P2
1 CONTINUE  
C
C 2) TEST DES FONCTIONS HA,HM,HT,HE,PA,PA4,PC
 IF(CAS.NE.2) GOTO 2
C
C TEST DE HA (POUR ABBATIALE)
 GOTO 21
 X1=-5 
 X2= 5
 Y1=-2
 Y2=2
 H=1
 TEST=HA(X,Y,X1,X2,Y1,Y2,H)
C  
C WRITE(6,100) X,Y,TEST
100  FORMAT(' X=',F6.2,'Y=',F6.2,' T=',F3.0)
21 CONTINUE
C
C TEST DE HM (MAISON-PARALLELEPIPEDE)
 GOTO 22
 XZ=5
 YZ=4
 BX=2
 BY=1
 H=1
 TEST=HM(X,Y,XZ,YZ,BX,BY,H)
 WRITE(6,100) X,Y,TEST
22 CONTINUE
C
C TEST DE HT (TOUR-CYLINDRE)
 GOTO 23
 XZ=5
 YZ=4
 R=1
 H=1
 TEST=HT(X,Y,XZ,YZ,R,H)
23 CONTINUE
C
C TEST DE HE (MUR D'ENCEINTE)
 GOTO 24
 XZ=0
 YZ=0
 R1=9
 R2=9.3
 H=1
 TEST=HE(X,Y,XZ,YZ,R1,R2,H)
24 CONTINUE
C
C TEST DE PA (MAISON AVEC TOIT EN PENTE)
C LIGNE DE FAITE SELON OY
 GOTO 25
 A=1
 B=2
 H1=2
 H2=1
 TEST=PA(X,Y,A,B,H1,H2)
25 CONTINUE
C
C TEST DE PA4 (TOUR CARREE AVEC 4 TOITS EN PENTE)
 GOTO 26
 A=0.5
 H1=3
 H2=0.5
 TEST=PA4(X,Y,A,H1,H2)
26 CONTINUE
C
C TEST DE PC (TOUR RONDE AVEC TOIT CONIQUE)
 GOTO 27
 R=0.5 
 H1=3 ; H2=2
 TEST=PC(X,Y,R,H1,H2)
27   CONTINUE
C
C TEST DE PCA8 (INTERSECTION DE 2 TOITS PERP)
  GOTO 28
  A=1 ; AL=3 ; H1=1 ; H2=1.1
 TEST=PCA8(X,Y,A,AL,H1,H2)
28 CONTINUE  
C
C TEST DE PCM (MOITIE DE TOUR RONDE AVEC TOIT CONIQUE)
 GOTO 29
 R=3
 H1=3 ; H2=2
 TEST=PCM(X,Y,R,H1,H2)
29   CONTINUE
C
C TEST DE PAM (MOITIE DE MAISON AVEC 1 TOIT EN PENTE)
C ON NE CONSERVE QUE LA PARTIE X>0
C GOTO 210
 A=2 ; B=3 ; H1=2 ; H2=1
 TEST=PAM(X,Y,A,B,H1,H2)
210  CONTINUE
C
2  CONTINUE  
CCCC
C 3) REPRESENTATION DE BEAULIEU
C  REMARQUE: LES DETAILS (COMME LES TOURS DE GUET)
C   N'APPARAISSENT QUE SI ON ZOOME DESSUS, AUTREMENT
C   IL FAUDRAIT
C UN DECOUPAGE TRES FIN POUR LES RENDRE VISIBLES
C (MAIS UN TEL DECOUPAGE NORCIT L'IMAGE)
 IF(CAS.NE.3) GOTO 3
C
C     MUR D'ENCEINTE + CHEMIN DE RONDE DERRIERE MUR
 XZ=0 ; YZ=0
 RA1=7.8 ; RA2=8.
 H=1.
 E1=HE(X,Y,XZ,YZ,RA1,RA2,H)
 RB1=7.4 ; RB2=RA1
 H=0.5
 E2=HE(X,Y,XZ,YZ,RB1,RB2,H)
 E=E1+E2
C
C     4 TOURS CARRES DEVANT LE REMPART
C     COTE DU CARRE=2A
C  5 MARS 2013: ON A MIS LES MURS // AU REMPART
 GOTO 310
 XZ=6
 YZ=SQRT(R1**2-XZ**2)
 BX=0.3
 BY=BX
 H=2
 T1=HM(X,Y,XZ,YZ,BX,BY,H)
 XZ=4
 YZ=-SQRT(R1**2-XZ**2)
 T2=HM(X,Y,XZ,YZ,BX,BY,H)
 XZ=-5
 YZ=-SQRT(R1**2-XZ**2)
 T3=HM(X,Y,XZ,YZ,BX,BY,H)
 XZ=-3
 YZ= SQRT(R1**2-XZ**2)
 T4=HM(X,Y,XZ,YZ,BX,BY,H)
 T=T1+T2+T3+T4
310  CONTINUE
 BX=0.4 ; BY=BX
 XZ=8+BX/2 ; YZ=0 ; H=2
C T1=HM(X,Y,XZ,YZ,BX,BY,H)
C
 A=0.1 ; H1=2.5 ; H2=0.8
C
 XZ=8+A ; YZ=0
C
C T1
C COORDONNEES (XRZ,YRZ) DE L'EMPLACEMENT CHOISI
 AL=40 ; AL=AL*3.14159/180 
 CO=COS(AL) ; SI=SIN(AL)
 XRZ=XZ*CO ; YRZ=XZ*SI
C ON CENTRE LE REPERE SUR CE POINT
 XX=X-XRZ ; YY=Y-YRZ
C ON FAIT UNE ROTATION DU MEME ANGLE AL
 XX2=XX*CO+YY*SI ; YY2=XX*SI-YY*CO
 T1=PA4(XX2,YY2,A,H1,H2)
C
C T2
C COORDONNEES (XRZ,YRZ) DE L'EMPLACEMENT CHOISI
 AL=110 ; AL=AL*3.14159/180 
 CO=COS(AL) ; SI=SIN(AL)
 XRZ=XZ*CO ; YRZ=XZ*SI
C ON CENTRE LE REPERE SUR CE POINT
 XX=X-XRZ ; YY=Y-YRZ
C ON FAIT UNE ROTATION DU MEME ANGLE AL
 XX2=XX*CO+YY*SI ; YY2=XX*SI-YY*CO
 T2=PA4(XX2,YY2,A,H1,H2)
C
C T3
C COORDONNEES (XRZ,YRZ) DE L'EMPLACEMENT CHOISI
 AL=230 ; AL=AL*3.14159/180 
 CO=COS(AL) ; SI=SIN(AL)
 XRZ=XZ*CO ; YRZ=XZ*SI
C ON CENTRE LE REPERE SUR CE POINT
 XX=X-XRZ ; YY=Y-YRZ
C ON FAIT UNE ROTATION DU MEME ANGLE AL
 XX2=XX*CO+YY*SI ; YY2=XX*SI-YY*CO
 T3=PA4(XX2,YY2,A,H1,H2)
C
C T4
C COORDONNEES (XRZ,YRZ) DE L'EMPLACEMENT CHOISI
 AL=300 ; AL=AL*3.14159/180 
 CO=COS(AL) ; SI=SIN(AL)
 XRZ=XZ*CO ; YRZ=XZ*SI
C ON CENTRE LE REPERE SUR CE POINT
 XX=X-XRZ ; YY=Y-YRZ
C ON FAIT UNE ROTATION DU MEME ANGLE AL
 XX2=XX*CO+YY*SI ; YY2=XX*SI-YY*CO
 T4=PA4(XX2,YY2,A,H1,H2)
C
 T=T1+T2+T3+T4
C
C         ABBATIALE
C FACTEUR D'AGRANDISSEMENT G -> LORSQUE J'AI MIS
C L'ABBATIALE SUR LE MEME DESSIN QUE LA CHAPELLE DES PENITENTS
C JE ME SUIS APERCU QUE POUR LA CHAPELLE J'AVAIS PRIS DES
C DIMENSIONS BIEN PLUS GRANDES QUE CELLES DE L'ABBATIALE
 G=2
C
C    CORPS PRINCIPAL
C LA BASE EST (BX)x(BY) SOIT (2*G)x(0.8*G)
 XZ=0*G
 YZ=0*G
 BX=2*G
 BY=0.8*G
 H=1.6*G
 A1=HM(X,Y,XZ,YZ,BX,BY,H)
C    PARTIE ARRONDIES COTE CHOEUR
 XZ=1.2*G
 YZ=0*G
 BX=0.4*G
 BY=0.4*G
 H=1.1*G
 A2=HM(X,Y,XZ,YZ,BX,BY,H)
C    TOUR CARREE COTE OUEST
 XZ=-0.8*G
 YZ=-0.2*G
 BX=0.4*G
 BY=0.4*G
 H=2.5*G
 A3=HM(X,Y,XZ,YZ,BX,BY,H)
C   MINCE TOURELLE SURMONTANT LA TOUR
 XZ=-0.8*G
 YZ=-0.2*G
 BX=0.12*G
 BY=0.12*G
 H=0.7*G
 A4=HM(X,Y,XZ,YZ,BX,BY,H)
C    TOUR RONDE SUR LE TOIT COTE CHOEUR
 XZ=0.6*G
 YZ=0*G
 R=0.20*G
 H=1.5*G
 A5=HT(X,Y,XZ,YZ,R,H)
C     CLOITRE (MAINTENANT DETRUIT) COTE NORD
 XZ=0.5*G
 YZ=0.9*G
 BX=1*G
 BY=1*G
 H=0.8*G
 A6A=HM(X,Y,XZ,YZ,BX,BY,H)
 BX=0.6*G
 BY=0.6*G
 A6B=HM(X,Y,XZ,YZ,BX,BY,H)
 A6=A6A-A6B
C           TOTAL
 A=A1+A2+A3+A4+A5+A6
C
 TEST=E+T+A
C WRITE(6,200) E,TEST
200  FORMAT(' E,TEST=',2F5.1)
3 CONTINUE
CCCC
C 4) EFFETS DE PERSPECTIVE
C    MAISONS RECTANGULAIRES LE LONG D'UNE RUE
 IF(CAS.NE.4) GOTO 4
 XZ=2
 BX=1
 YZ=2.5
 BY=2.
 H=2
 T1=HM(X,Y,XZ,YZ,BX,BY,H)
 H=3
 YZ=7.5
 T2=HM(X,Y,XZ,YZ,BX,BY,H)
C
C GOTO 41
 XZ=8
 YZ=5
 XX=X-XZ
 YY=Y-YZ
 A=1.1
 B=2.5
 H1=2
 H2=1
 T3=PA(XX,YY,A,B,H1,H2)
41  CONTINUE
C
 TEST=T1+T2+T3
4 CONTINUE
C
CCCC
C 5) CHAPELLE DES PENITENTS: SURFACE AU SOL-> -5<X<5 -10<Y<5
 IF(CAS.NE.5) GOTO 5
C
C POUR L'INTEGRATION DANS LE PLAN DE BEAULIEU
C L'EMPLACEMENT DE LA CHAPELLE EST (XC,YC)
 XC=11.5 ; YC=2
C POUR ZOOMER SUR LA CHAPELLE DESACTIVER LA LIGNE SUIVANTE
C X=X-XC ; Y=Y-YC
C
C TOITS FORMANT CROIX
* LE TOUT EST COMPRIS DS UN CARRE DE COTE 2*AL ET CENTRE (0,0),
* DONC COUVRANT -2<X<2, -2<Y<2
* LA LARGEUR DE CHAQUE BRANCHE EST 2*AA
* LES H SONT LES HAUTEURS
 AA=1 ; AL=2 ; HH1=2 ; HH2=1
 T1=PCA8(X,Y,AA,AL,HH1,HH2)
C PROLONGEMENT DE LA CROIX DS DIRECTION OY 
C PA() FAIT UN TOIT -B2<YY<B2, SOIT -B2+YZ<Y<B2+YZ, SOIT:
C -AL-4*AA<Y<-AL, SOIT ICI -6<Y<-2
 A2=AA ; B2=2*A2 
 XZ=0 ; YZ=-AL-B2
 XX=X-XZ ; YY=Y-YZ
 T2=PA(XX,YY,A2,B2,HH1,HH2)
C 
C CLOCHER
 A3=0.6*AA ; B3=0.4
 YZ=YZ-B2-B3
 XX=X-XZ ; YY=Y-YZ
 H1=4 ; H2=1
 T3=PA(XX,YY,A3,B3,H1,H2)
C 
C DEUX CONTREFORTS DE CHAQUE COTE DU CLOCHER
 A4=0.3*A3 ; B4=B3 ; H1=2.6 ; H2=0.4
 XZ=A3
 XX=X-XZ ; YY=Y-YZ
 T4=PAM(XX,YY,A4,B4,H1,H2)
 XZ=-XZ
 XX=X-XZ ; YY=Y-YZ
 T4B=PAMB(XX,YY,A4,B4,H1,H2)
 T4=T4+T4B
C
C MOITIE DE TOUR A L'OPPOSE DU CLOCHER
 R=AA 
 XZ=0 ; YZ=AL
 XX=X-XZ ; YY=Y-YZ
 T5=PCM(XX,YY,R,HH1,HH2)
C
C MAISON SITUEE A G DU CLOCHER
 A6=0.8 ; H1=1.3 ; H2=1
 XZ=-3 ; YZ=-9 
 XX=X-XZ ; YY=Y-YZ
 T6=PA4(XX,YY,A6,H1,H2)
C 
C TOTAL
 TEST=T1+T2+T3+T4+T5+T6
5  CONTINUE
C
 QTEST=TEST
C
C WRITE(6,202) CAS,X,Y,TEST
202  FORMAT(' CAS=',F10.3,' X,Y=',2F10.2,' TEST=',F10.3)
C
 END
CCCCCCCCCCCCCCCCCCCCCCCC
 FUNCTION HA(X,Y,X1,X2,Y1,Y2,H)
C LA FCT VAUT H DS LE RECTANGLE X1,X2,Y1,Y2 
C ET 0 AILLEURS
C POUR AVOIR H IL FAUT DONC QUE:
C X1<X<X2 ET Y1<Y<Y2
 PX=(X-X1)*(X-X2)
 PY=(Y-Y1)*(Y-Y2)
 IF(PX.GT.0.) TX=0
 IF(PX.LT.0.) TX=1
 IF(PY.GT.0.) TY=0
 IF(PY.LT.0.) TY=1
 A=TX*TY
 IF(A.EQ.0) HA=0
 IF(A.EQ.1) HA=H
C WRITE(6,100) X,Y,X1,X2,Y1,Y2,HA
100  FORMAT(' X,Y=',2F4.1,' X12,Y12=',4F4.1,' FA=',F3.0)
 END
CCCCCCCCCCCCCCCCCCCCCCCC
 FUNCTION HM(X,Y,XZ,YZ,BX,BY,H)
C FCT POUR REPRESENTER UNE MAISON PAR
C UN PARALLELEPIPEDE DE CENTRE (XZ,YZ),
C DE BASE BXxBY ET DE HAUTEUR H
 X1=XZ-BX/2 
 X2=XZ+BX/2
 Y1=YZ-BY/2
 Y2=YZ+BY/2
 HM=HA(X,Y,X1,X2,Y1,Y2,H)
 END
CCCCCCCCCCCCCCCCCCCCCCCCC
 FUNCTION HT(X,Y,XZ,YZ,R,H)
C FCT POUR REPRESENTER UNE TOUR PAR UN
C CYLINDRE DE CENTRE (XZ,YZ), DE RAYON R
C ET DE HAUTEUR H
 C=(X-XZ)**2+(Y-YZ)**2-R**2
 IF(C.LE.0) HT=H
 IF(C.GT.0) HT=0
 END
CCCCCCCCCCCCCCCCCCCCCCCCC
 FUNCTION HE(X,Y,XZ,YZ,R1,R2,H)
C FCT POUR REPRESENTER LE MUR D'ENCEINTE PAR
C UN ANNEAU CIRCULAIRE DE CENTRE (XZ,YZ), DE
C RAYON INTERIEUR R1, DE RAYON EXTERIEUR R2 
C ET DE HAUTEUR H
 C1=(X-XZ)**2+(Y-YZ)**2-R1**2
 C2=(X-XZ)**2+(Y-YZ)**2-R2**2
C
C IL FAUT QUE C1>0 ET C2<0
 IF(C1.LE.0) E1=0
 IF(C1.GT.0) E1=1
 IF(C2.LE.0) E2=1
 IF(C2.GT.0) E2=0
 A=E1*E2
 IF(A.EQ.0) HE=0
 IF(A.EQ.1) THEN
 HE=H
C WRITE(6,100) X,Y,HE
100   FORMAT(' X,Y=',2(F5.1),' HE=',F5.1)
 ENDIF
 END
CCCCCCCCCCCCCCCCCCCCCCCCC
 FUNCTION PA(X,Y,A,B,H1,H2)
C MAISON AVEC TOIT EN PENTE
C LIGNE DE FAITE SELON OY
C SECTION AU SOL: RECTANGLE DE CENTRE (0,0)
C LARGEUR SELON OX:  -A,A
C LONGUEUR SELON OY: -B,B
C HAUTEUR DU PARALLELEPIPEDE: H1
C HAUTEUR DU TOIT SOUS FAITE: H2
C
C PENTE DU TOIT
C P=H2/(A/SQRT(2.))
  P=H2/A
C POUR AVOIR UNE HAUTEUR H1 IL FAUT QUE:
C -A<X<A ET -B<Y<B
 PX=(X+A)*(X-A)
 PY=(Y+B)*(Y-B)
 IF(PX.GT.0.) TX=0
 IF(PX.LT.0.) TX=1
 IF(PY.GT.0.) TY=0
 IF(PY.LT.0.) TY=1
 QT=TX*TY
 IF(QT.EQ.0) PA=0
C
 IF(QT.EQ.1) THEN
  IF(X.GE.0) PA=-P*X+H2
  IF(X.LT.0) PA= P*X+H2
 PA=PA+H1
 ENDIF
C
 END
CCCCCCCCCCCCCCCCCCCCCCCCC
 FUNCTION PAM(X,Y,A,B,H1,H2)
C MOITIE DE MAISON AVEC TOIT EN PENTE
C ON CONSERVE LA PARTIE X>0
C LIGNE DE FAITE SELON OY
C SECTION AU SOL: RECTANGLE DE CENTRE (0,0)
C LARGEUR SELON OX:  -A,A
C LONGUEUR SELON OY: -B,B
C HAUTEUR DU PARALLELEPIPEDE: H1
C HAUTEUR DU TOIT SOUS FAITE: H2
C
C PENTE DU TOIT
 P=H2/(A/SQRT(2.))
C POUR AVOIR UNE HAUTEUR H1 IL FAUT QUE:
C -A<X<A ET -B<Y<B
 PX=(X+A)*(X-A)
 PY=(Y+B)*(Y-B)
 IF(PX.GT.0.) TX=0
 IF(PX.LT.0.) TX=1
 IF(PY.GT.0.) TY=0
 IF(PY.LT.0.) TY=1
 QT=TX*TY
 IF(QT.EQ.0) PA=0
C
 IF(QT.EQ.1) THEN
  IF(X.GE.0) PA=-P*X+H2
  IF(X.LT.0) PA= P*X+H2
 PA=PA+H1
 ENDIF
C
 IF(X.LE.0) PA=0
C
 PAM=PA
 END
CCCCCCCCCCCCCCCCCCCCCCCCC
 FUNCTION PAMB(X,Y,A,B,H1,H2)
C MOITIE DE MAISON AVEC TOIT EN PENTE
C ON CONSERVE LA PARTIE X<0
C LIGNE DE FAITE SELON OY
C SECTION AU SOL: RECTANGLE DE CENTRE (0,0)
C LARGEUR SELON OX:  -A,A
C LONGUEUR SELON OY: -B,B
C HAUTEUR DU PARALLELEPIPEDE: H1
C HAUTEUR DU TOIT SOUS FAITE: H2
C
C PENTE DU TOIT
 P=H2/(A/SQRT(2.))
C POUR AVOIR UNE HAUTEUR H1 IL FAUT QUE:
C -A<X<A ET -B<Y<B
 PX=(X+A)*(X-A)
 PY=(Y+B)*(Y-B)
 IF(PX.GT.0.) TX=0
 IF(PX.LT.0.) TX=1
 IF(PY.GT.0.) TY=0
 IF(PY.LT.0.) TY=1
 QT=TX*TY
 IF(QT.EQ.0) PA=0
C
 IF(QT.EQ.1) THEN
  IF(X.GE.0) PA=-P*X+H2
  IF(X.LT.0) PA= P*X+H2
 PA=PA+H1
 ENDIF
C
 IF(X.GE.0) PA=0
C
 PAMB=PA
 END
CCCCCCCCCCCCCCCCCCCC
 FUNCTION PA4(X,Y,A,H1,H2)
C TOUR OU MAISON CARREE SURMONTEE D'UN TOIT A 4 PANS
C SECTION AU SOL: CARRE DE CENTRE (0,0) ET COTE 2A
C HAUTEUR DE LA TOUR SANS LE TOIT: H1
C HAUTEUR DU TOIT: H2
C VERSION AVEC MURS PARALLELES A OX,OY
C
C PENTE DU TOIT
 P=H2/A
C
 PA=0
 XA=ABS(X) ; YA=ABS(Y)
 TX=0 ; TY=0
 IF(XA.LE.A) TX=1 ; IF(YA.LE.A) TY=1
 QT=TX*TY
C
 IF(QT.EQ.1) THEN
C
C QUADRANT DE DROITE
 T1=0 ; T2=0
 IF(X.GE.0) T1=1 ; IF(XA.GE.YA) T2=1
 Q=T1*T2
 IF(Q.EQ.1) PA=H1+H2-P*XA
C QUADRANT VERS LE HAUT
 T1=0 ; T2=0
 IF(Y.GE.0) T1=1 ; IF(XA.LE.YA) T2=1
 Q=T1*T2
 IF(Q.EQ.1) PA=H1+H2-P*YA
C QUADRANT DE GAUCHE
 T1=0 ; T2=0
 IF(X.LE.0) T1=1 ; IF(XA.GE.YA) T2=1
 Q=T1*T2
 IF(Q.EQ.1) PA=H1+H2-P*XA
C QUADRANT VERS LE BAS
 T1=0 ; T2=0
 IF(Y.LE.0) T1=1 ; IF(XA.LE.YA) T2=1
 Q=T1*T2
 IF(Q.EQ.1) PA=H1+H2-P*YA
C
 ENDIF
 PA4=PA
 END
CCCCCCCCCCCCCCCCCCCC
 FUNCTION PCA8(X,Y,A,AL,H1,H2)
C EDIFICE EN FORME DE CROIX SYMETRIQUE
C AVEC TOITS EN V INVERSE SUR CHAQUE BRANCHE.
C AVEC DONC INTERSECTION DE 2 TOITS PERPEND.
C CELA DONNE UN TOIT A 8 PANS
C SECTION AU SOL: CROIX DE CENTRE (0,0) 
C POUR LES 4 BRANCHES DE LA CROIX:
C LARGEUR 2A, LONGUEUR DES BRANCHES 2*AL 
C HAUTEUR SANS LE TOIT: H1
C HAUTEUR DU TOIT: H2
C
C CETTE CROIX POURRA ENSUITE ETRE PROLONGEE
C SUR CHAQUE DIRECTION POUR FAIRE UNE CROIX AVEC
C UNE BRANCHE PLUS LONGUE QUE L'AUTRE COMME
C LA CHAPELLE DES PENITENTS
C
C ON VA TRAITER LES 8 PANS SUCCESSIVEMENT EN COMMENCANT
C PAR LE PAN (X>0,Y>0,X>Y), PUIS EN TOURNANT DANS LE SENS
C TRIGONOMETRIQUE
C
C ON SUPPOSE QUE L>A
 IF(AL.LT.A) THEN 
 WRITE(6,200) AL,A
200  FORMAT(' PCA8: AL<A, AL,A=',2(F7.1))
 PC=0
 GOTO 1
 ENDIF 
C
C PENTE DU TOIT 
 P=H2/A
C
C TEST POUR LE CARRE (-AL,AL)
 PX=(X+AL)*(X-AL)
 PY=(Y+AL)*(Y-AL)
 TX=0
 TY=0
 IF(PX.LT.0.) TX=1
 IF(PY.LT.0.) TY=1
 QT=TX*TY
C ON MET A 0 HORS DU CARRE (-AL,AL)
 IF(QT.EQ.0) THEN
 PC=0
 GOTO 1
 ENDIF
C
C  PANS 1X,1Y
 TX=0 ; TY=0
 IF(X.GE.0) TX=1 ; IF(Y.GE.0) TY=1
 QT=TX*TY
 IF(QT.EQ.1) THEN
C     
 XA=ABS(X) ; YA=ABS(Y)
C      PAN 1X, CAD ABS(X)>ABS(Y)
 IF(XA.GE.YA) THEN 
 PC=0
 IF(YA.LE.A) THEN ; PC=H1+H2-P*YA ; ENDIF
 ENDIF 
C      PAN 1Y, CAD ABS(X)<ABS(Y)
 IF(XA.LE.YA) THEN 
 PC=0
 IF(XA.LE.A) THEN ; PC=H1+H2-P*XA ; ENDIF
 ENDIF
C
 ENDIF
C
C PANS 2X,2Y
 TX=0 ; TY=0
 IF(X.LE.0) TX=1 ; IF(Y.GE.0) TY=1
 QT=TX*TY
 IF(QT.EQ.1) THEN
C
 XA=ABS(X) ; YA=ABS(Y)
C      PAN 2X
 IF(XA.GE.YA) THEN 
 PC=0
 IF(YA.LE.A) THEN ; PC=H1+H2-P*YA ; ENDIF
 ENDIF 
C      PAN 2Y
 IF(XA.LE.YA) THEN 
 PC=0
 IF(XA.LE.A) THEN ; PC=H1+H2-P*XA ; ENDIF
 ENDIF
C
 ENDIF
C
C PANS 3X,3Y
 TX=0 ; TY=0
 IF(X.LE.0) TX=1 ; IF(Y.LE.0) TY=1
 QT=TX*TY
 IF(QT.EQ.1) THEN
C
 XA=ABS(X) ; YA=ABS(Y)
C      PAN 3X
 IF(XA.GE.YA) THEN 
 PC=0
 IF(YA.LE.A) THEN ;  PC=H1+H2-P*YA ; ENDIF
 ENDIF 
C      PAN 3Y
 IF(XA.LE.YA) THEN 
 PC=0
 IF(XA.LE.A) THEN ; PC=H1+H2-P*XA ; ENDIF
 ENDIF
C
 ENDIF
C
C PANS 4X,4Y
 TX=0 ; TY=0
 IF(X.GE.0) TX=1 ; IF(Y.LE.0) TY=1
 QT=TX*TY
 IF(QT.EQ.1) THEN
C
 XA=ABS(X) ; YA=ABS(Y)
C      PAN 4X
 IF(XA.GE.YA) THEN 
 PC=0
 IF(YA.LE.A) THEN ;  PC=H1+H2-P*YA ; ENDIF
 ENDIF 
C      PAN 4Y
 IF(XA.LE.YA) THEN 
 PC=0
 IF(XA.LE.A) THEN ; PC=H1+H2-P*XA ; ENDIF
 ENDIF
C
 ENDIF
C
1 CONTINUE
 PCA8=PC
C
 END
CCCCCCCCCCCCCCCCCCCCCCCCC
 FUNCTION PC(X,Y,R,H1,H2)
C TOUR CENTREE EN (0,0), DE RAYON R
C SURMONTE D'UN TOIT CONIQUE
C  H1=HAUTEUR DU CYLINDRE
C  H2=HAUTEUR DU TOIT CONIQUE
C VOIR FORMULE ET DESSIN DS SCHAUM P. 51
 XZ=0
 YZ=0
 C=(X-XZ)**2+(Y-YZ)**2-R**2
C
 IF(C.GT.0) PC=0
C
 IF(C.LE.0) THEN
 PC=H1+H2-(H2/R)*SQRT(X**2+Y**2)
 ENDIF
C
 END
CCCCCCCCCCCCCCCCCCCCCCCCC
 FUNCTION PCM(X,Y,R,H1,H2)
C MOITIE (Y>0) DE TOUR CENTREE EN (0,0), DE RAYON R
C SURMONTE D'UN TOIT CONIQUE
C  H1=HAUTEUR DU CYLINDRE
C  H2=HAUTEUR DU TOIT CONIQUE
C VOIR FORMULE ET DESSIN DS SCHAUM P. 51
 XZ=0
 YZ=0
 C=(X-XZ)**2+(Y-YZ)**2-R**2
C
 IF(C.GT.0) PC=0
C
 IF(C.LE.0) THEN
 PC=H1+H2-(H2/R)*SQRT(X**2+Y**2)
 ENDIF
C
 IF(Y.LE.0) PC=0
C
 PCM=PC
 END