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 c8e791f8ac2671df72fcd938ce56668fad9c2247..d467efa05d2d16c7783100541e613f5a8c036cbf 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
@@ -23,7 +23,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("bulk_mesh_node_ids")
+    point_ids.SetName("OriginalSubsurfaceNodeIDs")
     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 cb3bbfdb57754307549bd46effe317ddb7189b5f..3bd267343b661c420a38e1468708061a24b8b463 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("bulk_mesh_node_ids")
+    point_ids.SetName("OriginalSubsurfaceNodeIDs")
     point_ids.SetNumberOfComponents(1)
     point_ids.SetNumberOfTuples(N_boundary)