From adb7d5fb499857cd5431c1b813440026b2f2cdfd Mon Sep 17 00:00:00 2001
From: henry <Henry Weller h.weller@opencfd.co.uk>
Date: Thu, 19 Jun 2008 10:42:35 +0100
Subject: [PATCH] Bug fix:
 http://openfoam.cfd-online.com/cgi-bin/forum/show.cgi?126/8089

---
 .../Schaeffer/SchaefferFrictionalStress.C                   | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/frictionalStressModel/Schaeffer/SchaefferFrictionalStress.C b/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/frictionalStressModel/Schaeffer/SchaefferFrictionalStress.C
index b3ddd44ff2..1e0ec17021 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/frictionalStressModel/Schaeffer/SchaefferFrictionalStress.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/frictionalStressModel/Schaeffer/SchaefferFrictionalStress.C
@@ -123,13 +123,13 @@ Foam::tmp<Foam::volScalarField> Foam::SchaefferFrictionalStress::muf
         )
     );
 
-    volScalarField& muf_ = tmuf();
+    volScalarField& muff = tmuf();
 
     forAll (D, celli)
     {
         if (alpha[celli] > alphaMax.value()-5e-2)
         {
-            muf_[celli] =
+            muff[celli] =
                 0.5*pf[celli]*sin(phi.value())
                /(
                     sqrt(1.0/6.0*(sqr(D[celli].xx() - D[celli].yy())
@@ -141,7 +141,7 @@ Foam::tmp<Foam::volScalarField> Foam::SchaefferFrictionalStress::muf
         }
     }
 
-    return muf_;
+    return tmuf;
 }