From 9a67ae84c4c13aaf6b5481423c687d2501e5364b Mon Sep 17 00:00:00 2001 From: Dave Goodwin Date: Thu, 15 Dec 2005 17:07:06 +0000 Subject: [PATCH] added statements to getEquilibriumConstants to force a T update on the next call. Array m_rkcn was being left containing incorrect data. Problem is now fixed. --- Cantera/src/GasKinetics.cpp | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/Cantera/src/GasKinetics.cpp b/Cantera/src/GasKinetics.cpp index e1ab80814..70dec73d1 100755 --- a/Cantera/src/GasKinetics.cpp +++ b/Cantera/src/GasKinetics.cpp @@ -139,6 +139,10 @@ namespace Cantera { for (i = 0; i < m_ii; i++) { kc[i] = exp(-rkc[i]*rrt + m_dn[i]*logStandConc); } + + // force an update of T-dependent properties, so that m_rkcn will + // be updated before it is used next. + m_kdata->m_temp = 0.0; } /** @@ -440,7 +444,7 @@ namespace Cantera { */ const vector_fp& m_rkc = m_kdata->m_rkcn; for (int i = 0; i < m_ii; i++) { - krev[i] *= m_rkc[i]; + krev[i] *= m_rkc[i]; } } }