- 21? - 



27050 

 27100 

 27150 

 27200 

 27250 

 27300 

 27350 

 27400 

 27450 

 27500 

 27550 

 27600 

 27650 

 27700 

 27750 

 27800 

 27850 

 27900 

 27950 

 28000 

 28050 

 28100 

 28150 

 28200 

 28250 

 28300 

 28350 

 28400 

 28450 

 28500 

 28550 

 28600 

 28650 

 28700 

 28750 

 28800 

 28850 

 28900 

 28950 

 29000 

 29050 

 29100 

 29150 

 29200 

 29250 

 29300 

 29350 

 29400 

 29450 

 29500 

 29550 

 29600 

 29650 

 29700 

 29750 

 29800 

 29850 

 29900 

 29950 

 30000 



DO 445 m=1 ,nuf f 

 rnaxfrq(m) = 1 

 egymax = 0.0 

 DO 450 J=1 ,nn 

 if (A(j + 1 ). 1 t .EGYMAX) GO TO 450 

 EGYMAX = A(j+1) 

 maxfrq(m) = j+1 

 Cont Inue 

 k = maxfrq(m) 



EXCEED(m) = (A(k)-avrg) / avrg * 100.0 

 PD(m) = 6.28318/freq(k) 

 A(k) = 0.0 

 Continue 



WRITE(4,945) 



Format (//,' Maximum Energy located at bands: 

 WRITE (4, 946) (pd(m) ,exceed(m) ,m=1 .nuf f ) 

 Format ( f 15 . 3 , ' exceeds average energy by', 

 f8.2. ' %' ) 



Return 

 End 



subroutine FCOEFS (ab . abs , ab i , d1 , dpi , d1 1 , d2 , dp2 . 



1 d2 i . d3 . dp3 , d3 1 , d4 , dp4 , d4 i , dxx . dpxx , dxx 1 , dyy , 



2 dpyy ,dyy i , wvnr . f req, thtal , thta2, thta3, 



3 aO.al ,a2,a3.b1 ,b2.b3,blk. icut) 



C Finds the first nine directional fourier coefs 



C for a few selected frequencies. 



C Tl^is uses only one possible equation for each coef , 



C but calculates each coeficient with all possible 



C combinations of gages for that equation. 



DIMENSION ab(2050) .abs (4 100) .abi (2050) . 



1 d1(2050),dp1(4100).d1i(2050). 



2 d2(2050) ,dp2(4100) ,d21 (2050) 

 DIMENSION d3(2050) , dp3(4 100) , d3 i ( 2050) . 



1 d4 ( 2050 ). dp4( 4 100),d4 1(2050). 



2 dxx (2050) ,dpxx(4100) . dxx i( 2050), 



3 dyy(2050) .dpyy(4100) ,dyyi (2050) 

 DIMENSION wvnr(2050),freq(2050) 

 DIMENSION thtal (400) . thta2(400) , thta3(400) 

 DIMENSION a0(400),a1(400),a2(40O) .a3(400). 



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

 INTEGER blk 

 DIMENSION a1d1(400) , a 1d3(400) . a2d1d2(400) . 



1 a2d1d4(400) ,a2d3d2(400) .a2d3d4(400) .a3d1(400). 



2 a3d3(4OO),b1d2(40O),b1d4(4O0) , 



3 b2d2d1(400), b2d2d3(400),b2d4d1(400), b2d4d3(400), 



4 b3d2(400), b3d4(400) 



pi = 3. 14159 

 twopi = 6.283185 



DO 510 1=2, lout 



