Please implement sigma-clipping when coadding spectra

The current code uses a 1/sigma^2 weighting, which is fine (actually optimal) if the errors are Gaussian. Unfortunately the are not, with problems like cosmic rays and bad sky subtraction leaking into the data.

Please implement an iterative sigma-clipping schema