[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

@@ -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);
     }