@@ -0,0 +1,816 @@
+from parcels import Grid%0Aimport numpy as np%0Aimport pytest%0A%0A%0A@pytest.mark.parametrize('xdim', %5B100, 200%5D)%0A@pytest.mark.parametrize('ydim', %5B100, 200%5D)%0Adef test_grid_from_data(xdim, ydim):%0A    lon = np.linspace(0., 1., xdim, dtype=np.float32)%0A    lat = np.linspace(0., 1., ydim, dtype=np.float32)%0A    depth = np.zeros(1, dtype=np.float32)%0A    time = np.zeros(1, dtype=np.float64)%0A    u, v = np.meshgrid(lon, lat)%0A%0A    grid = Grid.from_data(u, lon, lat, v, lon, lat, depth, time)%0A    u_t = np.transpose(u).reshape((lat.size, lon.size))%0A    v_t = np.transpose(v).reshape((lat.size, lon.size))%0A    assert len(grid.U.data.shape) == 3  # Will be 4 once we use depth%0A    assert len(grid.V.data.shape) == 3%0A    assert np.allclose(grid.U.data%5B0, :%5D, u_t, rtol=1e-12)%0A    assert np.allclose(grid.V.data%5B0, :%5D, v_t, rtol=1e-12)%0A

Add a baseline test for creating grids from data