revert AvsC grid accuracy guide

Author: reint-fischer

…amples/tutorial_peninsula_AvsCgrid.ipynb …les_v3/tutorial_peninsula_AvsCgrid.ipynbdocs/user_guide/examples/tutorial_peninsula_AvsCgrid.ipynb renamed to docs/user_guide/examples_v3/tutorial_peninsula_AvsCgrid.ipynb docs/user_guide/examples/tutorial_peninsula_AvsCgrid.ipynb renamed to docs/user_guide/examples_v3/tutorial_peninsula_AvsCgrid.ipynb

@@ -27,7 +27,6 @@
     "import numpy as np\n",
     "import scipy\n",
     "import xarray as xr\n",
-    "import xgcm\n",
     "from matplotlib.lines import Line2D\n",
     "\n",
     "import parcels"
@@ -56,29 +55,25 @@
     "    R = 0.32 * domainsizeX / 2\n",
     "\n",
     "    x, y = np.meshgrid(lon, lat, sparse=True, indexing=\"xy\")\n",
-    "    P_data = u0 * R**2 * y / ((x - x0) ** 2 + y**2) - u0 * y\n",
-    "    landpoints = P_data >= 0.0\n",
-    "    P_data[landpoints] = 0.0\n",
+    "    P = u0 * R**2 * y / ((x - x0) ** 2 + y**2) - u0 * y\n",
+    "    landpoints = P >= 0.0\n",
+    "    P[landpoints] = 0.0\n",
     "\n",
     "    if grid_type == \"A\":\n",
-    "        U_data = u0 - u0 * R**2 * ((x - x0) ** 2 - y**2) / (((x - x0) ** 2 + y**2) ** 2)\n",
-    "        V_data = -2 * u0 * R**2 * ((x - x0) * y) / (((x - x0) ** 2 + y**2) ** 2)\n",
-    "        U_data[landpoints] = 0.0\n",
-    "        V_data[landpoints] = 0.0\n",
+    "        U = u0 - u0 * R**2 * ((x - x0) ** 2 - y**2) / (((x - x0) ** 2 + y**2) ** 2)\n",
+    "        V = -2 * u0 * R**2 * ((x - x0) * y) / (((x - x0) ** 2 + y**2) ** 2)\n",
+    "        U[landpoints] = 0.0\n",
+    "        V[landpoints] = 0.0\n",
     "    elif grid_type == \"C\":\n",
-    "        U_data = np.zeros(P.shape)\n",
-    "        V_data = np.zeros(P.shape)\n",
-    "        V_data[:, 1:] = (P_data[:, 1:] - P_data[:, :-1]) / (lon[1] - lon[0])\n",
-    "        U_data[1:, :] = -(P_data[1:, :] - P_data[:-1, :]) / (lat[1] - lat[0])\n",
+    "        U = np.zeros(P.shape)\n",
+    "        V = np.zeros(P.shape)\n",
+    "        V[:, 1:] = (P[:, 1:] - P[:, :-1]) / (lon[1] - lon[0])\n",
+    "        U[1:, :] = -(P[1:, :] - P[:-1, :]) / (lat[1] - lat[0])\n",
     "\n",
+    "    data = {\"U\": U, \"V\": V, \"P\": P}\n",
     "    dimensions = {\"lon\": lon, \"lat\": lat}\n",
-    "    grid = parcels.XGrid(xgcm.grid())\n",
-    "    U = parcels.Field(\"U\", U_data, grid=, interp_method=)\n",
-    "    V = parcels.Field(\"V\", V_data, grid=, interp_method=)\n",
-    "    UV = parcels.VectorField(\"UV\", U, V, vector_interp_method=)\n",
-    "    P = parcels.Field(\"P\", P_data, grid=, interp_method=)\n",
     "\n",
-    "    fieldset = parcels.FieldSet([U, V, P])\n",
+    "    fieldset = parcels.FieldSet.from_data(data, dimensions, mesh=\"flat\")\n",
     "    fieldset.U.interp_method = interp_method\n",
     "    fieldset.V.interp_method = interp_method\n",
     "    return fieldset"
@@ -97,8 +92,8 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "def SampleP(particles, fieldset):\n",
-    "    particles.p = fieldset.P[particles]\n",
+    "def SampleP(particle, fieldset, time):\n",
+    "    particle.p = fieldset.P[time, particle.lon, particle.lat, particle.depth]\n",
     "\n",
     "\n",
     "def DeleteParticle(particle, fieldset, time):\n",
@@ -363,7 +358,7 @@
  ],
  "metadata": {
   "kernelspec": {
-   "display_name": "test-notebooks",
+   "display_name": "parcels",
    "language": "python",
    "name": "python3"
   },
@@ -377,7 +372,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.11.0"
+   "version": "3.12.3"
   }
  },
  "nbformat": 4,