Calculate integrated intensity

[low-sky]

We'd like to calculate the integrated intensity of W11, W22 and W33 for the model

• Get a model spectrum out (??)
• np.sum(ModelSpectrum) * size_of_channel
• Then sum up significant parts of the spectrum.

[low-sky]

Good so far!

• Calculate the W11 empirically over the observed data where the model is > than, say 1e-6
• Calculate the typical root-mean-squared error in a spectrum using np.nanstd or np.std ignoring negative values or the median absolute deviation code here. This gives the RMS in a channel which propagates to W11 empirical.
• Store output data in an astropy Table.
• Identify bad fits based differences between W11_theory and W11_empirical.