36050 



36 1 GO 

 36150 

 36200 

 36250 

 3G300 

 3G350 

 36400 

 36450 

 36500 

 36550 

 36600 

 36650 

 36700 

 36750 

 36800 

 36850 

 36900 

 36950 

 37000 

 37050 



37 100 

 37150 

 37200 

 37250 

 37300 

 37350 

 37400 

 37450 

 37500 

 37550 

 37600 

 37650 

 37700 

 37750 

 37800 

 37850 

 37900 

 37950 

 38000 

 38050 



38 100 

 38150 

 38200 

 38250 

 38300 

 38350 

 38400 

 38450 

 38500 

 38550 

 38600 

 38650 

 38700 

 38750 

 38800 

 38850 

 38900 

 38950 

 39000 



CALL 



BLOCK 



Icut 



CALL 



BLOCK 



1cut 



CALL 



BLOCK 



Icut 



CALL 



BLOCK 



Icut 



CALL 



BLOCK 



icut 



CALL 



BLOCK 



icut 



CALL 



BLOCK 



icut 



CALL 



BLOCK 



(cut 



CALL 



BLOCK 



icut 



CALL 



BLOCK 



icut 



CALL 



BLOCK 



icut 



CALL 



BLOCK 



icut 



CALL 



BLOCK 



icut 



CALL 



BLOCK 



icut 



CALL 



BLOCK 



icut 



CALL 



BLOCK 



icut 



CALL 



BLOCK 



icut 



CALL 



BLOCK 



icut 



ICUT 



= MAXt 



)1l< 



Return 





End 







, ab , maxbl k , bl 

 , aO. maxbl k , bl 

 , a 1 , maxbl k . bl 

 , a2 .maxbl k .bl 

 , a3 . maxbl k , bl 

 , b1 . maxbl k .bl 

 , b2 .maxbl k . bl 

 , b3 . maxbl k .bl 

 , thta 1 , maxbl k 

 , thta2.maxblk 

 , thta3 . maxbl k 

 f req. maxbl k . 

 dl .maxblk.bl 

 d2.maxblk.bl 

 d3.maxblk .bl 

 d4 .maxblk .bl 

 dxx , maxbl k . b 

 dyy . maxbl k ,b 



k) 



k) 



k) 



k) 



k) 



k) 



k) 



k) 



.blk) 



.blk) 



.blk) 



blk) 



k) 



k) 



k) 



k) 



tk) 



Ik) 



SUBROUTINE FFT( FR , F I . K . ICO ) 

 DIMENSION FR(4100) .FI(4100) 

 N=2**K 



IF(ICO.EO.O) GO TO 10 

 DO 8 1 = 1 .N 

 8 FI(I)=-FI(I) 

 10 Continue 

 MR=0 

 NN=N- 1 

 DO 2 M-1.NN 

 L = N 



1 L=L/2 

 IF(MR+L.GT.NN) GO TO 1 

 MR=MOD(MR.L)+L 

 IF(MR.LE M) GO TO 2 

 TR=FR(M+1 ) 



FR(M+1 )=FR(MR+1 ) 

 FR(MR+1 )-TR 

 TI=FI(M+1 ) 

 FI(M+1 )=FI(MR+1 ) 

 FI(MR+1 )=TI 



2 Continue 

 L=1 



3 IF(L.GE .N) GO TO 7 

 ISTEP=2*L 



EL = L 



DO 4 M= 1 ,L 



A = 3. 14 15926535* FLOAT( 1-M)/EL 



WR=COS(A) 



WI=SIN(A) 



DO 4 I=M,N. ISTEP 



d = I + L 



