diff --git a/src/coupled/problem/PhysiologicalProblems.hpp b/src/coupled/problem/PhysiologicalProblems.hpp
index 79ecb2038b0d94a5bf4c735dbabbd5c0f28c15ce..f01d585a5c018ba3dfc8cb48865ca03bcf9f1b7a 100644
--- a/src/coupled/problem/PhysiologicalProblems.hpp
+++ b/src/coupled/problem/PhysiologicalProblems.hpp
@@ -6,6 +6,7 @@
 
 class ContractingVentricleProblem : public CoupledProblem {
   using ElasticityProblem::verbose;
+  using ElasticityProblem::meshesName;
   // Robin Boundary
   //std::array<double, 8> robinParams{10.0, 0.0, 5.0, 0.0, 1500.0, 1e-4, 1.0, 0.0};
   std::array<double, 8> robinParams{0.2, 0.0, 5.0e-3, 0.0, 1500.0, 1e-4, 1.0, 0.0};
@@ -19,17 +20,21 @@ class ContractingVentricleProblem : public CoupledProblem {
 
   VentricleProblem ventricleProblem;
 public:
-  explicit ContractingVentricleProblem() :
-      ventricleProblem() {
+  explicit ContractingVentricleProblem(std::string const &defaultMesh = "") :
+      ventricleProblem(defaultMesh) {
     Config::Get("TractionK", k);
     Config::Get("TractionC", c);
 
+    ElasticityProblem::meshesName = ventricleProblem.GetMeshName();
+    ElphyProblem::meshesName = ventricleProblem.GetMeshName();
   }
 
   ContractingVentricleProblem(const string &matName, const vector<double> &matPar) :
       ventricleProblem(matName, matPar) {
     Config::Get("TractionK", k);
     Config::Get("TractionC", c);
+    ElasticityProblem::meshesName = ventricleProblem.GetMeshName();
+    ElphyProblem::meshesName = ventricleProblem.GetMeshName();
   }
 
   string Name() const override {