diff --git a/dumux/porousmediumflow/nonequilibrium/thermal/localresidual.hh b/dumux/porousmediumflow/nonequilibrium/thermal/localresidual.hh
index 4fda272afa66ce6af8676ca60d7e2a7c8ae73efb..8b28d2983402048ef67e15f5b48826180441b564 100644
--- a/dumux/porousmediumflow/nonequilibrium/thermal/localresidual.hh
+++ b/dumux/porousmediumflow/nonequilibrium/thermal/localresidual.hh
@@ -79,8 +79,8 @@ public:
         // add up parts on the one energy equation
         storage[energyEq0Idx] += volVars.porosity()
                                  * volVars.density(phaseIdx)
-                                 * volVars.saturation(phaseIdx)
-                                 * (volVars.internalEnergy(phaseIdx) - gravityPotential);
+                                 * volVars.internalEnergy(phaseIdx)
+                                 * volVars.saturation(phaseIdx);
 
     }