diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index 357c2457ac235cc7a7fda5d7a7132e4dd75cb5cf..9f0d6f4974093a9620e811645a06d0c2bd33156c 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -196,7 +196,10 @@ int main (int argc, char *argv[]) try
     for (int i=0; i<dirichletNodes.size(); i++)
         dirichletNodes[i] = dirichletBoundary.containsVertex(i);
 
-    sampleOnBitField(*complex.continua_["continuum"].grid_, dirichletValues[0], dirichletValues.back(), dirichletNodes);
+    sampleOnBitField(*complex.continua_["continuum"].grid_, 
+                     dirichletValues[0], 
+                     complex.continua_["continuum"].dirichletValues_, 
+                     dirichletNodes);
     
     /////////////////////////////////////////////////////////////////////
     //  Create the two interface boundary patches
@@ -459,7 +462,6 @@ int main (int argc, char *argv[]) try
                                                                                                   &rodLocalStiffness,
                                                                                                   &rodSolver,
                                                                                                   &stiffnessMatrix3d,
-                                                                                                  &dirichletValues.back(),
                                                                                                   solver,
                                                                                                   &localAssembler);