211 - 



09050 

 09100 

 09150 

 09200 

 09250 

 09300 

 09350 

 09400 

 09450 

 09500 

 09550 

 09600 

 09650 

 09700 

 09750 

 09800 

 09850 

 09900 

 09950 

 10000 

 10050 

 10100 

 10150 

 10200 

 10250 

 10300 

 10350 

 10400 

 10450 

 10500 

 10550 

 10600 

 10650 

 10700 

 10750 

 10800 

 10850 

 10900 

 10950 

 1 1000 

 1 1050 

 1 1 100 

 1 1 150 

 1 1200 

 1 1250 

 1 1300 

 1 1350 

 1 1400 

 1 1450 

 1 1500 

 1 1550 

 1 1600 

 1 1650 

 1 1700 

 1 1750 

 11800 

 11850 

 1 1900 

 1 1950 

 12000 



g3 = gamma * dxinS 

 g4 = gamma * dx1n4 



TILT1 = atan2(t1.g1) * 360./twop1 



TILT2 = atan2(t2.g2) * 360./twopi 



TILT3 = atan2(t3.g3) * 360./twopi 



TILT4 = atan2(t4.g4) ♦ 360./twopi 



WRITE(4.119) t11t1,t11t2,tnt3.tnt4 

 119 Format(/,' POSSIBLE INSTRUMENT TILT (deg)',/. 



1 5x. 'DPI ' .6x, 'DP2' ,6x, 'DP3' .Sx, 'DP4' ,/.4f9.3,/) 



; Create acceleration terms 

 DO 56 1=1 .N 



DPXX(i) = (DP3(i) - DP1(I))/arm13 

 DPYY(i) = (DP4(i) - DP2( 1 ) )/arm24 

 56 Continue 



DO 7034 1=1,600 



WRITE( 14,7033) dpxx( 1 ) , dpyy ( i ) 

 7034 Continue 

 7033 Format(2f 16.9) 



CALL FFT(dp1 ,d1i .k,0) 



CALL FFT(dp2,d2i ,k.O) 



CALL FFT(dp3.d3i .k,0) 



CALL FFT(dp4,d4i ,k,0) 



CALL FFT(dpxx,dxxi ,k,0) 



CALL FFT(dpyy,dyyi ,k,0) 



CALL FFT(abs,abi ,k,0) 



Generate Pressure Response Function and Wavenumber 



using the calculated water depth. Sensors are average 

 height of 3.79 feet above bottom. 



bo\ tcm 

 DO 762 



+ 3.79 



(i-1) 



depth 

 =2 , icut 

 Sigma = twop1/(N*dT) 

 f req( 1 ) = Sigma 

 call WVLEN (bottom, Sigma, www) 

 wvnr(i) = www/bottom 

 PRSP( 1 ) -cosh(wvnr( i )* (bottom- sgage) ) 

 /cosh(wvnr( i )*bottom) 



Divide Fourier coefficients by pressure response 

 function and seawater specific gravity. 



CALL HYDRO ( 1 cu t , prsp , wvnr , abs . ab i ,dp1 ,d1 1 , dp2 , d2 i 

 dp3 , d3 i , dp4 . d4 i , dpxx , dxx i , dpyy , dyy i , gamma ) 



Create power spectra 

 DO 67 1=2, icut 



AB( 1 ) = 2.0*(abs( i 



Dl(i) = 

 D2(1) = 

 D3(i) = 

 D4(1) = 

 DXX( i ) 



2.0* (dpi ( i 

 2.0*(dp2( i 

 2.0*(dp3( i 

 2.0*(dp4( i 

 = 2.0*(dpxx 



**2 



**2 

 **2 

 **2 

 **2 

 1)*' 



ab1( i )**2) 

 d11( i )**2) 

 d2i( i )**2) 

 d31 ( i )**2) 

 d4i( i )**2) 

 + dxxi ( 1 )**2) 



