C [DEM10] C DEMONSTRATION FOR THE GRAPHIC OF A TRIPLE OF SURFACES, C MADE UP OUT POF PLANES AND CONES. C [14-APR-74] C C [THD] EDER 30-NOV-76, COPIED FROM DECUS 10-228 C EXTERNAL PLTCA DIMENSION S1(51,51),S2(51,51),S3(51,51) AS(I,J)=0.02*FLOAT(MOD(I+J,2)) NX=51 NY=51 I0=NX/2 J0=NY/2 C1=2. C2=0.5 C3=0.75 EF=1.025 DO 10 I=1,NX DO 10 J=I,NY X=0.1*FLOAT(I-I0) Y=01.*FLOAT(J-J0) RR1=(X+1.0)**2+(Y-1.2)**2 RR2=(X-1.1)**2+(Y+0.25)**2 RR3=(X-1.7)**2+(Y-1.05)**2 R1=SQRT(RR1) R2=SQRT(RR2) R3=SQRT(RR3) IF(R1.LE.C1) P1=3.5*(1.0-R1/C1) IF((R1.GT.EF*C1).AND.(R1.LT.1.66*C1)) P1=3.5*(1.-EF) IF(R1.GT.1.66*C1) P1=AS(I,J) IF(R2.LE.EF*C2) P2=3.0*(1.0-R2/C2) IF((R2.GT.EF*C2).AND.(R2.LT.2.0*C2)) P2=3.0*(1.0-EF) IF(R2.GT.2.0*C2) P2=AS(I,J) IF(R3.LE.EF*C3) P3=4.0*(1.0-R3/C3) IF((R3.GT.EF*C3).AND.(R3.LT.2.0*C3)) P3=4.0*(1.0-EF) IF(R3.GT.2.0*C3) P3=AS(I,J) S1(I,J)=TANH(0.25*P1-0.3) S2(I,J)=TANH(-0.25*P2) S3(I,J)=TANH(-0.25*P3+0.3) 10 CONTINUE CALL PFNAME('DEM10.PLO') CALL PLT00 CALL VISTR(-0.5,S1,S2,S3,0.5,NX,NX,NY,NY,0.2,0.2,0.2,-1,-1, * PLTCA) CALL PLTEND END