Support encoding in to_xarray : notably unlimited_dims

[pp-mo]

Direct loading via xr.open_dataset() sets ds.encoding['unlimited_dims']
But currently, ncdata.xarray.to_xarray does not do this, meaning that unlimited dims get lost.

( Maybe surprisingly though, this does not affect Dataset.identical() checking : encodings do not count there )