SUBROUTINE S11 (NUMB,NCD,M) DIMENSION F0(2),ND(2) DIMENSION NDP (4,5,15), SPOT (4) COMMON/AR/BRAY (4,600),STRESS (200,4) COMMON/SIHI/IVECT(1800) DATA F0(1),F0(2)/3HNDP,4H SRC/ DO 30 I=1,NUMB SPOT (1)=STRESS (I,1) SPOT (2)=STRESS (I,2) SPOT (3)=SPOT (1) SPOT (4)=SPOT (2) SS=STRESS(I,3) C SCALED FROM ZMIN TO ZMAX IN RANGE 1 TO 10. C L=IFIX(SS) ND(1)=L/10 ND(2)=L-ND(1)*10 DO 30 JK=1,2 L=ND(JK) IF (L.EQ.0) L=10 IF (L.LT.0) L=11 IF (L.GT.15) L=10 E=FLOAT (M) EF=FLOAT(-JK) EF=E*EF DO 30 J=1,5 NCP=0 DO 29 N=1,4 IF (NDP(N,J,L).NE.0)NCP=2 29 CONTINUE IF (NCP.EQ.0) GO TO 31 K=K+1 DO 30 N=1,4 F=FLOAT (NDP(N,J,L)) BRAY (N,K)=SPOT (N)+E*F+EF 30 CONTINUE 31 CONTINUE NCD=K-1 IR=1 ITAG=11 CALL VECTOR(1,BRAY,1,NCD,IVECT,IR,LA,7,0,.FALSE.) CALL DISPLY(6,ITAG,IVECT,IR,LA) C NCD RETURNS THE SIZE OF THE 4XNCD BRAY ARRAY RETURN END