From a375d6b2edd944aee2d0df73fdcf9f58614f658f Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Tue, 25 Oct 2016 18:33:59 +0100
Subject: [PATCH] solidRegionDiffNo: Include the boundary faces in the
 calculation of the diffusion number

Patch contributed by Mattijs Janssens
Resolves bug-report http://bugs.openfoam.org/view.php?id=2303
---
 .../heatTransfer/chtMultiRegionFoam/solid/solidRegionDiffNo.C   | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/solidRegionDiffNo.C b/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/solidRegionDiffNo.C
index 57a024f430..de5307c11b 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/solidRegionDiffNo.C
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/solidRegionDiffNo.C
@@ -46,7 +46,7 @@ Foam::scalar Foam::solidRegionDiffNo
       / fvc::interpolate(Cprho)
     );
 
-    DiNum = gMax(kapparhoCpbyDelta.primitiveField())*runTime.deltaT().value();
+    DiNum = max(kapparhoCpbyDelta).value()*runTime.deltaT().value();
 
     meanDiNum = (average(kapparhoCpbyDelta)).value()*runTime.deltaT().value();