diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 895da3fdc9c1786c5723c2652c9aaeff7eb71537..a2b5b42bf8aabb46d96363319fa36da019e9d8f8 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -131,6 +131,7 @@ int main(int argc, char *argv[]) {
     auto const solver_tolerance = parset.get<double>("solver.tolerance");
 
     auto const refinements = parset.get<size_t>("grid.refinements");
+    double const refinement_factor = std::pow(2, refinements);
 
     auto const verbose = parset.get<bool>("verbose");
     Solver::VerbosityMode const verbosity =
@@ -235,9 +236,11 @@ int main(int argc, char *argv[]) {
         stiffnessMatrix.mmv(u4, b4);
         // Apply Dirichlet condition
         for (int i = 0; i < finestSize; ++i)
-          if (ignoreNodes[i].count() == dim)
+          if (ignoreNodes[i].count() == dim) {
             functions.get("dirichletCondition")
                 .evaluate(h * run, u4_diff[i][0]);
+            u4_diff[i][0] /= refinement_factor;
+          }
 
         for (int state_fpi = 0;
              state_fpi < parset.get<int>("solver.tnnmg.fixed_point_iterations");
diff --git a/src/one-body-sample.py b/src/one-body-sample.py
index 803f86f6cc0077ebcb30c0faab36271a4dd58668..81d09658a8eebbe6a01479ab4805917c7627e848 100644
--- a/src/one-body-sample.py
+++ b/src/one-body-sample.py
@@ -14,9 +14,9 @@ class neumannCondition:
 class dirichletCondition:
     def __call__(self, x):
         # return 0
-        fst = 3e-5
-        snd = 1e-5
-        trd = 3e-5
+        fst = 5e-4
+        snd = 2e-4
+        trd = 5e-4
         if x < 1.0/3:
             return fst
         elif x < 2.0/3: