opening a ROMS zarr store

open_zarr_store_ROMS.py

import xarray as xr
from xgcm import Grid

CCS2 = xr.open_zarr(CCS2_store, consolidated=True)

# Create xgcm grid object
grid_ccs2 = Grid(CCS2, coords={'X': {'center': 'xh', 'outer': 'xq'},
                               'Y': {'center': 'yh', 'outer': 'yq'},
                               'Z': {'center': 's_rho', 'outer': 's_w'}},
                 periodic=False)

# Add depths of layers and interfaces to dataset
CCS2 = compute_depth_layers(CCS2, grid_ccs2)