lie - 



24050 

 24 100 

 24150 

 24200 

 24250 

 24300 

 24350 

 24 400 

 24450 

 24500 

 24550 

 24600 

 24650 

 24700 

 24750 

 24800 

 24850 

 24900 

 24950 

 25000 

 25050 

 25100 

 25150 

 25200 

 25250 

 25300 

 25350 

 25400 

 25450 

 25500 

 25550 

 25600 

 25650 

 25700 

 25750 

 25800 

 25850 

 25900 

 25950 

 26000 

 26050 

 26100 

 26150 

 26200 

 26250 

 26300 

 26350 

 26400 

 26450 

 26500 

 26550 

 26600 

 26650 

 26700 

 26750 

 26800 

 26850 

 26900 

 26950 

 27000 



If (dp4 

 if (dp4 

 tanu = 

 tanv = 

 tanw = 

 tanz = 

 D12 = 

 D34 = 

 D23 = 

 D14 = 

 snum = 

 denom 



DIRKj ) 

 znum = 

 denom 



DIR2(j) 

 znutn = 

 denom 



DIR3(J ) 

 znum = 

 denom 



DIR4( j ) 



( J ) . ne . O 



( j ) . eq . O 



tan2 - 



tan4 - 



tan2 - 



tan4 - 



sqrt(d1( 



sqrt(d3( 



sqrt(d3( 



sqrt(d1( 



2*D12 

 = d1(J) 

 =ORIENT- 



2*D34 

 = d3(j) 

 =ORIENT- 



2*D23 

 = d3(j) 

 =ORIENT- 



2*D14 

 = d1(J) 

 =ORIENT- 



tan4=atan2(d4i (J ) .dp4( j )) 

 jn4 = 0.0 



d2(j) 

 ( 180./3 . 14 159)*0.5*atan2(znum.denom) 



d4(J) 

 ( 180. /3. 14159)*0.5*atan2(znum, denom) 



d2(j) 

 ( 180./3. 14159)*0.5*atan2(2num,denorn) 



d4(j) 

 ( 180. /3 . 14159)*0. 5 *atan2( znum, denom) 



IF(d1r1( j ) . 1 t 

 IF(d1r2(J ) . It 

 IF(d1r3( j ). It 

 IF(d1r4( j ) . It 



Cent inue 



0. )d1r1( j )=dir1(J )+360.0 

 0. )dir2( j )=dir2( j ) +360.0 

 0. )dir3( j )=d1r3( j ) +360.0 

 0. )dir4( j )=dir4( j ) + 360.0 



CALL BLOCK( 1cut ,dir1 .maxblk.blk) 

 CALL BLOCK( 1cut.dir2. maxblk.blk) 

 CALL BLOCK( 1cut,dir3, maxblk.blk) 

 CALL BLOCK( 1cut,d1r4. maxblk.blk) 



Return 

 End 



SUBROUTINE MAXSIG ( R , icut , maxf rq, n. dT , nuf f . f req) 

 C Locates "nuff" number of peak energy bands 

 C and reports their location and relative magnitude. 



Dimension R(2050), maxf rq( 20) . f req(2050) . PD(20) , 

 1 exceed(20).A(2050) 



NN = Icut - 1 

 TOT = 0.0 



C Identify average energy in spectra 

 DO 440 1=1, Icut 

 A(l) = R(1) 

 440 TOT = A( 1 ) + tot 

 AVRG = TOT/nn 



C Identify "nuff" number of peak bands. Calculate each 



C peak band's magnitude beyond the average. Set band = 



C after identification as a maximum in the temporary array. 



