From b8a5d2ef595571ba7c0f8d63b7dba2f01c60fe18 Mon Sep 17 00:00:00 2001
From: andy <a.heather@opencfd.co.uk>
Date: Fri, 10 Oct 2008 14:55:46 +0100
Subject: [PATCH] minor cosmetics change

---
 .../parcels/Templates/ReactingParcel/ReactingParcel.C        | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
index 5fe3aa16333..d7b2f2b9695 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
@@ -57,6 +57,7 @@ void Foam::ReactingParcel<ParcelType>::calcCoupled
     // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
     const vector U0 = this->U_;
     const scalar mass0 = this->mass();
+    const scalar cp0 = this->cp_;
     const scalar np0 = this->nParticle_;
     const scalar T0 = this->T_;
 
@@ -116,11 +117,11 @@ void Foam::ReactingParcel<ParcelType>::calcCoupled
     T1 -=
         td.constProps().Ldevol()
        *sum(dMassMT)
-       /(0.5*(mass0 + mass1)*this->cp_);
+       /(0.5*(mass0 + mass1)*cp0);
 
     // Add retained enthalpy from surface reaction to particle and remove
     // from gas
-    T1 += dhRet/(0.5*(mass0 + mass1)*this->cp_);
+    T1 += dhRet/(0.5*(mass0 + mass1)*cp0);
     dhTrans -= dhRet;
 
     // Correct dhTrans to account for enthalpy of evolved volatiles
-- 
GitLab