- 22^ 



48050 

 48100 

 48150 

 48200 

 48250 

 48300 

 48350 

 48400 

 48450 

 48500 

 48550 

 48600 

 48650 

 48700 

 48750 

 48800 

 48850 

 48900 

 48950 

 49000 

 49050 

 49100 

 49150 

 49200 

 49250 

 49300 

 49350 

 49400 

 49450 

 49500 

 49550 

 49600 

 49650 

 49700 

 49750 

 49800 

 49850 

 49900 

 49950 

 50000 

 50050 

 50100 

 50150 

 50200 

 50250 

 50300 

 50350 

 50400 

 50450 

 50500 

 50550 

 50600 

 50650 

 50700 

 50750 

 50800 

 50850 

 50900 

 50950 

 51000 



Return 

 End 



Subrout 1 ne SCAN( X , N , NGAGE , A VG . DX INCH ) 



Dimension X(4100). HIGH( 1000) 

 Real LOW(IOOO), LOLIM 

 Integer HISTOP 



C Calculates the mean and std. deviat ion (N-1 weighting) 

 C for the record. 



RUN =0.0 

 SUM = O.O 

 DO 1 19 1 = 1 .N 

 1 19 SUM = X( i ) + sum 

 AVG = SUM/N 



00 130 1=1 ,N 

 130 RUN = (X( 1 )-AVG)**2 + RUN 

 DEV = sqrt(RUN/(N-1 )) 



C Reduce those points beyond 3 std . dev i at ions to 

 C 3 std. deviations beyond record mean or the 

 C average of neighboring points. 



HILIM = AVG + 3.0*DEV 



LOLIM = AVG - 3.0*DEV 



Khi = O 



Klo = 



X( 1 ) = AVG 



DO 300 I=2.N 



IF(X( 1 ) .LE .HILIM and. X( i ). GT . LOLIM) GO TO 300 

 IF(X( i ) .GT.HILIM) GO TO 150 

 IF(X( 1 ) .LT. LOLIM) Klo = Klo + 1 

 LOW(Klo) = X( i ) 

 IF(X( i-i ) . 1 t .LOLIM. or. X( 1+1 ) . It LOLIM) 



X( 1 ) = LOLIM 

 IF(X(1-1 ) . ge. LOLIM. or. X( 1 + 1 ) .ge. LOLIM) 



X(i) = (X(i-1)+X(1+1))/2.0 

 GO TO 300 

 Khi = Khi + 1 

 HIGH(Kh1 ) = X( i ) 

 IF(X( 1-1) .gt .HILIM.or.XC i + 1 ) .gt.HILIM) 



X( 1 ) = HILIM 

 IF(X( 1- D.le. HILIM. or. X( i+1 ).le. HILIM) 

 X( 1 ) = (X( i-1 ) + X( i + 1))/2.0 

 Cont inue 



Print out truncated points 



