From b889cc58839cb477f3e07ffcd01b76c46c47645e Mon Sep 17 00:00:00 2001
From: Ian Bell <ian.h.bell@gmail.com>
Date: Sat, 22 Oct 2016 17:38:43 -0600
Subject: [PATCH] Now errors if T < Tmin and p > ptriple; closes #1165

---
 src/Backends/Helmholtz/HelmholtzEOSMixtureBackend.cpp | 9 +++++++++
 1 file changed, 9 insertions(+)

diff --git a/src/Backends/Helmholtz/HelmholtzEOSMixtureBackend.cpp b/src/Backends/Helmholtz/HelmholtzEOSMixtureBackend.cpp
index c3f24df0..71f4366c 100644
--- a/src/Backends/Helmholtz/HelmholtzEOSMixtureBackend.cpp
+++ b/src/Backends/Helmholtz/HelmholtzEOSMixtureBackend.cpp
@@ -1436,6 +1436,15 @@ void HelmholtzEOSMixtureBackend::p_phase_determination_pure_or_pseudopure(int ot
         {
             case iT:
             {
+                if (_T < Tmin()){
+                    if (get_config_bool(DONT_CHECK_PROPERTY_LIMITS)){
+                        _phase = iphase_liquid;
+                    }
+                    else{
+                        throw ValueError(format("For now, we don't support T [%g K] below Tmin(saturation) [%g K]", _T, Tmin()));
+                    }
+                }
+                
                 CoolPropDbl T_vap =  0.1 + static_cast<double>(_TVanc);
                 CoolPropDbl T_liq = -0.1 + static_cast<double>(_TLanc);