o o o o
469
IF(Z.GT.-ACCUR(1).AND.Z.LT.(H(1)+ACCUR(1))) GO TO 230 BI
200 WRITE(NOUT,9110) BI
GO TO 400 BI
210 IF(LAYER.LT.NLAYS) GO TO 220 BI
IF(Z. GT (H(M) -ACCUR( 1) ) ) GO TO 230 BI
GO TO 200 BI
220 IFCZ.GT.(H(LAYER-1)-ACCUR(1)).AND.Z.LT.(H (LAYER)+ACCUR(1)BI
& )) GO TO 230 BI
GO TO 200 BI
230 RADI=RADIUS(I) BI
LOAD=LDSTRS(I) BI
HOSTRS=HOSTR(I) BI
RLOW=R.LT.ACCUR(1) BI
ST=SIN(THETA) BI
CT=COS(THETA) BI
c BI
C CONPNT DETERMINES FOR EACH POINT-LOAD BI
C CONFIGURATION WHICH INTEGRALS HAVE TO BE BI
C CALCULATED. BI
c BI
CALL CONPNT(R,HOSTRS,LOAD,Z,N2,L2) BI
IF(LAYER.NE.1) GO TO 250 BI
CZ = DBLE(Z) BI
IF(Z.LT.ACCUR(1).AND.ABS(R-1.0).LT.ACCUR(1)) GO TO 240 BI
c BI
C ASYMPT DETERMINES THE LIPSCHITZ-HANKEL BI
C INTEGRALS NEEDED FOR THE ASYMPTOTIC PART BI
C OF THE INTEGRALS,FOR POINTS IN THE TOP- BI
C LAYER ONLY. BI
c BI
CALL ASYMPT(R,ACCUR(1)) BI
GO TO 250 BI
C BI
C FOR POINTS AT THE RIM OF THE LOAD THE BI
C LIPSCHITZ-MANKEL INTEGRALS CAN BE GIVEN BI
C DIRECTLY. BI
C BI
240 F10M1 = 0.63662 BI
F100 = 0.5 BI
F11M1 =0.5 BI
F11M2 = 0.424413 BI
FIO 1 = 0.0 BI
F110 = 0.0 BI
Fill = 0.0 BI
BI
COMPUTATION OF THE REQUIRED INTEGRALS BY BI
CALLING IN GENDAT AND INGRAL BI
BI
250 INTT =0 BI
DO 260 J = 1,17 BI
INT(J) = 0.0 BI
260 CONTINUE BI
DO 270 J = 1,10 BI
INTV(J) =0 BI
K = IVER2CJ) BI
IF(.NOT.EPS(K)) GO TO 270 BI