From 8f99f732f09ba43765436921440ee0528241aa05 Mon Sep 17 00:00:00 2001
From: ignis <ignis@postech.ac.kr>
Date: Wed, 23 Aug 2017 21:32:38 +0900
Subject: [PATCH] remove selective species update function

---
 solvers_post/SLFMFoam/SLFMFoam.C       |  8 +-------
 solvers_post/SLFMFoam/createFields.H   |  5 -----
 solvers_post/SLFMFoam/preIntegration.H |  9 ++++-----
 solvers_post/SLFMFoam/updateYi.H       | 10 ++++------
 4 files changed, 9 insertions(+), 23 deletions(-)

diff --git a/solvers_post/SLFMFoam/SLFMFoam.C b/solvers_post/SLFMFoam/SLFMFoam.C
index 2a293f2..ec55e8e 100644
--- a/solvers_post/SLFMFoam/SLFMFoam.C
+++ b/solvers_post/SLFMFoam/SLFMFoam.C
@@ -95,16 +95,10 @@ int main(int argc, char *argv[])
 
         turbulence->correct();
 
-            #include "updateYi.H" //update species
         if(runTime.write() == true)
         {
-            rho.write();
             #include "updateYi.H" //update species
-            forAll(postSpecieIndices, yi)
-            {
-                const label y = postSpecieIndices[yi];
-                Y[y].write();
-            }
+            rho.write();
         }
 
         Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
diff --git a/solvers_post/SLFMFoam/createFields.H b/solvers_post/SLFMFoam/createFields.H
index 8443f0a..4ef7655 100644
--- a/solvers_post/SLFMFoam/createFields.H
+++ b/solvers_post/SLFMFoam/createFields.H
@@ -11,10 +11,6 @@ thermo.validate(args.executable(), "h", "e");
 
 basicSpecieMixture& composition = thermo.composition();
 PtrList<volScalarField>& Y = composition.Y();
-forAll (Y, i)
-{
-    Y[i].writeOpt() = IOobject::NO_WRITE;
-}
 
 word inertSpecie(thermo.lookup("inertSpecie"));
 
@@ -303,7 +299,6 @@ forAll(postSpecieNames, i)
     const word &spName (postSpecieNames[i]);
     const label spI (composition.species()[spName]);
     postSpecieIndices[i] = spI;
-    Y[spI].writeOpt() = IOobject::AUTO_WRITE;
 }
 scalar stoiMF(readScalar(SLFMdict.lookup("stoiMF")));
 scalar stretchFac(readScalar(SLFMdict.lookup("stretchFactor")));
diff --git a/solvers_post/SLFMFoam/preIntegration.H b/solvers_post/SLFMFoam/preIntegration.H
index a5c4b05..2c3f295 100644
--- a/solvers_post/SLFMFoam/preIntegration.H
+++ b/solvers_post/SLFMFoam/preIntegration.H
@@ -26,9 +26,9 @@ scalarFieldArray2d Rtable
 
 scalarFieldArray3d Ytable
 (
-    postSpecieNames.size(),
+    Y.size(),
     scalarFieldArray2d
-        (
+    (
         etamax,
         scalarFieldArray1d
         (
@@ -58,10 +58,9 @@ for(label j=0 ; j<etamax ; j++)
         {
             Ttable[j][v][n] = bpdf.betaIntegrate(T_SLFM[n]);
             Rtable[j][v][n] = bpdf.betaIntegrate(Rgas_SLFM[n]);
-            forAll(postSpecieIndices, yi)
+            forAll(Y, yi)
             {
-                const label y = postSpecieIndices[yi];
-                Ytable[yi][j][v][n] = bpdf.betaIntegrate(Y_SLFM[y][n]);
+                Ytable[yi][j][v][n] = bpdf.betaIntegrate(Y_SLFM[yi][n]);
             }
         }
     }
diff --git a/solvers_post/SLFMFoam/updateYi.H b/solvers_post/SLFMFoam/updateYi.H
index c653619..09a38c7 100644
--- a/solvers_post/SLFMFoam/updateYi.H
+++ b/solvers_post/SLFMFoam/updateYi.H
@@ -31,9 +31,8 @@ forAll(rho, cellI)
     scalar jfac = jfc[cellI], vfac = vfc[cellI], nfac = nfc[cellI];
 
 
-    forAll(postSpecieIndices, yi)
+    forAll(Y, yi)
     {
-        const label y = postSpecieIndices[yi];
         scalar y00 = Ytable[yi][jl][vl][nl]*(1-jfac)+Ytable[yi][jh][vl][nl]*jfac;
         scalar y01 = Ytable[yi][jl][vl][nh]*(1-jfac)+Ytable[yi][jh][vl][nh]*jfac;
         scalar y10 = Ytable[yi][jl][vh][nl]*(1-jfac)+Ytable[yi][jh][vh][nl]*jfac;
@@ -44,15 +43,14 @@ forAll(rho, cellI)
 
         scalar ygas = y0*(1-nfac)+y1*nfac;
 
-        Y[y][cellI] = ygas;
+        Y[yi][cellI] = ygas;
     }
 }
 
 //correct processor boundary value of T and rho
-forAll(postSpecieIndices, yi)
+forAll(Y, yi)
 {
-    const label y = postSpecieIndices[yi];
-    Y[y].correctBoundaryConditions();
+    Y[yi].correctBoundaryConditions();
 }