The moment 0 should apply the `BUNIT` of the file as well since doesn't have the weights in the denominator.  Thus, I'd expect `cube.moment(axis=0,order=0)` to return a Quantity with units of `K km/s` or `K GHz` or `Jy/Beam km/s` and their many, many equivalents.