From 6d265b0a5f2a5062e69f541cbae6957d7a1423ae Mon Sep 17 00:00:00 2001
From: Graham <graham@noisy.opencfd.co.uk>
Date: Wed, 18 Mar 2009 10:29:05 +0000
Subject: [PATCH] Modifying dsmcFields to allow negative rhoNMean values, only
 cheching for min(mag(rhoNMean)) < VSMALL

---
 .../functionObjects/utilities/dsmcFields/dsmcFields.C      | 7 +++----
 1 file changed, 3 insertions(+), 4 deletions(-)

diff --git a/src/postProcessing/functionObjects/utilities/dsmcFields/dsmcFields.C b/src/postProcessing/functionObjects/utilities/dsmcFields/dsmcFields.C
index 3aaa9c514a4..273c88f3699 100644
--- a/src/postProcessing/functionObjects/utilities/dsmcFields/dsmcFields.C
+++ b/src/postProcessing/functionObjects/utilities/dsmcFields/dsmcFields.C
@@ -137,7 +137,7 @@ void Foam::dsmcFields::write()
             iDofMeanName
         );
 
-        if (min(rhoNMean).value() > VSMALL)
+        if (min(mag(rhoNMean)).value() > VSMALL)
         {
             Info<< "Calculating dsmcFields." << endl;
 
@@ -223,10 +223,9 @@ void Foam::dsmcFields::write()
         }
         else
         {
-            Info<< "Small or negative value (" << min(rhoNMean)
+            Info<< "Small value (" << min(mag(rhoNMean))
                 << ") found in rhoNMean field. "
-                << "Not calculating dsmcFields to avoid division by zero "
-                << "or invalid results."
+                << "Not calculating dsmcFields to avoid division by zero."
                 << endl;
         }
     }
-- 
GitLab