WRITE(I0UT,570) 

 570 FORMATC ROUGH PHASE SP£CTRUM(DEG) ' ) 

 GO TO 542 



541 WRITE(I0UT,543) 



543 FORMATC ROUGH PHASE SPECTRUM(DATA INTERVALS)') 



542 IF(ILIST.NE.l) WRITE(I0UT,81)(PHASE( J),J = 1,JX ) 

 558 IF(K0UNT-2)202,203,204 



202 WRITE (lOUT, 75) 



75 F0RMAT(28H AMP.SPECT.OF ORIGINAL DATA ) 



77 IF(ILIST.NE.I) WRITE(I0UT,81)(BDATA( J),J=1,JX ) 

 81 FORMAT (lOE 11. 4) 

 GO TO 206 



203 WRITE (lOUT, 76) 



76 F0RMAT(30H AMP.SPECT.OF PREWHITENED DATA) 

 GO TO 77 



204 WRITE(I0UT,91) CUT(1),H(1),N(1) 



91 FORMATC SPECT.OF PREWHITENING FILTER WITH PHASE REVERSALS CUTOF 

 1F=',F5.4,' H=',F5.4,' N=',I3) 



00 92 J=1,JX 



92 BDATA(J) =BDATA( J)/ANORM 



IF(ILIST.NE.l) WRITE(I0UT,81)(BDATA(J),J=1,JX) 

 C CORRECT SPECTRUM OF FILTERED DATA FOR PREWHITENING FILTER 

 00 93 J=1,JX 

 IF(ABS(BDATA(J)).LT.l.E-20)BDATA(J)=BDATA(J)+l.E-20 



93 8DATA(J)=ABS(CDATA(J)/BDATA(J)) 

 WRITE(I0UT,94) 



94 FORMATC FINAL ROUGH CORRECTED SPECTRUM WITHOUT PHASE REVERSALS') 

 IF(ILIST.NE.l) WRITE(I0UT,81)(BDATA(J).J=1,JX) 



60 TO 206 



95 CONTINUE 



C SMOOTH FINAL ROUGH CORRECTED SPECTRUM 

 C COMPUTE LOW PASS WEIGHTS 



K=2 



IF(N(2).EQ.O) GO TO 330 

 C STORE ROUGH SPECTRUM IN DATA 



DO 97 J=1,JX 

 97 DATA(J)=BDATA(J) 



GO TO 98 

 83 00 99 1=1, NA 

 99 BDATA(I)=BDATA(I)+DELTA/X 

 C CONVOLVE LOW PASS SMOOTHING WEIGHTS WITH ROUGH SPECTRUM 



NB=JX-N(K) 



GO TO 235 

 206 IF(JCOOE.EQ.l.AND.N(2).GT.O) GO TO 311 



IF(JCOOE.LT.O) GO TO 311 



X=AL0G10 ( DELK/X0RG)*3.0125 



TPP=BDATA(2)/Y0RG 



IF(TPP.LT.l.E-20) TPP=YORG 



Y=-0.5+AL0G10(TPP)*l. 35714 



CALL PL0T(X,Y,3) 



DO 42 J=3,JX-1 



XJ=J-1 



X=AL0G10((XJ*DELK)/X0RG)*3.0125 



189 



