Encode coords for zarr router

[mpiannucci]

See #175

Previously, the coordinates of each DataArray were not encoded, only the dims were. So we couldnt plot this how it is shown below.

Tested with:

# We can access our API using fsspec's HTTPFileSystem
fs = HTTPFileSystem()

# The http mapper gives us a dict-like interface to the API
http_map = fs.get_mapper("http://0.0.0.0:8090/datasets/dbofs")

ds = xr.open_zarr(http_map, consolidated=True)
ds.temp.isel(ocean_time=0, s_rho=0).plot(x='lon_rho', y='lat_rho', cmap='jet')

[abkfenris]

It would also be nice to get a test that captures the current loss of the coords. It could also help comparing if the change in extracting attrs from the un/encoded data array is going to cause issues.

[mpiannucci]

The only thing missing is a test case. I will circle back to that, probably have to modify the create_dataset function in the test utils to make something with differing coords and dims.

[mpiannucci]

Added a simple test case, lemme know if we need something more comprehensive

[abkfenris]

That looks good to me. Thanks Matt!