; PURPOSE: hc = histogram cut
;          used to remove extreme values found in volumes
;          useful to set min values in VolCorr
; INPUT:   file input, and the upper/lower bounds % of what you want 
;          to keep of a file, and threshold the rest out
;
; AUTHOR:  11/24/03 J.LeBarron
;
;-------------input----------------
FR G
[infile]5sec001	        ; File to be "cleaned"
FR G
[outfile]th5sec001	; Output,cleaned file

x50= 0.5		; Lowerbound, decimal
x51= 0.8		; Upperbound, decimal
;----------end input --------------

IF (x50.GT.1) THEN
   VM
   echo " Enter % as a decimal"
   EN 
ENDIF

IF (x51.GT.1) THEN
   VM
   echo " Enter % as a decimal"
   EN
ENDIF

FS x12,x11	; Get max/min of file
[infile]

;---Compute values for max/min
x30=(x12-x11)*x50+x11	; min
x31=(x12-x11)*x51+x11	; max

TH		; Masking/thresholding
[infile]
[outfile]
C		; Both upper and lower
x31,x30		; Upper/lower
	
EN