BUNIT exponent likely wrong in FluxCalibrate.py

[BMazzilliCiraulo]

Hi,

The header of the flux-calibrated cubes (i.e. icubes.fits) says:
BUNIT = '1e+16 erg / (Angstrom cm2 s)'

In FluxCalibrate.py:

flam16_u = 1.e16 * u.erg / (u.angstrom * u.cm ** 2 * u.s)
self.action.args.ccddata.unit = flam16_u
self.action.args.ccddata.uncertainty.unit = flam16_u

I think it should be 1e-16 instead of 1e+16. This would make the numerical values physically reasonable.

A pull request can be prepared if desired.

Thanks!