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



Y=ABS(SMPH(J)/30.0) 



IF(IPHA.GT.0)Y=Y*30.0/APHA 



IBB=4 



IF(Y.LT.8.0.AND.SMPH(J).GE.0.0) GO TO 550 



IF(Y.LT.8.0.AND.SMPH(J).LT.0.0) IBB=2 



IF(Y.GE.8.0) Y=8.5 



CALL SYMB0L(X,Y,.05,IBB,0.,-2) 



GO TO 574 

 550 CALL PL0T(X,Y,2) 

 574 CONTINUE 

 C PLOT PHASE AXIS 



579 IF(JCODE.LT.O) GO TO 581 

 C CALL NEWPEN(l) 



IF(IPHA.GT.O) GO TO 547 



CALL AXES(0.,0.,14HABS PHASE(DEG) ,14,6. ,90. ,1. ,0. ,30. ,-1) 



GO TO 558 

 547 CALL AXES(0.,0.,20HABS PHASE(DATA INT. ),20,8. ,90. ,1. ,0. ,APHA,-1) 

 558 CALL AXES(0.,0.,ITG,+30,12.2,0.,3.0125,10.,0.,-1) 



CALL SYMB0L(2.0,-.2,.14,33HTRIANGLE INDICATES NEGATIVE PHASE, 0.0, 

 1 33) 

 581 CONTINUE 



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

 C COMPUTE AMPLITUDE SPECTRUM OF SMOOTHING FILTER 



DELX= 0.01 



DO 239 J=l,51 



AJ=J-1 



XJ=AJ*DELX 



SUM=0.0 



DO 240 1=2, NA 



AI=I-1 



240 SUM=SUM+2*B0ATA(I)* C0S(2.*PI*AI*XJ ) 

 DATA(J)= SUM+BDATA(1) 



239 CDATA(J)=XJ 



WRITE(I0UT,241) CUT(2) ,H(2),N(2) 



241 FORMATC AMP.SPECT.OF SMOOTHING FILTER CUTOFF=' ,F5.4, 

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



IFdLIST.NE.'l) WRlfE(I0UT,242)(CDATA(I),DATA(I),I=l,51) 



242 F0RMAT(8(F5.2,F8.3)) 

 331 WRITE(I0UT,243) JSET 



243 FORMATC END OF DATA SET NO. ',13) 

 C END FILE 12 



C END FILE 12 



IF(JSET.EQ.NSETS) GO TO 244 



JSET=JSET+1 



IF(JCOOE.LT.O) GO TO 245 

 C CALL PLOT (15. 0,0. 0,-3) 



GO TO 245 



244 IF(JCOOE.LT.O) GO TO 487 

 CALL PLOT (0.0, 0.0, 999) 



487 END FILE 14 

 END FILE 14 

 STOP 

 END 



196 



