diff --git a/Cantera/src/equil/vcs_solve_TP.cpp b/Cantera/src/equil/vcs_solve_TP.cpp index 36d727eb9..5c768bf6f 100644 --- a/Cantera/src/equil/vcs_solve_TP.cpp +++ b/Cantera/src/equil/vcs_solve_TP.cpp @@ -4411,8 +4411,8 @@ namespace VCSnonideal { if (tPhMoles_ptr[iph] > 0.0) { feSpecies[kspec] = m_SSfeSpecies[kspec] + log(actCoeff_ptr[kspec] * VCS_DELETE_MINORSPECIES_CUTOFF) - - tlogMoles[m_phaseID[kspec]] - m_lnMnaughtSpecies[kspec]; - + m_chargeSpecies[kspec] * m_Faraday_dim * m_phasePhi[iphase]; ; + - tlogMoles[m_phaseID[kspec]] - m_lnMnaughtSpecies[kspec] + + m_chargeSpecies[kspec] * m_Faraday_dim * m_phasePhi[iphase]; } else { feSpecies[kspec] = m_SSfeSpecies[kspec] - m_lnMnaughtSpecies[kspec] + m_chargeSpecies[kspec] * m_Faraday_dim * m_phasePhi[iphase];