223 - 



45050 

 45100 

 45150 

 45200 

 45250 

 45300 

 45350 

 45400 

 45450 

 45500 

 45550 

 45G00 

 45650 

 45700 

 45750 

 45800 

 45850 

 45900 

 45950 

 46000 

 46050 

 46100 

 46150 

 46200 

 46250 

 46300 

 46350 

 46400 

 46450 

 46500 

 46550 

 46600 

 46650 

 46700 

 46750 

 46800 

 46850 

 46900 

 46950 

 47000 

 47050 

 47100 

 47150 

 47200 

 47250 

 47300 

 47350 

 47400 

 47450 

 47500 

 47550 

 47600 

 47650 

 47700 

 47750 

 47800 

 47850 

 47900 

 47950 

 48000 



DO 601 NN=1 ,nuff 

 save = 0.0 

 I = maxfrq(nn) 

 PD(nn) = 6.2831/freq( i ) 



WRITE (4, 439) pd( nn ) .aO( 1 ) , a 1 ( 1 ) . a2( 1 ) , a3( 1 ) , 

 1 b1( 1 ),b2( 1 ) .b3( 1 ) 

 439 Format (f 6. 2, ' sees ' , 7f 16 . 4 ) 



DO 604 m=1 , 181 

 deg(m) = save 

 as = aO( 1 ) + a1(1)*w1 * cosd(deg(iii) ) 



1 + a2(1)*w2 * cosd(2.0*deg(n!) ) 



2 + a3(1)*w3 * cosd(3.0*deg(m) ) 

 bs = b1(1)*w1 * s1nd(deg(ni)) 



1 + b2(1)*w2 * s1nd(2.0*deg(ni)) 



2 + b3(1)*w3 * s1nd(3.0*deg(m) ) 



604 

 601 



S ( nn , m ) 



save = 

 Cont Inue 

 Cont Inue 



save + deginc 



C Determine the direction of greatest energy for 



C each band and convert to true north heading from 



C which waves propogate. 



WRITE(4,G002) 

 DO 613 NN=1 ,NUFF 



mBIG = O 



BIG = 0.0 



DO G12 m=1, 181 



IF(S(nn,m) .gt .BIG)mBIG = m 



612 IF(S(NN,M).gt.BIG)BIG=S(NN,M) 

 BIGDEG = mBIG * deginc - deginc 

 BIGDEG = ORIENT - BIGDEG 



IF (BIGDEG. 1 t O. )BIGDEG=BIGDEG+360.0 

 WRITE(4,976) nn, BIGDEG 



613 Continue 



976 Format(' Band', 12,': Greatest Energy at',f7.2, 

 1 ' deg' ) 



C Record the directional spectra for each band. 



C Directions correspond to angle from which waves 



C propagate with respect to true north. 



WRITE(3.969) (pd(1), 1=1.nuff) 



DO 615 m=1 . 181 



DEG(m) = ORIENT - DEG(m) 

 IF(DEG(m) . It .0. )DEG(m)=DEG(m)+360 .0 

 WRITE (3, 970) deg(m) . (S(nn,m) . nn=1 .nuff ) 

 615 Continue 



969 Format (//.' DIRECTIONAL - FREQUENCY SPECTRUM'.//, 

 1 ' dirxn' ,9x,6f 17.3, 'sees' ,/) 



970 Format ( f 8 . 1 . ' deg' , 3x,6f 17 . 3) 

 6002 Format (/) 



