fix Tmin higher than 0

EcoExtreML · Jun 30, 2022 · a50de3e · a50de3e
1 parent 5dc0a82
commit a50de3e
Showing 1 changed file with 7 additions and 2 deletions.
diff --git a/src/HT_frez.m b/src/HT_frez.m
@@ -8,8 +8,13 @@
                     else
                         TT_CRIT(MN)=T0+g*T0*HHH(MN)/L_f/1e4/2;
                     end
-                    if TT_CRIT(MN)<=T0+Tmin
-                        TT_CRIT(MN)=T0+Tmin;  % unit K
+                    if Tmin>0
+                        TTmin=0;
+                    else
+                        TTmin=Tmin;
+                    end
+                    if TT_CRIT(MN)<=T0+TTmin
+                        TT_CRIT(MN)=T0+TTmin;  % unit K
                     elseif TT_CRIT(MN)>=T0
                         TT_CRIT(MN)=T0;
                     else