Add support for 2D spectra in `ApplyFluxCal`

[migueldvb]

Sometimes it is helpful to apply the sensitivity function from a standard target to a 2D spectrum with one axis being wavelength and the spatial direction to study extended sources.