From 94e86ddd977bc58e05b720e89f839d5c70239105 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Mon, 5 Sep 2011 19:12:42 +0100
Subject: [PATCH] MULES: correct the handling of boundary conditions for
 limiting the sum

---
 .../fvMatrices/solvers/MULES/MULESTemplates.C | 21 +++++++++++--------
 1 file changed, 12 insertions(+), 9 deletions(-)

diff --git a/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C b/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
index eaa64aebf5..55478e59aa 100644
--- a/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
+++ b/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
@@ -670,17 +670,20 @@ void Foam::MULES::limitSum(SurfaceScalarFieldList& phiPsiCorrs)
 
     forAll(phiPsiCorrs[0].boundaryField(), patchi)
     {
-        UPtrList<scalarField> phiPsiCorrsPatch(phiPsiCorrs.size());
-        forAll(phiPsiCorrs, phasei)
+        if (phiPsiCorrs[0].boundaryField()[patchi].coupled())
         {
-            phiPsiCorrsPatch.set
-            (
-                phasei,
-                &phiPsiCorrs[phasei].boundaryField()[patchi]
-            );
-        }
+            UPtrList<scalarField> phiPsiCorrsPatch(phiPsiCorrs.size());
+            forAll(phiPsiCorrs, phasei)
+            {
+                phiPsiCorrsPatch.set
+                (
+                    phasei,
+                    &phiPsiCorrs[phasei].boundaryField()[patchi]
+                );
+            }
 
-        limitSum(phiPsiCorrsPatch);
+            limitSum(phiPsiCorrsPatch);
+        }
     }
 }