diff --git a/include/utils/parameters.h b/include/utils/parameters.h
index db2722ab02..f05628be14 100644
--- a/include/utils/parameters.h
+++ b/include/utils/parameters.h
@@ -37,7 +37,8 @@
 namespace libMesh
 {
 /**
- * Helper functions for printing scalar, vector and vector<vector> types.  Called from Parameters::Parameter<T>::print(...).
+ * Helper functions for printing scalar, vector, vector<vector> and vector<vector<vector>>
+ * types.  Called from Parameters::Parameter<T>::print(...).
  */
 template<typename P>
 void print_helper(std::ostream & os, const P * param);
@@ -48,6 +49,9 @@ void print_helper(std::ostream & os, const std::vector<P> * param);
 template<typename P>
 void print_helper(std::ostream & os, const std::vector<std::vector<P>> * param);
 
+template<typename P>
+void print_helper(std::ostream & os, const std::vector<std::vector<std::vector<P>>> * param);
+
 template<typename P1, typename P2, typename C, typename A>
 void print_helper(std::ostream & os, const std::map<P1, P2, C, A> * param);
 
@@ -566,6 +570,16 @@ void print_helper(std::ostream & os, const std::vector<std::vector<P>> * param)
       os << p << " ";
 }
 
+//non-member vector<vector<vector>> print function
+template<typename P>
+void print_helper(std::ostream & os, const std::vector<std::vector<std::vector<P>>> * param)
+{
+  for (const auto & pvv : *param)
+    for (const auto & pv : pvv)
+      for (const auto & p : pv)
+        os << p << " ";
+}
+
 //non-member map print function
 template<typename P1, typename P2, typename C, typename A>
 void print_helper(std::ostream & os, const std::map<P1, P2, C, A> * param)