diff --git a/src/tpx/Sub.cpp b/src/tpx/Sub.cpp index bfc646f95..b6c9b0663 100644 --- a/src/tpx/Sub.cpp +++ b/src/tpx/Sub.cpp @@ -299,7 +299,7 @@ void Substance::update_sat() double gf = hp() - T*sp(); if (i==0) { - Rho = pp*MolWt()/(8314.0*T); // trial value = ideal gas + Rho = pp*MolWt()/(GasConstant*T); // trial value = ideal gas } else { Rho = Rhv; }