When astrodata reads a FITS file with a wavelength axis, the units (from CUNITj) are stored but not changed. However, DRAGONS assumes that the units of wavelength axes are nm, without checking.
So we have to implement one of two things:
- Convert units at read time to nm
- Wherever a wavelength is determined from the WCS, its units have to be checked and converted to nm if necessary
Clearly 1 is easiest; 2 is a pain (which is why we didn't do this originally).