diff --git a/src/transport/MultiTransport.cpp b/src/transport/MultiTransport.cpp
index 88b41d243..b09cb9e00 100644
--- a/src/transport/MultiTransport.cpp
+++ b/src/transport/MultiTransport.cpp
@@ -542,6 +542,11 @@ void MultiTransport::updateThermal_T()
      *       The original Dixon-Lewis paper subtracted 1.5 here.
      */
     const vector_fp& cp = ((IdealGasPhase*)m_thermo)->cp_R_ref();
+
+    vector_fp cp_R_all(m_thermo->nSpecies());
+    m_thermo->getCp_R_ref(&cp_R_all[0]);
+    const vector_fp& cp = cp_R_all;
+    
     for (size_t k = 0; k < m_nsp; k++) {
         m_cinternal[k] = cp[k] - 2.5;
     }