pro setup_MRL, deconv_info, conv_cut

;Frank Varosi NASA/GSFC 1992.

	Lag = round_off( deconv_info.psf_stats.fwhm ) > 1
	text=""
	read," select residual auto-correlation Lag max, min  (default=" + $
					strtrim( Lag,2 ) + " to 2) ",text
	if strlen( text ) GT 0 then Lag = fix( get_words( text ) ) > 0
	deconv_info.Lag_chisq = Lag(0)
	if N_elements( Lag ) GT 1 then deconv_info.Lag_min = Lag(1) $
				  else deconv_info.Lag_min = 2
	deconv_info.Lag_chisq = deconv_info.Lag_chisq > deconv_info.Lag_min
	print," auto-correlation Lag range (max,min) =", $
				deconv_info.Lag_chisq, deconv_info.Lag_min

	sigma = deconv_info.image_stats.sigma_noise
	read," enter st.dev. (sigma) of noise (default=" + $
					strtrim( sigma,2 ) + ") : ",text
	if strlen( text ) GT 0 then sigma = float( text ) > 0
	deconv_info.image_stats.sigma_noise = sigma
	print," using st.dev. of noise =",sigma

	conv_cut = sigma/deconv_info.image_stats.rms
	read," convergence ratio cutoff (default=" + $
					strtrim( conv_cut,2 ) + ") ? ",text
	if strlen( text ) GT 0 then conv_cut = float( text ) > 0
	print," using convergence ratio cutoff =",conv_cut

;	epsilon = 1
;	print," 1 = use differentiable positivity function"
;	print," 0 = just truncate image > 0"
;	read," select type of positivity (default=1) : ",text
;	if strlen( text ) GT 0 then epsilon = fix( text ) > 0
;	if (epsilon LE 0) then print," truncating image > 0" $
;			  else print," using differentiable positivity"

	deconv_info.positivity_eps = 1.e-7
return
end