diff --git a/src/transport/GasTransport.cpp b/src/transport/GasTransport.cpp
index e0fb03e96..213abf900 100644
--- a/src/transport/GasTransport.cpp
+++ b/src/transport/GasTransport.cpp
@@ -136,7 +136,12 @@ bool GasTransport::initGas(GasTransportParams& tr)
 }
 
 void GasTransport::update_T(void) {
-    m_temp = m_thermo->temperature();
+    double T = m_thermo->temperature();
+    if (T == m_temp) {
+        return;
+    }
+
+    m_temp = T;
     m_kbt = Boltzmann * m_temp;
     m_sqrt_kbt = sqrt(Boltzmann*m_temp);
     m_logt = log(m_temp);
diff --git a/src/transport/MultiTransport.cpp b/src/transport/MultiTransport.cpp
index d198a8e32..e819e723a 100644
--- a/src/transport/MultiTransport.cpp
+++ b/src/transport/MultiTransport.cpp
@@ -307,6 +307,7 @@ void MultiTransport::getSpeciesFluxes(size_t ndim, const doublereal* const grad_
                                       size_t ldf, doublereal* const fluxes)
 {
     // update the binary diffusion coefficients if necessary
+    update_T();
     updateDiff_T();
 
     // If any component of grad_T is non-zero, then get the
@@ -456,6 +457,7 @@ void MultiTransport::getMassFluxes(const doublereal* state1, const doublereal* s
     m_thermo->getMoleFractions(DATA_PTR(m_molefracs));
 
     // update the binary diffusion coefficients if necessary
+    update_T();
     updateDiff_T();
 
     // If there is a temperature gradient, then get the
@@ -576,6 +578,7 @@ void MultiTransport::getMultiDiffCoeffs(const size_t ld, doublereal* const d)
     update_C();
 
     // update the binary diffusion coefficients
+    update_T();
     updateDiff_T();
 
     // evaluate L0000 if the temperature or concentrations have
@@ -652,6 +655,7 @@ void MultiTransport::updateThermal_T()
     }
     // we need species viscosities and binary diffusion coefficients
     updateSpeciesViscosities();
+    update_T();
     updateDiff_T();
 
     // evaluate polynomial fits for A*, B*, C*