Pyvista: use point_data instead of point_arrays (resolve warning from 0.32)

Author: MuellerSeb

examples/00_misc/04_herten.py

@@ -70,7 +70,7 @@ def generate_transmissivity():
 
     print("Loading Herten data with pyvista")
     mesh = pv.read(VTK_PATH)
-    herten = mesh.point_arrays["facies"].reshape(mesh.dimensions, order="F")
+    herten = mesh.point_data["facies"].reshape(mesh.dimensions, order="F")
     # conductivity values per fazies from the supplementary data
     cond = 1e-4 * np.array(
         [2.5, 2.3, 0.61, 260, 1300, 950, 0.43, 0.006, 23, 1.4]

gstools/tools/export.py

@@ -81,7 +81,7 @@ def to_vtk_structured(pos, fields):  # pragma: no cover
     x, y, z, fields = _vtk_structured_helper(pos=pos, fields=fields)
     if pv is not None:
         grid = pv.RectilinearGrid(x, y, z)
-        grid.point_arrays.update(fields)
+        grid.point_data.update(fields)
     else:
         raise ImportError("Please install PyVista to create VTK datasets.")
     return grid
@@ -155,7 +155,7 @@ def to_vtk_unstructured(pos, fields):  # pragma: no cover
     x, y, z, fields = _vtk_unstructured_helper(pos=pos, fields=fields)
     if pv is not None:
         grid = pv.PolyData(np.c_[x, y, z]).cast_to_unstructured_grid()
-        grid.point_arrays.update(fields)
+        grid.point_data.update(fields)
     else:
         raise ImportError("Please install PyVista to create VTK datasets.")
     return grid

setup.cfg

@@ -74,7 +74,7 @@ doc =
     sphinx-rtd-theme>=1
 plotting =
     matplotlib>=3
-    pyvista>=0.29
+    pyvista>=0.32
 rust =
     gstools_core>=0.2.0,<1
 test =