diff --git a/ntrfc/meshquality/nondimensionals.py b/ntrfc/meshquality/nondimensionals.py
index a53f16ae52e5a5e06fe1a71aecb90f98d79a286e..3bb72bc76d33cd9ca55e9514ec041cd1e9592ae2 100644
--- a/ntrfc/meshquality/nondimensionals.py
+++ b/ntrfc/meshquality/nondimensionals.py
@@ -45,7 +45,7 @@ def cellDirections(cellUMean, wallNorm):
 def closestWallNormalPoint(point, surfacenormals):
     surfacepoint_id = surfacenormals.find_closest_point(point)
     wallpoint = surfacenormals.points[surfacepoint_id]
-    return surfacenormals.point_arrays["Normals"][surfacepoint_id], wallpoint - point
+    return surfacenormals.point_data["Normals"][surfacepoint_id], wallpoint - point
 
 
 def calcWallNormalVectors(surfaceMesh, volmesh):
diff --git a/ntrfc/turbo/probegeneration.py b/ntrfc/turbo/probegeneration.py
index ba0b2c9213e66740c7eb177911858bc9459a0254..fcfac0ded838c01acd1dba577c7d982435aefd71 100644
--- a/ntrfc/turbo/probegeneration.py
+++ b/ntrfc/turbo/probegeneration.py
@@ -70,7 +70,7 @@ def create_stagnationpointprobes(length, nop, sortedpoly, ind_vk, u_inlet, midsp
     vk_point = sortedpoly.points[ind_vk]
     angle = vecAngle(u_inlet, np.array([1, 0, 0]) * 180 / np.pi)
     stagnationLine = pv.Line((0, 0, 0), (-length, 0, 0), nop - 1)
-    stagnationLine.rotate_z(angle)
+    stagnationLine.rotate_z(angle,inplace=False)
     stagnationLine.translate(vk_point)
     x_probes = stagnationLine.points[::, 0]
     y_probes = stagnationLine.points[::, 1]
diff --git a/tests/turbo/test_ntrfc_probegeneration.py b/tests/turbo/test_ntrfc_probegeneration.py
index e8263df8664d2a93527c0a5bc7e99cfd7369b52d..6a7d8d6cbea8eba7a423fbed89af9701f8fc98a2 100644
--- a/tests/turbo/test_ntrfc_probegeneration.py
+++ b/tests/turbo/test_ntrfc_probegeneration.py
@@ -30,7 +30,7 @@ def test_create_midpassageprobes():
     res = 240
     xs, ys = naca(naca_code, res, half_cosine_spacing=False)
     sorted_poly = pv.PolyData(np.stack([xs, ys, np.zeros(len(xs))]).T)
-    sorted_poly.rotate_z(angle)
+    sorted_poly.rotate_z(angle,inplace=False)
     sorted_extracted_poly, psPoly, ssPoly, ind_vk, ind_hk, midsPoly, beta_leading, beta_trailing, camber_angle = extract_geo_paras(
         sorted_poly, 1)
     nop = 40