From eb088793e6864c6e7c5b65f63d5e45a0f28c8ca6 Mon Sep 17 00:00:00 2001
From: Ray Speth <speth@mit.edu>
Date: Sun, 26 Apr 2015 11:42:21 -0400
Subject: [PATCH] [Kinetics] Fix indices used in calculating coverage-dependent
 rates

The coverage dependent rate calculation is passed the array of coverages,
which include only the interface species, so the species index in the
interface phase is what is needed here.
---
 src/kinetics/InterfaceKinetics.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/kinetics/InterfaceKinetics.cpp b/src/kinetics/InterfaceKinetics.cpp
index 949dc5a15..bdfaa0db8 100644
--- a/src/kinetics/InterfaceKinetics.cpp
+++ b/src/kinetics/InterfaceKinetics.cpp
@@ -966,7 +966,7 @@ void InterfaceKinetics::addReaction(shared_ptr<Reaction> r_base)
     for (map<string, CoverageDependency>::const_iterator iter = r.coverage_deps.begin();
          iter != r.coverage_deps.end();
          ++iter) {
-        size_t k = kineticsSpeciesIndex(iter->first);
+        size_t k = thermo(reactionPhaseIndex()).speciesIndex(iter->first);
         rate.addCoverageDependence(k, iter->second.a, iter->second.m, iter->second.E);
     }