diff --git a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
index ee77cbec5967781b73f331c4470f20e1076c0ce4..f34fa905abf27d1c15bf09677a19ab4c053cd310 100644
--- a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
+++ b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
@@ -1105,7 +1105,7 @@ iterateWithContact(std::map<std::pair<std::string,std::string>, RigidBodyMotion<
     totalX3d = multigridStep->getSol();
 
     // Separate 3d solution vector
-    std::vector<VectorType> x3d;
+    std::vector<VectorType> x3d(continuumName.size());
     contactAssembler->postprocess(totalX3d, x3d);
     
     // the subdomain solutions in canonical coordinates, stored in a map