From 813cc36f8e25f1819c04a8f4731ae60a08cca7ea Mon Sep 17 00:00:00 2001 From: Ray Speth Date: Sat, 23 May 2015 22:33:21 -0400 Subject: [PATCH] Fix undefined variables in ElectrodeKinetics with DEBUG_MODE enabled --- src/kinetics/ElectrodeKinetics.cpp | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/kinetics/ElectrodeKinetics.cpp b/src/kinetics/ElectrodeKinetics.cpp index 609ea8144..1413c1b6c 100644 --- a/src/kinetics/ElectrodeKinetics.cpp +++ b/src/kinetics/ElectrodeKinetics.cpp @@ -547,7 +547,7 @@ double ElectrodeKinetics::calcForwardROP_BV(size_t irxn, size_t iBeta, double io if (ro_fwd == 0) { throw CanteraError("ElectrodeKinetics::calcForwardROP_BV()", "forward orders pointer is zero ?!?"); } - double tmp = exp(- m_beta[iBeta] * m_deltaG0[irxn] * rrt); + double tmp = exp(- m_beta[iBeta] * m_deltaG0[irxn] / rt); double tmp2 = 1.0; const std::vector& kinSpeciesIDs = ro_fwd->kinSpeciesIDs_; const std::vector& kinSpeciesOrders = ro_fwd->kinSpeciesOrders_; @@ -562,12 +562,13 @@ double ElectrodeKinetics::calcForwardROP_BV(size_t irxn, size_t iBeta, double io // // Calculate the chemical reaction rate constant // + double iorc = m_rfn[irxn] * m_perturb[irxn]; double kf = iorc * tmp; // // Calculate the electrochemical factor // double eamod = m_beta[iBeta] * deltaElectricEnergy_[irxn]; - kf *= exp(- eamod * rrt); + kf *= exp(- eamod / rt); // // Calculate the forward rate of progress // -> get the pointer for the orders