diff --git a/Cantera/src/thermo/HMWSoln.cpp b/Cantera/src/thermo/HMWSoln.cpp index 51f8f76ee..47afa9db9 100644 --- a/Cantera/src/thermo/HMWSoln.cpp +++ b/Cantera/src/thermo/HMWSoln.cpp @@ -1454,6 +1454,7 @@ namespace Cantera { int leng = m_kk; m_electrolyteSpeciesType.resize(m_kk, cEST_polarNeutral); m_speciesSize.resize(leng); + m_speciesCharge_Stoich.resize(leng, 0.0); m_Aionic.resize(leng, 0.0); m_expg0_RT.resize(leng, 0.0); diff --git a/Cantera/src/thermo/HMWSoln_input.cpp b/Cantera/src/thermo/HMWSoln_input.cpp index 53a70b797..e5af2774c 100644 --- a/Cantera/src/thermo/HMWSoln_input.cpp +++ b/Cantera/src/thermo/HMWSoln_input.cpp @@ -1001,6 +1001,16 @@ namespace Cantera { */ m_waterProps = new WaterProps(dynamic_cast(m_waterSS)); + /* + * Fill in parameters for the calculation of the + * stoichiometric Ionic Strength + * + * The default is that stoich charge is the same as the + * regular charge. + */ + for (k = 0; k < m_kk; k++) { + m_speciesCharge_Stoich[k] = m_speciesCharge[k]; + } /* * Go get all of the coefficients and factors in the @@ -1067,18 +1077,6 @@ namespace Cantera { } - - /* - * Fill in parameters for the calculation of the - * stoichiometric Ionic Strength - * - * The default is that stoich charge is the same as the - * regular charge. - */ - m_speciesCharge_Stoich.resize(m_kk, 0.0); - for (k = 0; k < m_kk; k++) { - m_speciesCharge_Stoich[k] = m_speciesCharge[k]; - } /* * First look at the species database. * -> Look for the subelement "stoichIsMods"