- 213 - 



15050 

 15100 

 15 150 

 15200 

 15250 

 15300 

 15350 

 15400 

 15450 

 15500 

 15550 

 15600 

 15650 

 15700 

 15750 

 15800 

 15850 

 15900 

 15950 

 16000 

 16050 

 16100 

 16150 

 16200 

 16250 

 16300 

 16350 

 16400 

 16450 

 16500 

 16550 

 16600 

 16650 

 16700 

 16750 

 16800 

 16850 

 16900 

 16950 

 17000 

 17050 

 17100 

 17150 

 17200 

 17250 

 17300 

 17350 

 17400 

 17450 

 17500 

 17550 

 17600 

 17650 

 17700 

 17750 

 17800 

 17850 

 17900 

 17950 

 18000 



914 



WRITE(17,914) freq(i), ab(1 ) , dl ( i ) , cl2( i ) . d3( i ) , 

 d4( i ).dxx( i ),dyy( 1 ) 

 Cont 1 nue 

 Format(f 1 1.3,f11.5.4F14.9.2F20.18) 



C Identify bands of highest energy 



Cal 1 MAXSIG(ab, 1cut , maxf rq. n , dT , nuf f , f req) 



C Generate directional spectra for highest energy bands 



Call DIRXSP (AO.A1 , A2, A3.B1 ,B2,B3. f req. n, dT . nuf f , 

 1 maxf rq, orient ) 



C GO TO 9999 



WRITE(4.910) 



DO 666 nn=1 ,ri\jf f 

 LKM = maxfrq(nn) 

 PD(nn) = twopi/f req( 1 km) 



IF(thta1(lkm) . 1 t .0. ) thta1 ( 1 km) = 360. +thta1 ( 1 km) 

 IF(thta2( 1km) . 1 t .0. )thta2( 1 km) = 360 . +thta2( 1 km) 

 IF( thta3( lkm).lt.O.)thta3( 1 km) =360 . +thta3( 1km) 

 WRITE(4,913) PD(nn) . thtaK 1km) , thta2(lkm), 



1 thta3(lkm) , dir 1( 1 km) .di r2( 1 km) ,dir3( 1 km) .di r4( 1 km) , 



2 angi (1km) ,ang2( 1km) ,ang3( 1 km) , ang4( 1 km) 

 666 Continue 



910 Format(/.5x, ' MEAN DIRECTION OF PROROGATION ', 6x ,' ESTIM' , 



1 'ATE W/OUT ABSOLUTE GAGE :', 12x , 'WITH 1 CHANNEL'. 



2 ' AND ABSOLUTE'./.' Per iod' . 5x . ' n= 1 ' . 7x , ' n=2 ' . 



3 7x. 'n=3' ,Sy. ' 1 -2 ' , 6x . ' 3-4 ' , 6x . ' 2-3 ' . 6x . ' 1-4 ' . 



4 1 1 X . ' dp 1 ' . 6x . ' dp2 ' , 6x . ' dp3 ' . Gx , ' dp4 ' ) 

 913 Format(f6.2.3f 10.2.4f9. 1.5x.4f9. 1) 



SUBROUTINE INPUT(dpv1 . dpv2 ,dpv3 , dpv4 . absv) 



DIMENSION dpv1(4100).dpv2(4100).dpv3(4100), 

 1 dpv4(4100) .absv(4100).M(4100) 



INTEGER dev.yr.day.hr. A( 10) 



WRITE(4.890) 



READ( 18. 900) dev.yr . mo. day , hr.min 



WRITE(4,901 ) dev.yr . mo . day , hr . min 



n=1 



DO 200 J = 1 .410 



READ(18,905) ( A( 1 ) . i = 1 . 10) . LAST 

 DO 100 1=1 . 10 

 DPVI(n) = A(i)/265.0 * 1.65 

 100 n=n+1 

 200 Continue 



