diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index c0d9dcd8c4bfd3aea87a6c5306e6b7c0c1bb4105..0137963c3ad7fa1fdf33dfb7521b27da59e65520 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -471,6 +471,10 @@ int main (int argc, char *argv[]) try
             
             rodContinuumSteklovPoincareStep.iterate(lambda);
             
+            // get the subdomain solutions
+            rodX = rodContinuumSteklovPoincareStep.rodSubdomainSolutions_["rod"];
+            x3d  = rodContinuumSteklovPoincareStep.continuumSubdomainSolutions_["continuum"];
+            
         } else
             DUNE_THROW(NotImplemented, ddType << " is not a known domain decomposition algorithm");