diff --git a/dune/gfe/nonplanarcosseratshellenergy.hh b/dune/gfe/nonplanarcosseratshellenergy.hh
index cfa82d662c2cfe22508d198cad580e304f978295..4901f8d5f6d1aaf08437c28ce1ca716d36496fac 100644
--- a/dune/gfe/nonplanarcosseratshellenergy.hh
+++ b/dune/gfe/nonplanarcosseratshellenergy.hh
@@ -135,6 +135,7 @@ public:
     return mu_ * Dune::GFE::sym(S).frobenius_norm2() + mu_c_ * Dune::GFE::skew(S).frobenius_norm2() + lambda_ * 0.5 * Dune::GFE::traceSquared(S);
   }
 
+  //For b1 = b2 = 1 and b3 = 1/3, this reduces to S.frobenius_norm2()
   RT W_curv(const Dune::FieldMatrix<field_type,3,3>& S) const
   {
     return mu_ * L_c_ * L_c_ * (b1_ * Dune::GFE::dev(Dune::GFE::sym(S)).frobenius_norm2()