From 523f2eb35a16d76f39d330e37d95e673d856b422 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Mon, 23 Apr 2012 15:44:42 +0100
Subject: [PATCH] Gauss grading bug fix, Zeljko Tukovic

---
 .../finiteArea/gradSchemes/gaussFaGrad/gaussFaGrad.C        | 6 ++----
 1 file changed, 2 insertions(+), 4 deletions(-)

diff --git a/src/finiteArea/finiteArea/gradSchemes/gaussFaGrad/gaussFaGrad.C b/src/finiteArea/finiteArea/gradSchemes/gaussFaGrad/gaussFaGrad.C
index b78686eeb..a336770e5 100644
--- a/src/finiteArea/finiteArea/gradSchemes/gaussFaGrad/gaussFaGrad.C
+++ b/src/finiteArea/finiteArea/gradSchemes/gaussFaGrad/gaussFaGrad.C
@@ -22,8 +22,6 @@ License
     along with OpenFOAM; if not, write to the Free Software Foundation,
     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
 
-Description
-    
 \*---------------------------------------------------------------------------*/
 
 #include "gaussFaGrad.H"
@@ -69,12 +67,12 @@ gaussGrad<Type>::grad
 
     GeometricField<GradType, faPatchField, areaMesh>& gGrad = tgGrad();
 
+    gGrad -= vsf*fac::edgeIntegrate(vsf.mesh().Le());
+
     // Remove component of gradient normal to surface (area)
     const areaVectorField& n = vsf.mesh().faceAreaNormals();
-//     const areaScalarField& K = vsf.mesh().faceCurvatures();
 
     gGrad -= n*(n & gGrad);
-//     gGrad -= K*(n*vsf);
     gGrad.correctBoundaryConditions();
 
     gGrad.rename("grad(" + vsf.name() + ')');