Preserve the spectral units when creating an output header

astropy's wcs will automatically convert to reduced units, e.g. Hz, m/s, etc. for the spectral axis (https://github.com/astropy/astropy/issues/2469).  This is undesirable for reading/writing in the same unit.  This PR changes the `Header` of the cube without touching the associated WCS.