@keflavich - actually, even easier, if WCS only ever 'fixes' units to other directly equivalent units, you could just do in the reader:

```
cube = SpectralCube(...)
cube._spectral_unit = ...
```

This would be very nice since then we'd have no surprises due to astropy.wcs.