Clarify what is happpening during saturations
Spit out a timeseries when we report a channel saturates to see if there are one or many points above the saturation threshold. It's unclear what needs to be passed in the command line argument --ignore-saturation to enter the ignore-saturated=False logic. Plot spectra for saturated channels regardless of their saturation status. Redefine "saturated" as fraction of points above saturation threshold?