A simple 1D sigma-clipping routine that returns a boolean array of indices to non-clipped data
sigma_clip()