diff --git a/dune/tectonic/nicefunction.hh b/dune/tectonic/nicefunction.hh
index cc779fc57b11c213897ed38a9f1f045708bd2cb7..fe2008d9c0c096e5af79b7366ed78fab3dcde6cf 100644
--- a/dune/tectonic/nicefunction.hh
+++ b/dune/tectonic/nicefunction.hh
@@ -35,11 +35,11 @@ class RuinaFunction : public NiceFunction {
       : coefficient(coefficient),
         a(a),
         mu(mu),
-        rho(exp(-mu / a)),
         eta(eta),
         normalStress(normalStress),
         compound_state(b * (state - log(eta * L))),
-        h(h) {}
+        h(h),
+        rho(exp(-(mu + compound_state) / a)) {}
 
   /*
     If mu and sigma_n denote the coefficient of friction and the