issue_comments

@jhamman thanks I'll add to the discussion there and close this issue.

Is there a way of appending an xarray dataset onto an existing zarr array?

@NickMortimer - No, not yet. It has been proposed already though in #2022.

@spencerkclark yes that helps very much and a great example of how to answer a question! I'm learning so much from this group. Is there a way of appending an xarray dataset onto an existing zarr array? That's why I've been accessing direct through zarr, what I'm trying to do is build a zarr file of all the Argo float profiles and add new ones as they arrive.

When writing datetime objects to disk, xarray encodes them following CF conventions (converting them to numerical values quantifying some units of time since a given reference date). This is to support other backends which do not support writing datetime64 objects directly to disk.

For accurate roundtripping of Datasets, xarray also includes logic to automatically decode datetimes stored following CF conventions. For that reason, instead of loading in the raw zarr store using zarr's open function, I recommend using xarray's open_zarr function, which will automatically decode the CF-encoded values to datetime64 objects.

See the following example: ``` In [1]: import numpy as np

In [2]: import xarray as xr

In [3]: da = xr.DataArray(np.datetime64('2000-01-01'), name='date')

In [4]: da Out[4]: <xarray.DataArray 'date' ()> array('2000-01-01T00:00:00.000000000', dtype='datetime64[ns]')

In [5]: da.to_dataset().to_zarr('example.zarr') Out[5]: <xarray.backends.zarr.ZarrStore at 0x1109ca190>

In [6]: ds = xr.open_zarr('example.zarr')

In [7]: ds.date Out[7]: <xarray.DataArray 'date' ()> array('2000-01-01T00:00:00.000000000', dtype='datetime64[ns]') Note if we open the zarr store directly, we'll find that the date was encoded with the integer `0` and appropriate units and calendar attributes: In [1]: import zarr

In [2]: z = zarr.open('example.zarr')

In [3]: z['date'][...] Out[3]: array(0)

In [4]: z['date'].attrs['units'] Out[4]: u'days since 2000-01-01 00:00:00'

In [5]: z['date'].attrs['calendar'] Out[5]: u'proleptic_gregorian' ``` Unlike xarray, zarr does not include logic for automatically decoding CF-encoded datetimes. I hope that helps.