Calculate T properly with HP inputs in two-phase; closes #762

CoolProp · Aug 23, 2015 · 44ab82b · 44ab82b
1 parent 5c2a518
commit 44ab82b
Showing 1 changed file with 7 additions and 1 deletion.
diff --git a/src/Backends/Tabular/TabularBackends.cpp b/src/Backends/Tabular/TabularBackends.cpp
@@ -364,6 +364,7 @@ CoolPropDbl CoolProp::TabularBackend::calc_p(void){
 }
 CoolPropDbl CoolProp::TabularBackend::calc_T(void){
     PhaseEnvelopeData & phase_envelope = dataset->phase_envelope;
+    PureFluidSaturationTableData &pure_saturation = dataset->pure_saturation;
     if (using_single_phase_table){
         switch (selected_table){
         case SELECTED_PH_TABLE: return evaluate_single_phase_phmolar(iT, cached_single_phase_i, cached_single_phase_j);
@@ -377,7 +378,12 @@ CoolPropDbl CoolProp::TabularBackend::calc_T(void){
             return phase_envelope_sat(phase_envelope, iT, iP, _p);
         }
         else{
-            return _T;
+            if (ValidNumber(_T)){
+                return _T;
+            }
+            else{
+                return pure_saturation.evaluate(iT, _p, _Q, cached_saturation_iL, cached_saturation_iV);
+            }
         }
     }
 }