C [DEMO3] C FLAT-BOTTOM CRATER ON HILL C [03-JUN-74] PROGRAM COPIED FROM DECUS DOCUMENTATION C 13-OCT-75 EDER EXTERNAL PLTCA DIMENSION WY(71,71) NX=71 NY=71 DO 10 I=1,NX DO 10 J=1,NY X1=0.1*FLOAT(I-40) Y1=0.1*FLOAT(J-60) RR1=X1*X1+2.*Y1*Y1+2.*X1*Y1 R1=SQRT(RR1) P1=RR1*RR1*EXP(-R1) ZW=0.2*P1+0.005*X1*X1*Y1 10 WY(I,J)=TANH(ZW) CALL PLT00('TST3') CALL PLTSE (-1.0,WY,1.0,NX,NX,NY) C CALL PLTSW (-1.0,WY,1.0,NX,NX,NY) CALL EXIT END