From 7443011217a470da512101bf26753da35d835fc5 Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Wed, 21 Aug 2019 11:47:10 -0700
Subject: [PATCH] BUG: Fixing solid specie index in omega function

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

diff --git a/src/thermophysicalModels/solidChemistryModel/pyrolysisChemistryModel/pyrolysisChemistryModel.C b/src/thermophysicalModels/solidChemistryModel/pyrolysisChemistryModel/pyrolysisChemistryModel.C
index 3b892483e5..b284e01682 100644
--- a/src/thermophysicalModels/solidChemistryModel/pyrolysisChemistryModel/pyrolysisChemistryModel.C
+++ b/src/thermophysicalModels/solidChemistryModel/pyrolysisChemistryModel/pyrolysisChemistryModel.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           |
+    \\  /    A nd           | Copyright (C) 2019 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
                             | Copyright (C) 2013-2017 OpenFOAM Foundation
@@ -264,7 +264,7 @@ Foam::pyrolysisChemistryModel<CompType, SolidThermo, GasThermo>::omega
     for (label s=0; s<Nl; s++)
     {
         label si = R.lhs()[s].index;
-        const scalar exp = R.lhs()[si].exponent;
+        const scalar exp = R.lhs()[s].exponent;
 
         kf *=
             pow(c1[si]/Ys0_[si][celli], exp)