221 



39050 

 39100 

 39150 

 39200 

 39250 

 39300 

 39350 

 39400 

 39450 

 39500 

 39550 

 39600 

 39650 

 39700 

 39750 

 39800 

 39850 

 39900 

 39950 

 40000 

 40050 

 40100 

 40150 

 40200 

 40250 

 40300 

 40350 

 40400 

 40450 

 40500 

 40550 

 40600 

 40650 

 40700 

 40750 

 40800 

 40850 

 40900 

 40950 

 41000 

 4 1050 

 41 100 

 4 1 150 

 4 1200 

 41250 

 4 1300 

 41350 

 4 1400 

 41450 

 4 1500 

 4 1550 

 4 1600 

 41650 

 4 1700 

 4 1750 

 4 1800 

 41850 

 4 1900 

 41950 

 4 2000 



IF( ICO. EO. 1 ) GO TO 11 

 TR-WR*FR( J)-WI*FI(d) 

 TI=WR*FI(d)+WI*FR(d) 

 GO TO 12 



11 TR=WR*FR(d)+WI*FI(d) 

 TI=WR*FI( J)-WI*FR( J) 



12 FR(d)=FR(I)-TR 

 FI(J) = FI(I )-TI 

 FR(I ) = FR(I)+TR 



4 FI(I )=FI( I) + TI 

 L-ISTEP 



GO TO 3 

 7 Continue 

 AN = N 



IFdCO.EO. 1 ) GO TO G 

 DO 5 1 = 1 ,N 

 FR(I ) = FR(I )/AN 



5 FI(I ) = -FI(I )/AN 



6 RETURN 

 END 



SUBROUTINE WVLEN(DPT , SIG. XKH) 



THIS SUBROUTINE CALCULATES LINEAR WAVELENGTH BY NEWTONS 

 METHOD 



PI=3. 14159265 

 TW0PI=2.0*PI 

 XKH0=SIG**2*DPT/32.2 

 IF(XKH0-6. 3) 2.1.1 

 XKH^XKHO 

 GO TO 9 

 XKH=SORT(XKHO) 

 SH=SINH(XKH) 

 CH=COSH(XKH) 

 EPS=XKHO-XKH*SH/CH 

 SLOPE =-XKH/CH**2-SH/CH 

 DXKH=-EPS/SLOPE 

 I F ( ABS ( DXKH/XKH ) -O . 000 1 ) 

 XKH=XKH+DXKH 

 GO TO 3 



XLENTH=TWOPI*DPT/XKH 

 RETURN 

 end 



subroutine HYDRO ( i cut , prsp, wvnr , abs , ab i , a ,b. c , d, 

 1 f . g, h.p.q, r . s .gamma ) 

 C Divides through by pressure response function 



C and gamma . 



C Operates on the a(n)' and b(n)' terms from 



C each of the five gages' FFT results. 



DIMENSION prsp(2050) .w.'nr(2050) 



DIMENSION abs(2050) . abi(2050). a( 4 100) . b( 2050) . 



1 c(4100).d(2050).e(4100),f(2050).g(4100),h(2050), 



2 p(4100).q(2050),r(4 100),s(2050) 



