- if ix is not None, only the flux from rInfo[ix].Fcolumn is used for the spectrum
- continuum remains disk integrated, as we do not have it for each column.
With a few more modifications, this could address issue #37
Also, issue #41 would be useful to address in this branch/MR.