From 98b9eaac986b118a42f268bc664993afddcf5fb4 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Mon, 19 May 2014 17:19:36 +0100
Subject: [PATCH] BUG: parcel phase change - bypass calc if phase fraction <
 SMALL - mantis #1298

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

diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
index 32427cba8a..8ae6f5fe09 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
@@ -60,14 +60,14 @@ void Foam::ReactingParcel<ParcelType>::calcPhaseChange
     typedef typename TrackData::cloudType::reactingCloudType reactingCloudType;
     PhaseChangeModel<reactingCloudType>& phaseChange = td.cloud().phaseChange();
 
-    if (!phaseChange.active())
+    if (!phaseChange.active() || (YPhase < SMALL))
     {
         return;
     }
 
     scalar Tvap = phaseChange.Tvap(YComponents);
 
-    if (T < Tvap || YPhase < SMALL)
+    if (T < Tvap)
     {
         return;
     }