Preserve spectral unit requested

At the moment, if I run ``with_spectral_unit(u.km/u.s)``, it returns m/s due to the limitations of astropy.wcs. However, I think this will confuse users, so how about storing internally the spectral unit required, and converting on-the-fly when accessing ``spectral_axis``?