|
25 | 25 |
|
26 | 26 | _NEMO_EXPECTED_COORDS = ["glamf", "gphif"] |
27 | 27 |
|
28 | | -_NEMO_DIMENSION_COORD_NAMES = ["x", "y", "time", "x", "x_center", "y", "y_center", "depth", "glamf", "gphif"] |
| 28 | +_NEMO_DIMENSION_COORD_NAMES = ["x", "y", "time", "x", "x_center", "y", "y_center", "depth", "depth_center", "glamf", "gphif"] |
29 | 29 |
|
30 | 30 | _NEMO_AXIS_VARNAMES = { |
31 | 31 | "x": "X", |
32 | 32 | "x_center": "X", |
33 | 33 | "y": "Y", |
34 | 34 | "y_center": "Y", |
35 | 35 | "depth": "Z", |
| 36 | + "depth_center": "Z", |
36 | 37 | "time": "T", |
37 | 38 | } |
38 | 39 |
|
@@ -297,7 +298,6 @@ def nemo_to_sgrid(*, fields: dict[str, xr.Dataset | xr.DataArray], coords: xr.Da |
297 | 298 | coords = coords.isel({dim: 0}) |
298 | 299 | if len(coords.dims) != 2: |
299 | 300 | raise ValueError("Expected coordsinates to be 2 dimensional") |
300 | | - |
301 | 301 | ds = xr.merge(list(fields.values()) + [coords]) |
302 | 302 | ds = _maybe_rename_variables(ds, _NEMO_VARNAMES_MAPPING) |
303 | 303 | ds = _maybe_create_depth_dim(ds) |
@@ -336,7 +336,7 @@ def nemo_to_sgrid(*, fields: dict[str, xr.Dataset | xr.DataArray], coords: xr.Da |
336 | 336 | sgrid.DimDimPadding("x_center", "x", sgrid.Padding.LOW), |
337 | 337 | sgrid.DimDimPadding("y_center", "y", sgrid.Padding.LOW), |
338 | 338 | ), |
339 | | - vertical_dimensions=(sgrid.DimDimPadding("z_center", "depth", sgrid.Padding.HIGH),), |
| 339 | + vertical_dimensions=(sgrid.DimDimPadding("depth_center", "depth", sgrid.Padding.HIGH),), |
340 | 340 | ).to_attrs(), |
341 | 341 | ) |
342 | 342 |
|
|
0 commit comments