From a748510054dfc27f378f5178dc1c00caab9d20e9 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Thu, 22 Nov 2012 13:14:41 +0000
Subject: [PATCH] ENH: Updated DESModelRegions function object

---
 .../utilities/DESModelRegions/DESModelRegions.C    | 14 +++++++++-----
 1 file changed, 9 insertions(+), 5 deletions(-)

diff --git a/src/postProcessing/functionObjects/utilities/DESModelRegions/DESModelRegions.C b/src/postProcessing/functionObjects/utilities/DESModelRegions/DESModelRegions.C
index 94e375304a..6d5c2f7eaa 100644
--- a/src/postProcessing/functionObjects/utilities/DESModelRegions/DESModelRegions.C
+++ b/src/postProcessing/functionObjects/utilities/DESModelRegions/DESModelRegions.C
@@ -130,7 +130,7 @@ void Foam::DESModelRegions::write()
             Info<< type() << " output:" << nl;
         }
 
-        tmp<volScalarField> result;
+        tmp<volScalarField> tresult;
 
         label DESpresent = false;
         if (mesh.foundObject<icoModel>("turbulenceModel"))
@@ -142,7 +142,7 @@ void Foam::DESModelRegions::write()
             {
                 const icoDESModel& des =
                     dynamic_cast<const icoDESModel&>(model);
-                result = des.LESRegion();
+                tresult = des.LESRegion();
                 DESpresent = true;
             }
         }
@@ -155,15 +155,17 @@ void Foam::DESModelRegions::write()
             {
                 const cmpDESModel& des =
                     dynamic_cast<const cmpDESModel&>(model);
-                result = des.LESRegion();
+                tresult = des.LESRegion();
                 DESpresent = true;
             }
         }
 
         if (DESpresent)
         {
+            const volScalarField& result = tresult();
+
             scalar prc =
-                gSum(result().internalField()*mesh.V())/gSum(mesh.V())*100.0;
+                gSum(result.internalField()*mesh.V())/gSum(mesh.V())*100.0;
 
             if (Pstream::master())
             {
@@ -174,9 +176,11 @@ void Foam::DESModelRegions::write()
             if (log_)
             {
                 Info<< "    LES = " << prc << " % (volume)" << nl
-                    << "    RES = " << 100.0 - prc << " % (volume)" << nl
+                    << "    RAS = " << 100.0 - prc << " % (volume)" << nl
                     << endl;
             }
+
+            result.write();
         }
         else
         {