For clipping, use astropy.stats.sigma_clip; use clipping in GMOSPrepare