456 NN = NN - 1 



457 tMPI = Kl + 1 



458 IF (K1 .NE. 3) GO TO 321 



459 NR = NRT 



460 GO TO 501 



461 321 READ 20, SOURCE 



462 320 NR = NR + 1 



463 READ 20, RECVRS(NR), FINAL, STEPP 



464 IF (NR.GT.50) GO TO 300 



465 350 IF (RECVRS(NR) . EO • . ) GO TO 300 



466 310 IF (FINAL .EO-O.) GO TO 320 



467 330 RECVRS(NR+1) = RECVRS(NR) + STEPP 



468 IF (RECVRS(NR + 1 ) . GT . FINAL) GO TO 320 



469 340 NR = NR + 1 



470 IF(NR .GT. 50) GO TO 300 



471 GO TO 330 



472 300 PRINT 303 



473 303 FORMAT (/21H SOURCE AND RECEIVERS ) 



474 PRINT 21 , (DEPTH( I ) , I = I.NR) 



475 21 FORMAT (8F10.2) 



476 C COMPUTE DEPTH FUNCTIONS 



477 DO 500 I = 1 ,NN 



478 LOC = 1 



479 DO 305 J = 1 ,NR 



480 IF ((J .EO. 1) .AND. (Kl .GT. -5)) GO TO 305 



481 LCTR = 



482 380 IF((DEPTH(J) . GE . Z ( LOC ) ) . AND. ( DEPTH( J ) .LT. Z( LOC+1 ) ) ) GO TO 360 



483 371 IF (LOC .GE. N) GO TO 385 



484 370 LOC = LOC + 1 



485 GO TO 380 



486 385 IF (DEPTH(J) .GE. Z ( LOC ) ) GO TO 360 



487 390 L0C=1 



488 LCTR=LCTR+1 



489 IF (LCTR .GT. 2) GO TO 305 



490 GO TO 380 



491 360 XI = CAY (LOC) - EIGEN (I) 



492 X2 = CAY (LOC) + EIGEN (I) 



493 X3 = CAYI(LOC) - EIGENI(I) 



494 X4 = CAYI(LOC) + EIGENI(I) 



495 TEMP = X1 * X2 - X3 * X4 

 4P6 TEWPI - XI + X4 + X3 * X2 



497 TEMDEN = G SO(LOC) **2 + G S0I(L0C)**2 



498 ZE = (TEMP * GSO(LOC) + TEMPI ♦ G SQKLOC)) / TEMDEN 



499 ZEI = (TEMPI * GSO(LOC) + TEMP ♦ GSOKLOC)) / TEMDEN 



500 TEM1 = ZE 



501 IF (ZE .GT. -7.5) GO TO 438 



502 S = CAY( LOC) 



503 T = CAY I (LOC) 



504 DO 437 K = 1 ,20 



505 TEMP = 5**2 + T**2 



506 TEMPI = (EIGENI(I) * S - EIGEN(I) • T) / TEMP 



507 TEMP = (EIGEN(I) * S + EIGENI(I) * T) / TEMP 



508 ZE = ((1. + TEMP) * (1. - TEMP) + TEMPI**2) * CON(LOC) 



509 ZEI = -2. * TEMPI * TEMP * CON(LOC) 



510 ZR = ZE / -7.5 



511 IF (DABS(ZR-1.) .LT. 1.D-3) GO TO 438 



512 S = EIGEN(I) + (S - EIGEN(I)) / ZR 



52 



