diff --git a/pre/BCs/inhomogeneous-mass-flux/prepare-inhomogeneous-mass-flux-BC.py b/pre/BCs/inhomogeneous-mass-flux/prepare-inhomogeneous-mass-flux-BC.py
index d23de9fc8cbe3155f008154106181007d95db6dd..659778babe5f7ad38bdaa327854a99758a2590dc 100755
--- a/pre/BCs/inhomogeneous-mass-flux/prepare-inhomogeneous-mass-flux-BC.py
+++ b/pre/BCs/inhomogeneous-mass-flux/prepare-inhomogeneous-mass-flux-BC.py
@@ -31,7 +31,7 @@ def do_enumerate_points():
     # point_ids = vtk.vtkIdTypeArray()  ## TODO why does that not work?
     # point_ids = vtk.vtkLongArray()
     point_ids = vtk.vtkUnsignedLongArray()
-    point_ids.SetName("OriginalSubsurfaceNodeIDs")
+    point_ids.SetName("bulk_node_ids")
     point_ids.SetNumberOfComponents(1)
     N = data.GetPoints().GetNumberOfPoints()
     point_ids.SetNumberOfTuples(N)
diff --git a/pre/BCs/map-boundary-mesh-to-bulk-mesh/map-boundary-to-bulk.py b/pre/BCs/map-boundary-mesh-to-bulk-mesh/map-boundary-to-bulk.py
index 3bd267343b661c420a38e1468708061a24b8b463..95907980096494f9605589b4093b270672fab498 100755
--- a/pre/BCs/map-boundary-mesh-to-bulk-mesh/map-boundary-to-bulk.py
+++ b/pre/BCs/map-boundary-mesh-to-bulk-mesh/map-boundary-to-bulk.py
@@ -32,7 +32,7 @@ def do_mapping():
     # point_ids = vtk.vtkIdTypeArray()  ## TODO why does that not work?
     # point_ids = vtk.vtkLongArray()
     point_ids = vtk.vtkUnsignedLongArray()
-    point_ids.SetName("OriginalSubsurfaceNodeIDs")
+    point_ids.SetName("bulk_node_ids")
     point_ids.SetNumberOfComponents(1)
     point_ids.SetNumberOfTuples(N_boundary)