From e7d92199f350424aa46a33658a04b15ff6fee1dd Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Tue, 4 Dec 2012 17:05:29 +0000
Subject: [PATCH] ENH: Pyrolysis - updated handling for phiGas

---
 .../pyrolysisModels/reactingOneDim/reactingOneDim.C          | 5 +----
 1 file changed, 1 insertion(+), 4 deletions(-)

diff --git a/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C b/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C
index 3588060c72..84fb91407f 100644
--- a/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C
+++ b/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C
@@ -251,15 +251,12 @@ void reactingOneDim::solveEnergy()
 
     tmp<volScalarField> alpha(solidThermo_.alpha());
 
-    const surfaceScalarField phiGas(fvc::interpolate(phiHsGas_));
-
     fvScalarMatrix hEqn
     (
         fvm::ddt(rho_, h_)
       - fvm::laplacian(alpha, h_)
      ==
         chemistrySh_
-      + fvc::div(phiGas)
     );
 
     if (regionMesh().moving())
@@ -334,7 +331,7 @@ reactingOneDim::reactingOneDim(const word& modelType, const fvMesh& mesh)
             "phiGas",
             time().timeName(),
             regionMesh(),
-            IOobject::NO_READ,
+            IOobject::READ_IF_PRESENT,
             IOobject::AUTO_WRITE
         ),
         regionMesh(),