Don't set t_soisno below freezing when FATES is on

olyson · Feb 20, 2024 · 04167dc · 04167dc
1 parent e4dde90
commit 04167dc
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/biogeophys/TemperatureType.F90 b/src/biogeophys/TemperatureType.F90
@@ -741,7 +741,11 @@ subroutine InitCold(this, bounds, &
                   end if
                end if
             else
-               this%t_soisno_col(c,1:nlevgrnd) = 272._r8
+               if (use_fates) then
+                  this%t_soisno_col(c,1:nlevgrnd) = 274._r8
+               else
+                  this%t_soisno_col(c,1:nlevgrnd) = 272._r8
+               end if
                if (use_excess_ice .and. (lun%itype(l) == istsoil .or. lun%itype(l) == istcrop)) then
                   this%t_soisno_col(c,1:nlevgrnd) = SHR_CONST_TKFRZ - 5.0_r8 !needs to be below freezing to properly initiate excess ice
                end if