From ac2a275a7518a96ab02004439fa015e40677e275 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Wed, 2 Dec 2015 11:47:11 +0000
Subject: [PATCH] wallHeatFlux: Corrected sign of radiative heat-flux Corrects
 feature-request http://www.openfoam.org/mantisbt/view.php?id=1856 Patch
 provided by Juho Peltola

---
 .../utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C b/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C
index 76550b9699..d1c37e70cc 100644
--- a/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C
+++ b/applications/utilities/postProcessing/wall/wallHeatFlux/wallHeatFlux.C
@@ -82,7 +82,7 @@ int main(int argc, char *argv[])
             if (isA<wallFvPatch>(mesh.boundary()[patchi]))
             {
                 scalar convFlux = gSum(magSf[patchi]*patchHeatFlux[patchi]);
-                scalar radFlux = gSum(magSf[patchi]*patchRadHeatFlux[patchi]);
+                scalar radFlux = -gSum(magSf[patchi]*patchRadHeatFlux[patchi]);
 
                 Info<< mesh.boundary()[patchi].name() << endl
                     << "    convective: " << convFlux << endl
@@ -135,7 +135,7 @@ int main(int argc, char *argv[])
             forAll(totalWallHeatFlux.boundaryField(), patchi)
             {
                 totalWallHeatFlux.boundaryField()[patchi] =
-                    patchHeatFlux[patchi] + patchRadHeatFlux[patchi];
+                    patchHeatFlux[patchi] - patchRadHeatFlux[patchi];
             }
 
             totalWallHeatFlux.write();