void TurbulentKinetics::update_rates_T() { doublereal T = thermo().temperature(); doublereal P = thermo().pressure(); m_logStandConc = log(thermo().standardConcentration()); doublereal logT = log(T); doublereal TempFluc = Tprime(); if (T != m_temp) { if (!m_rfn.empty()) { m_rates.updateTurb(T, logT, &m_rfn[0],TempFluc); } if (!m_rfn_low.empty()) { m_falloff_low_rates.updateTurb(T, logT, &m_rfn_low[0],TempFluc); m_falloff_high_rates.updateTurb(T, logT, &m_rfn_high[0],TempFluc); } if (!falloff_work.empty()) { m_falloffn.updateTemp(T, &falloff_work[0]); } updateKc(); m_ROP_ok = false; } if (T != m_temp || P != m_pres) { if (m_plog_rates.nReactions()) { m_plog_rates.updateTurb(T, logT, &m_rfn[0],TempFluc); m_ROP_ok = false; } if (m_cheb_rates.nReactions()) { m_cheb_rates.updateTurb(T, logT, &m_rfn[0],TempFluc); m_ROP_ok = false; } } m_pres = P; m_temp = T; }
void T(void) { affiche_balise_ouvrante(__func__, TRACE_XML); F(); Tprime(); affiche_balise_fermante(__func__, TRACE_XML); }