diff --git a/solvers_post/LagrangianCMCFoam/CMCintegration.H b/solvers_post/LagrangianCMCFoam/CMCintegration.H index cc12438..d1e7201 100644 --- a/solvers_post/LagrangianCMCFoam/CMCintegration.H +++ b/solvers_post/LagrangianCMCFoam/CMCintegration.H @@ -5,7 +5,7 @@ { scalarField Ytemp(Ysize, 0.0); - for(label i=0 ; i Y[o2index][celli]; + Switch fuelIncrease = Ytemp[fuelindex] > Y[fuelindex][celli]; + if( ignoreOFIncrease && (oxygenIncrease || fuelIncrease)) { - Y[i][celli] = Ytemp[i]; + // Do not update Yi + } + else + { + for(label i = 0 ; i < Ysize ; i++) + { + Y[i][celli] = Ytemp[i]; + } } } }