From 7194bf96d4026c74372432b17c57c9f4f6d1bf8f Mon Sep 17 00:00:00 2001
From: Ray Speth <speth@mit.edu>
Date: Tue, 10 Jul 2012 16:23:32 +0000
Subject: [PATCH] Fixed calculation of some multicomponent transport properties

For some properties, the internal temperature-dependent properties
were being updated in the wrong order.
---
 src/transport/GasTransport.cpp   | 2 ++
 src/transport/MultiTransport.cpp | 1 -
 2 files changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/transport/GasTransport.cpp b/src/transport/GasTransport.cpp
index 213abf900..f56c06161 100644
--- a/src/transport/GasTransport.cpp
+++ b/src/transport/GasTransport.cpp
@@ -212,6 +212,7 @@ void GasTransport::updateViscosity_T()
 
 void GasTransport::updateSpeciesViscosities()
 {
+    update_T();
     if (m_mode == CK_Mode) {
         for (size_t k = 0; k < m_nsp; k++) {
             m_visc[k] = exp(dot4(m_polytempvec, m_visccoeffs[k]));
@@ -229,6 +230,7 @@ void GasTransport::updateSpeciesViscosities()
 
 void GasTransport::updateDiff_T()
 {
+    update_T();
     // evaluate binary diffusion coefficients at unit pressure
     size_t ic = 0;
     if (m_mode == CK_Mode) {
diff --git a/src/transport/MultiTransport.cpp b/src/transport/MultiTransport.cpp
index e819e723a..9b8cfd53c 100644
--- a/src/transport/MultiTransport.cpp
+++ b/src/transport/MultiTransport.cpp
@@ -655,7 +655,6 @@ void MultiTransport::updateThermal_T()
     }
     // we need species viscosities and binary diffusion coefficients
     updateSpeciesViscosities();
-    update_T();
     updateDiff_T();
 
     // evaluate polynomial fits for A*, B*, C*