From 3d1fc441c2cfdec3025232e5b5c6f25d0465b3aa Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Thu, 24 Mar 2011 15:08:10 +0000
Subject: [PATCH] ENH: moved rho update

---
 applications/solvers/lagrangian/reactingParcelFoam/pEqn.H       | 2 --
 .../solvers/lagrangian/reactingParcelFoam/reactingParcelFoam.C  | 2 ++
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/applications/solvers/lagrangian/reactingParcelFoam/pEqn.H b/applications/solvers/lagrangian/reactingParcelFoam/pEqn.H
index 3c65ca9fb9..5599e861bb 100644
--- a/applications/solvers/lagrangian/reactingParcelFoam/pEqn.H
+++ b/applications/solvers/lagrangian/reactingParcelFoam/pEqn.H
@@ -88,8 +88,6 @@ else
 #include "rhoEqn.H"
 #include "compressibleContinuityErrs.H"
 
-rho = thermo.rho();
-
 U -= rAU*fvc::grad(p);
 U.correctBoundaryConditions();
 
diff --git a/applications/solvers/lagrangian/reactingParcelFoam/reactingParcelFoam.C b/applications/solvers/lagrangian/reactingParcelFoam/reactingParcelFoam.C
index b0a82278f0..fb89974a62 100644
--- a/applications/solvers/lagrangian/reactingParcelFoam/reactingParcelFoam.C
+++ b/applications/solvers/lagrangian/reactingParcelFoam/reactingParcelFoam.C
@@ -98,6 +98,8 @@ int main(int argc, char *argv[])
 
             turbulence->correct();
 
+            rho = thermo.rho();
+
             if (finalIter)
             {
                 mesh.data::remove("finalIteration");