399 TEMP = (RH0(I-1) / RHO(I)) / TEMDI 



400 TtMl =US (1-1) * G(I-1) + ZBI(I-I) * GI(I-lj) / TEMDR 



401 * -(ZT (I) * G(I) + ZTI(I) * GI(I)) * TEMP 



402 TEM1I =(ZBI(I-1) * G(I-1) - 2B (1-1) * GI(I-I)) / TEMDR 



403 * -(ZTI(I) * G(I) - ZT (I) * GI(I)) ♦ TEMP 



404 TEMRSP = AF( 1-1 ,NN)*B( 2*1-1 ,2) - AG( 1-1 ,NN)*BI( 2»I-1 ,2) ■( 



405 1 BF{ 1-1 ,NN)*B( 2*1-1 ,1) - BG( 1-1 ,NN)*BI( 2*1-1 ,1) 



406 TEMISP = AG( 1-1 ,NN)*B( 2*1-1 ,2) + AF( 1-1 ,NN)*BI( 2*1-1 ,2) ^ 



407 1 BG( 1-1 ,NN)*B{ 2*1-1 ,1) + BF( 1-1 ,NN)*BI( 2*1-1 ,1) 



408 AX2 = TEMRSP*TEMRSP - TEMISP*TEMISP 



409 AX2I = TEMRSP * TEMISP 



410 AX2I = AX2I + AX2I 



411 TEMDR = RHO(I-I) / (G CU(I-1)**2 + G CUI(I-1)**2) 

 41? TEMDI = RHO(I ) / (G CU( I )**2 + G CUI(I )**2) 



413 TEM2 = G CU( 1-1 ) * TEMDR - G CU(I) * TEMDI 



414 TEM2I = G CUI(I) * TEMDI - G CUI(I-I) * TEMDR 



415 TEMR1 = AX1+TEM1 - AX1I*TEM1I 



416 TEMI1 = AX1I*TEM1 + AX1*TEM1I 



417 TEMR2 = AX2 * TEM2 - AX2I * TEM2I 



418 TEMI2 = AX2I * TEM2 + AX2 * TEM2I 



419 D(NN) = D(NN) + TEMR1 / RH0(I-1 ) + TEMR2 



420 rH(NM) -r nT(NN) + TEMI1 / RHO(I-I) + TEMI2 



421 111 CONTINUE 



422 IF ( K1 .GT. 3) DA(NN) = DSQRT ( (D(NN) ♦*2 + DI(NN)**2) ♦ FREQ / 



423 * PHASE V(NN) ) 



424 EIGEN(NN) = LAMBDA 



425 EIGENI(NN) = LAMBDI 



426 IF (KB .GT. 2) GO TO 131 



427 L = 



428 K = 24 



429 DO 112 I = 1 ,N 



430 L = L + 1 



431 COL(L) = SNGL(ZT(I)) * 100. 



432 L = L + 1 



433 COL(L) = SNGL(ZTI(I)) * 1000. 



434 K = K + 1 



435 COL(K) = SNGL(ZB(I)) * 100. 



436 K = K + 1 



437 COL(K) = SNGL(ZBI(I)) * 1000. 



438 112 CONTINUE 



439 PRINT 130, (COL(I), 1=1, L) 



440 PRINT 130, (COL(I), 1=25, K) 



441 130 FORMAT (4H Z = , 11(16.15)) 



442 M = N + N 



443 PRINT 132, (RATIO(I), I = 1,M) 



444 132 FORMAT ( 1 1 ( 1 X , 2F5 . 3 ) ) 



445 131 DB LOSS(NN) = - LAMBDI * 8686. DO 



446 PHINV = V * PHASE V(NN-1) /((V - PHASE V(NN-1))* FREQ) 



447 PRINT 109, NN , E IGEN( NN) , E IGENI ( NN ), D( NN) , DI ( NN ), PHINV , 08 LOSS(NN) 



448 109 FORMAT ( 3H N=,I5,10H LAMBDA =,2E15.7,4H D= ,2E15.7, 



449 * 12H INT RANGE = ,F8.0, 6H L/K =, F8.5) 



450 IF (K3 .EO. 0) GO TO 50 



451 CALL RCOEF (K3) 



452 50 CONTINUE 



453 C READ IN SOURCE AND RECEIVERS DEPTHS 



454 301 NRT = NR 



455 NR = 



51 



