diff --git a/dune/tectonic/nicefunction.hh b/dune/tectonic/nicefunction.hh
index fe2008d9c0c096e5af79b7366ed78fab3dcde6cf..83ae7def1102c9c677d2ee9ac849aa24c5a19615 100644
--- a/dune/tectonic/nicefunction.hh
+++ b/dune/tectonic/nicefunction.hh
@@ -58,7 +58,7 @@ class RuinaFunction : public NiceFunction {
   void virtual evaluate(double const &x, double &y) const {
     double const arg = std::max(eta * x / h, rho);
     double const r = arg * (a * (std::log(arg) - 1) + mu + compound_state);
-    y = 1 / eta * normalStress * r + a * rho + mu + compound_state;
+    y = 1 / eta * normalStress * (r + a * rho + mu + compound_state);
     y *= coefficient;
     y *= h;
   }