From 3612abb1172299b3abe1c828eb997dd187a4c5c2 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Thu, 26 May 2011 17:00:34 +0100
Subject: [PATCH] ENH: Updated parcel properties for coneNozzleInjector nozzle

---
 .../InjectionModel/ConeNozzleInjection/ConeNozzleInjection.C  | 4 +---
 1 file changed, 1 insertion(+), 3 deletions(-)

diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeNozzleInjection/ConeNozzleInjection.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeNozzleInjection/ConeNozzleInjection.C
index 28ca032f2e..1339545c23 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeNozzleInjection/ConeNozzleInjection.C
+++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeNozzleInjection/ConeNozzleInjection.C
@@ -349,9 +349,7 @@ void Foam::ConeNozzleInjection<CloudType>::setProperties
     scalar massFlowRate =
         this->massTotal()*volumeFlowRate_().value(t)/this->volumeTotal();
 
-    scalar rho = this->owner().constProps().rho0();
-//    scalar Umag = massFlowRate/(parcel.rho()*Cd_().value(t)*(Ao - Ai));
-    scalar Umag = massFlowRate/(rho*Cd_().value(t)*(Ao - Ai));
+    scalar Umag = massFlowRate/(parcel.rho()*Cd_().value(t)*(Ao - Ai));
     parcel.U() = Umag*dirVec;
 
     // set particle diameter