471













 INTV(J) = K BI
 INTT = INTT+1 BI
 270 CONTINUE BI
 IF(INTT.EQ.0) GO TO 280 BI
 IF(NLAYS.NE.1) CALL GENDAT(1,NZEROS,R,ACC) BI
 CALL INGRAL(2,INTV,INTT,INT) BI
 280 INTT = 0 BI
 DO 290 J = 1,7 BI
 INTV(J) = 0 BI
 K = IVER1(J) BI
 IF(.NOT.EPS(K)) GO TO 290 BI
 INTV(J) = K BI
 INTT = INTT+1 BI
 290 CONTINUE BI
 IF(INTT.EQ.0) GO TO 300 BI
 IF(NLAYS.NE.1) CALL GENDAT(0,NZEROS,R,ACC) BI
 CALL INGRAL(1,INTV,INTT,INT) BI
 300 PSIO = PSI(I) BI
C-----------------------------------------------------------------------BI
C CALC COMPUTES AND OUTPUTS THE STRESSES, BI
C STRAINS AND DISPLACEMENTS,INDUCED BY EACH BI
C LOAD SEPARATELY. BI
C------------------------------------------------------------------------ B
 CALL CALC(INT,V,R,MU,RADI,FT,LOAD,HOSTRS,PSIO,Z) BI
 IF(.NOT.N2) GO TO 330 BI
C------------------------------------------------------------------------BI
C COMPUTATION AND SUMMATION OF CARTESIAN BI
C COORDINATES.THE USED COORDINATE SYSTEM IS BI
C THE ONE WHEREIN POINTCOORDINATES WERE BI
C STATED. BI
C-----------------------------------------------------------------------BI
 UZ =UZ+V(3) BI
 IF(ABS(RADDIS).LT.ACCUR(1)) GO TO 310 BI
 CT = (AX(IPOS)-X(I))/RADDIS BI
 ST = (AY(IPOS)-Y(I))/RADDIS BI
 GO TO 320 BI
 310 CT =COS(ETA) BI
 ST =SIN(ETA) BI
 320 CT2 =CT*CT BI
 ST2 =ST*ST BI
 STCT =ST*CT BI
 A(1,1)=A(1,1)+V(4)*CT2+V(5)*ST2-2.0*V(7)*STCT BI
 A(1,2)=A(1,2)+V(7)*(CT2-ST2)+(V(4)-V(5))*STCT BI
 A(1,3)=A(1,3)+V(8)*CT-V(9)*ST BI
 A(2,1)=A(1,2) BI
 A(2,2)=A(2,2)+V(4)*ST2+V(5)*CT2+2.0*V(7)*STCT BI
 A(2,3)=A(2,3)+V(8)*ST+V(9)*CT BI
 A(3,1)=A(1,3) BI
 A(3,2)=A(2,3) BI
 A(3,3)=A(3,3)+V(6) BI
 UX =UX+V(1)*CT-V(2)*ST BI
 UY =UY+V(1)*ST+V(2)*CT BI
 330 CONTINUE BI
 TRACE=A(1,1)+A(2,2)+A(3,3) BI
 AB =(1.0+MU)/E(LAYER) BI
 AC =MU*TRACE/E(LAYER) BI