https://nbviewer.org/gist/rsignell-usgs/07bc53e5ef7768ccc345c7d03fcb85b8
import xbitinfo as xb
# Explore bit information in all dimensions
bi = xb.get_bitinformation(ds)
fig = xb.plot_bitinformation(bi.sel(dim="latitude"))
kb = xb.get_keepbits(bi)
xb.xr_bitround(ds, 10).to_compressed_zarr("file.zarr")