513 



437 



514 



438 



515 





516 





517 





518 





519 





520 



439 



521 





522 





523 





524 





525 



442 



526 



302 



527 





528 





529 





530 





531 



305 



532 





5?T 





534 



451 



535 



270 



536 





537 



432 



538 





539 





540 





541 



431 



542 





543 





544 





545 



433 



546 





547 





548 





549 





550 





551 



435 



552 





553 



436 



554 





555 



452 



556 



C PI 



557 



420 



558 



210 



559 





560 



220 



561 





562 





563 





564 



240 



565 





566 





567 



250 



568 



225 



569 



260 



CONTINUE 



IF (G(LOC) .LT. 0.) GO TO 439 



ZE = G(LaC) * (DEPTH(d) - Z(LOC)) + TEM1 



IF (ZE .GT. -7.5) GO TO 442 



F(J) = 1.D-12 



FI(J) = 0. 



GO TO 305 



ZE = G(LOC) * (DEPTH(d) - Z ( LOC ) ) + ZE 



IF (ZE .GT. SLIM) GO TO 442 



F(iJ) = 1 .D-12 



FI(d) = 0. 



GO TO 305 



ZEI =GI(LOC) * (DEPTH(lI) - Z(LOC)) + ZEI 



CALL HANKEL(ZE,ZEI , 1 ) 



F(J) =(AF(LOC, I)*H1R - AG( LOC , I ) *H1 I + BF( LOC , I )*H2R - BG(LOC,I) 

 1 *H2I) * RHO(LOC) 



FI(J) =(AG(LOC, I)+H1R + AF ( LOC , I ) *H1 I + BG(L0C,I)*H2R + BF(LOC,I) 

 1 *H2I ) * RHO( LOC) 



CONTINUE 



IF (K1 .EQ. 2) GO TO 451 



GO TO 432 



PRINT 270, DEPTH(NR) 



FORMAT (7H1 DEPTH ,F5.1 ,6X,3HE-8, 1 7X , 3HE-6 , 1 7X , 3HE-4 , 1 7X , 3HE-2, 

 + 17X,3HE ) 



IF (K1 . LT. 4) GO TO 431 



IF (K1 .GT. 5) GO TO 433 



SRES(I) = (F(1)**2 + FI(1)**2) / DA(I) 



GO TO 500 



TEMDEN = D(I)+D(I) + DI(I)*DI(I) 



TEMRE = F(1)*D(I) + FI(1)*DI(I) 



FD = TEMRE/TEMDEN 



FDI = (D(I) ♦ FI(1) - DI(I) * F(1)) / TEMDEN 



DO 400 K = 2,NR 



d = K - 1 



L=d*NN-NN+I 



IF (K1 . LT. 6) GO TO 435 



FF = SRES(I) ♦ (F(K)**2 + FI(K)**2) / DA(I) 



GO TO 436 



FF = FD * F (K) - FDI * FI(K) 



FFI = FD * FI(K) + FDI * F(K) 



UU(L) = FF 



UUI(L) = FFI 



GO TO (400,410,420,400,400,400,400,400), K1P1 

 PLOT DEPTH FUNCTIONS 



DO 210 II = 1 , 120 



COL(II)= 1H 



DO 220 11= 20, 100, 20 



COL(II)= 1HI 



FE = FF * FF + FFI * FFI 



IF ( (FE.GT. lE-20) .AND. (FE.LT. 10000. )) GO TO 240 



GO TO 250 



INT = 100. DO + 2.17147D0 * DLOG(FE) 



COL(INT) = 1H* 



GO TO 225 



C0L(2) = 1H* 



PRINT 260, COL 



FORMAT (120A1) 



53 



