CALL AXES(0.,0.,12HDISTANCE(NM),-12,XL,0.,1.,0.,DV,-1) 



EV=AZAVE-(3.0/OSCALE) 



DV=1.0/0SCALE 



CALL AXES(0.,-3.,13H0RIGINAL DATA, !#,(t. ,().,!. ,EV,DV,1) 



CALL SYMB0L(XL+.75,-.25,.5,12HINPUT SIGNAL, 0.0, 12) 



IF (C2MAX.lt, XL) C2MAX=XL 



CALL PL0T(0.0,(AZ(1)-AZAVE)*0SCALE,3) 



DO 26 I=1,JCT 



AJ=I-1 



C=AJ*.02 



D=(AZ(I)-AZAVE)*OSCALE 



CALL PL0T(C,D,2) 



26 CONTINUE 



DO 27 I=1,JCT 



27 Z(I)=AZ(I) 

 IFdPLOT.EQ.DCALL PLOT(0. 0,1. 0,-3) 

 DO 237 IRUN=1,10 



XRUN=IRUN 



ENVMAX=0.0 



IFdPLOT.EQ.l.AND.IBANO.EQ.DCALL PLOT(0. 0.4. 0,-3) 

 C BAND PASS FILTER DATA AT SPECIFIED INTERVAL 



CALL FILTER(AZ,XT,FRL0W(IRUN),FILSLP,NUMFIL,1,WGT) 



CALL FILTER(AZ,JCT,FRHIGH(IRUN),FILSLP,NUMFIL,0,W6T) 



IF(IPLOT.EQ.l) CALL PL0T(0.0,0.0,3) 



WGTMN=0.0 



XN0RM=((FRHIGH(IRUN)+.2)-FRL0W(IRUN))*(l./GRIDKM) 

 C XN0RM=1. 



NTWICE=2*NUMFIL 



DO 31 I=NTWICE,JCT-NTWICE 

 C — 1 



IF(IPLOT.NE.l.OR.IBAND.NE.l) GO TO 31 

 C — 1 



AJ=I-1 



C=AJ*.02 



B=(AZ(I)/XNORM)*ESCALE 



IF(I.EQ.NTWICE) CALL PL0T(C,8,3) 



CALL PL0T(C,B,2) 



31 AZ(I)=(AZ(I))/XNORM 

 C — 1 



IF(IPLOT.NE.l.OR.IBAND.NE.l) GO TO 32 

 C — 1 



EV=-1.5/ESCALE 



DV=1./ESCALE 



CALL PL0T(0.0,0.0,3) 



CALL AXES(0.,0.,1H ,1,XL,0. ,1. ,0. ,0. ,-2) 



CALL AXES(XL,-1.5,13HFILTERED DATA,-13,3.,90. ,1.5,EV,DV,0) 



CALL SYMB0L(XL+.75,-.25,.5,9HBAND PASS, 0.0, 9) 



CALL NUMBER(XL+5.75,-.25,.5,XRUN,0.0,-l) 



CALL PL0T(0.0,0.0,3) 



CALL NEWPEN(2) 

 C COMPUTE APPROXIMATION TO ENVELOPE OF BAND PASSED DATA 



32 CALL ENVEL(AZ,XT) 



152 



