diff --git a/packages/muelu/src/Operators/MueLu_MultiPhys_def.hpp b/packages/muelu/src/Operators/MueLu_MultiPhys_def.hpp
index 6155dd7f092f..5a04f8ebe20f 100644
--- a/packages/muelu/src/Operators/MueLu_MultiPhys_def.hpp
+++ b/packages/muelu/src/Operators/MueLu_MultiPhys_def.hpp
@@ -142,7 +142,9 @@ void MultiPhys<Scalar, LocalOrdinal, GlobalOrdinal, Node>::compute(bool reuse) {
   // Generate the (iii,iii) Hierarchy
 
   for (int iii = 0; iii < nBlks_; iii++) {
-    arrayOfParamLists_[iii]->sublist("user data").set("Coordinates", arrayOfCoords_[iii]);
+    if (arrayOfCoords_ != Teuchos::null) {
+      arrayOfParamLists_[iii]->sublist("user data").set("Coordinates", arrayOfCoords_[iii]);
+    }
 
     bool wantToRepartition = false;
     if (paramListMultiphysics_->isParameter("repartition: enable"))
@@ -312,6 +314,8 @@ void MultiPhys<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
   }
   oss << std::endl;
 
+  out << oss.str();
+
   for (int i = 0; i < nBlks_; i++) {
     arrayOfHierarchies_[i]->describe(out, GetVerbLevel());
   }