Merge pull request #72 from EcoExtreML/fix_Tmin

Fix tmin and grasslands
EcoExtreML · Jun 30, 2022 · dea100f · dea100f
2 parents 5dc0a82 + 5d9761a
commit dea100f
Show file tree

Hide file tree

Showing 4 changed files with 9 additions and 4 deletions.
diff --git a/exe/STEMMUS_SCOPE b/exe/STEMMUS_SCOPE
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

diff --git a/src/Initial_root_biomass.m b/src/Initial_root_biomass.m
@@ -25,7 +25,7 @@
     beta = 0.943; 
 elseif strcmp(IGBP_veg_long(1:14)', 'Woody Savannas') 
     beta = 0.943; 
-elseif strcmp(IGBP_veg_long(1:10)', 'Grasslands') 
+elseif strcmp(IGBP_veg_long(1:9)', 'Grassland') 
     beta = 0.943; 
 else 
     beta = 0.943; 

diff --git a/src/STEMMUS_SCOPE.m b/src/STEMMUS_SCOPE.m
@@ -269,7 +269,7 @@
     V(10).Val= [9];
     V(11).Val= [0];
     V(28).Val= [0.03];
-elseif strcmp(IGBP_veg_long(1:10)', 'Grasslands')  
+elseif strcmp(IGBP_veg_long(1:9)', 'Grassland')  
     V(14).Val= [0.2 0.3 288 303 328];
     if isequal(sitename1,{'AR-SLu'})||isequal(sitename1,{'AU-Ync'})||isequal(sitename1,{'CH-Oe1'})||isequal(sitename1,{'DK-Lva'})||isequal(sitename1,{'US-AR1'})||isequal(sitename1,{'US-AR2'})||isequal(sitename1,{'US-Aud'})||isequal(sitename1,{'US-SRG'})
         V(9).Val= [120];