Adding unit test for vertical 1D field

Note that this unit test currently breaks

Author: erikvansebille

tests/test_xgrid.py

@@ -4,14 +4,17 @@
 import numpy as np
 import pytest
 import xarray as xr
+import xgcm
 from numpy.testing import assert_allclose
 
+from parcels import Field
 from parcels._core.index_search import (
     LEFT_OUT_OF_BOUNDS,
     RIGHT_OUT_OF_BOUNDS,
     _search_1d_array,
 )
 from parcels._core.xgrid import (
+    _DEFAULT_XGCM_KWARGS,
     XGrid,
     _transpose_xfield_data_to_tzyx,
 )
@@ -134,6 +137,15 @@ def test_invalid_depth():
         XGrid.from_dataset(ds)
 
 
+def test_vertical1D_field():
+    nz = 11
+    ds = xr.Dataset({"z1d": (["depth"], np.linspace(0, 10, nz))}, coords={"depth": np.linspace(0, 1, nz)})
+    grid = XGrid(xgcm.Grid(ds, **_DEFAULT_XGCM_KWARGS))
+    field = Field("z1d", ds["z1d"], grid)
+
+    assert field.eval(np.timedelta64(0, "s"), 0.45, 0, 0) == 4.5
+
+
 @pytest.mark.parametrize(
     "ds",
     [