diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index d560e9ea415736974fbd8766c117c7d79af0531e..fd37631381bf05df7441c4bc3b0d332436086cf6 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -435,15 +435,17 @@ int main(int argc, char *argv[]) {
     std::cout << "sup |u1 - u4| = " << diff4.infinity_norm() << ", "
               << "|u1 - u4| = " << diff4.two_norm() << std::endl;
 
-    // Print displacement on frictional boundary
-    boost::format const formatter("u1[%02d] = %+3e, "
-                                  "%|40t|u2[%02d] = %+3e, "
-                                  "%|80t|u3[%02d] = %+3e, "
-                                  "%|120t|u4[%02d] = %+3e");
-    for (size_t i = 0; i < frictionalNodes.size(); ++i)
-      if (frictionalNodes[i][0])
-        std::cout << boost::format(formatter) % i % u1[i] % i % u2[i] % i %
-                         u3[i] % i % u4[i] << std::endl;
+    if (parset.get<bool>("printFrictionalBoundary")) {
+      // Print displacement on frictional boundary
+      boost::format const formatter("u1[%02d] = %+3e, "
+                                    "%|40t|u2[%02d] = %+3e, "
+                                    "%|80t|u3[%02d] = %+3e, "
+                                    "%|120t|u4[%02d] = %+3e");
+      for (size_t i = 0; i < frictionalNodes.size(); ++i)
+        if (frictionalNodes[i][0])
+          std::cout << boost::format(formatter) % i % u1[i] % i % u2[i] % i %
+                           u3[i] % i % u4[i] << std::endl;
+    }
   }
   catch (Dune::Exception &e) {
     Dune::derr << "Dune reported error: " << e << std::endl;
diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset
index d5f9d7a6ca4d7ddb25e0c5b766e4be0a897a1ea2..7f89be797905faac228fe4a1b2c1099e882213c1 100644
--- a/src/one-body-sample.parset
+++ b/src/one-body-sample.parset
@@ -2,6 +2,7 @@
 timesteps = 100
 
 verbose = false
+printFrictionalBoundary = false
 
 useNonlinearGS = false # true
 useGS = false