diff --git a/examples/ScalableProblems/MittelmannBndryCntrlDiri3D.cpp b/examples/ScalableProblems/MittelmannBndryCntrlDiri3D.cpp
index 0532bab3a..9ccf7b226 100644
--- a/examples/ScalableProblems/MittelmannBndryCntrlDiri3D.cpp
+++ b/examples/ScalableProblems/MittelmannBndryCntrlDiri3D.cpp
@@ -906,7 +906,9 @@ void MittelmannBndryCntrlDiriBase3D::finalize_solution(
 
     for (Index i=0; i<=N_+1; i++) {
     for (Index j=0; j<=N_+1; j++) {
-    fprintf(fp, "y[%6d,%6d] = %15.8e\n", i, j, x[y_index(i,j)]);
+    for (Index k=0; k<=N_+1; k++) {
+    fprintf(fp, "y[%6d,%6d,%6d] = %15.8e\n", i, j, k, x[y_index(i,j,k)]);
+    }
     }
     }
 
diff --git a/examples/ScalableProblems/MittelmannBndryCntrlDiri3D_27.cpp b/examples/ScalableProblems/MittelmannBndryCntrlDiri3D_27.cpp
index a69bdf53a..e7fa774eb 100644
--- a/examples/ScalableProblems/MittelmannBndryCntrlDiri3D_27.cpp
+++ b/examples/ScalableProblems/MittelmannBndryCntrlDiri3D_27.cpp
@@ -1074,7 +1074,9 @@ void MittelmannBndryCntrlDiriBase3D_27::finalize_solution(
 
     for (Index i=0; i<=N_+1; i++) {
     for (Index j=0; j<=N_+1; j++) {
-    fprintf(fp, "y[%6d,%6d] = %15.8e\n", i, j, x[y_index(i,j)]);
+    for (Index k=0; k<=N_+1; k++) {
+    fprintf(fp, "y[%6d,%6d,%6d] = %15.8e\n", i, j, k, x[y_index(i,j,k)]);
+    }
     }
     }
 
diff --git a/examples/ScalableProblems/MittelmannBndryCntrlDiri3Dsin.cpp b/examples/ScalableProblems/MittelmannBndryCntrlDiri3Dsin.cpp
index f3fbb9101..201a6e913 100644
--- a/examples/ScalableProblems/MittelmannBndryCntrlDiri3Dsin.cpp
+++ b/examples/ScalableProblems/MittelmannBndryCntrlDiri3Dsin.cpp
@@ -1001,7 +1001,9 @@ void MittelmannBndryCntrlDiriBase3Dsin::finalize_solution(
 
     for (Index i=0; i<=N_+1; i++) {
     for (Index j=0; j<=N_+1; j++) {
-    fprintf(fp, "y[%6d,%6d] = %15.8e\n", i, j, x[y_index(i,j)]);
+    for (Index k=0; k<=N_+1; k++) {
+    fprintf(fp, "y[%6d,%6d,%6d] = %15.8e\n", i, j, k, x[y_index(i,j,k)]);
+    }
     }
     }