diff --git a/HEU_MET_FAST_003/HEU_MET_FAST_003.py b/HEU_MET_FAST_003/HEU_MET_FAST_003.py
index 2660938..79efe11 100644
--- a/HEU_MET_FAST_003/HEU_MET_FAST_003.py
+++ b/HEU_MET_FAST_003/HEU_MET_FAST_003.py
@@ -132,7 +132,7 @@
     fflist = phys.GetScalarFieldFunctionList(only_scalar_flux=False)
     vtk_basename = "HEU_MET_FAST_003"
     # export only the flux of group g (first []), moment 0 (second [])
-    FieldFunctionGridBased.ExportMultipleToVTK(
+    FieldFunctionGridBased.ExportMultipleToPVTU(
         [fflist[g][0] for g in range(num_groups)],
         vtk_basename
     )
diff --git a/OpenSn_Logo_CAD/opensn.py b/OpenSn_Logo_CAD/opensn.py
index 7fdf6ad..279b90a 100644
--- a/OpenSn_Logo_CAD/opensn.py
+++ b/OpenSn_Logo_CAD/opensn.py
@@ -70,4 +70,4 @@
 
     # Export results to VTK
     fflist = phys.GetScalarFieldFunctionList()
-    FieldFunctionGridBased.ExportMultipleToVTK(fflist, "flux")
+    FieldFunctionGridBased.ExportMultipleToPVTU(fflist, "flux")
diff --git a/Six_1g_spherical_benchmarks/Problem_2.py b/Six_1g_spherical_benchmarks/Problem_2.py
index 3905d0f..95e4111 100644
--- a/Six_1g_spherical_benchmarks/Problem_2.py
+++ b/Six_1g_spherical_benchmarks/Problem_2.py
@@ -90,7 +90,7 @@
     """
         vtk_basename = "Problem_2_"
         # export only the scalar flux of group g
-        FieldFunctionGridBased.ExportMultipleToVTK(
+        FieldFunctionGridBased.ExportMultipleToPVTU(
             [fflist[g] for g in range(num_groups)],
             vtk_basename
         )
diff --git a/Urban_Source/urban_source.py b/Urban_Source/urban_source.py
index 0a3a592..d2230d3 100644
--- a/Urban_Source/urban_source.py
+++ b/Urban_Source/urban_source.py
@@ -73,4 +73,4 @@
 
     # Export results to VTK
     fflist = phys.GetScalarFieldFunctionList()
-    FieldFunctionGridBased.ExportMultipleToVTK(fflist, "flux")
+    FieldFunctionGridBased.ExportMultipleToPVTU(fflist, "flux")