- 222 - 



42050 



42 100 

 42150 

 42200 

 42250 

 4 2300 

 42350 

 42400 

 42450 

 4 2500 

 42550 

 42600 

 42650 

 42700 

 42750 

 42800 

 42850 

 42900 

 42950 

 4 3000 

 43050 



43 100 

 43150 

 43200 

 43250 

 43300 

 43350 

 43400 

 43450 

 43500 

 43550 

 43600 

 43650 

 43700 

 43750 

 43800 

 43850 

 43900 

 43950 

 44000 

 44050 



44 100 

 44150 

 44200 

 44250 

 44300 

 44350 

 44400 

 44450 

 44500 

 44550 

 44600 

 44650 

 44700 

 44750 

 44800 

 44850 

 44900 

 44950 

 45000 



DO 58 1=2, icut 



PK = prsp( 1 ) * gamma 



abs 



1) 



= abs( 



) / 



PK 



abi 



1) 



= ab 



(1) / 



PK 



a(1 



= 



a(l) 



/ 



PK 





b(1 



= 



b(1) 



/ 



PK 





c(1 



= 



c(1) 



/ 



PK 





dd 



= 



d(1) 



/ 



PK 





e(1 



= 



e(1) 



/ 



PK 





f(1 



= 



f(1) 



/ 



PK 





g(i 



= 



g(i) 



/ 



PK 





h(1 



= 



h(1) 



/ 



PK 





P(1 



= 



P{1) 



/ 



PK 





q(1 



= 



q(1) 



/ 



PK 





r(i 



= 



r(1) 



/ 



PK 





s(1 



= 



s(1) 



/ 



PK 





Cont Ini 



je 











Return 













End 













Subroutine DIRXSP ( aO. a 1 , a2 , a3 . 

 1 b1, b2, b3, freq,n,dT,nuff, ma xfrq, orient) 



We weight the fourier series directional coefs with 

 cosine**2s factors. Here, s=2. Then we create a 



matrix of the frequencies of Interest analyzed at 



incremental degrees. 



DIMENSION aO(400) ,a1 (400) ,a2(400) ,a3(400) .b1 (400) . 



1 b2(400),b3(400) ,freq(2050) ,maxfrq(20), 



2 S(20, 186),deg( 186) ,pd(20) 



DEGINC is the increment in degrees for which the 

 directional spectra is calculated over. 

 DEGINC = 2.0 



Weighting factors to smooth directional spectra and 

 eliminate negative side energy lobes: 

 W1 = 66. /77. 

 W2 = 15. /28. 

 W3 = 5./2 1 . 



Weighting factors for non-smoothed directional spectra: 

 W1 = 1 . 

 W2 = 1 . 

 W3 = O. 



Print out the directional Fourier coefficients, 

 if desired. The header follows: 



WRITE(4,438) 

 38 Formate/,' PERIOD '. I6x ,' AO' , 14x ,' A 1 ', 14x ,' A2 ', 14x . 



1 'A3' , 14x, 'B1 ' , 14x, 'B2' . 14x, 'B3' ./, ' ',16x, 



2 ' ',14x,' ',14x,'__'.14x,'__',14x,' ',14x, 



3 ' ',14x,' ',/) 



D5 Continue 



