diff --git a/src/AtmosForcingMod.F90 b/src/AtmosForcingMod.F90
index 94451a74..9209a2dd 100644
--- a/src/AtmosForcingMod.F90
+++ b/src/AtmosForcingMod.F90
@@ -154,7 +154,7 @@ subroutine ProcessAtmosForcing(noahmp)
     ! wet-bulb scheme (Wang et al., 2019 GRL), C.He, 12/18/2020, R. Abolafia-Rosnezweig, 02/01/2024
     if ( OptRainSnowPartition == 5 ) then
 
-        if ( TemperatureAirRefHeight >= (ConstFreezePoint+10) ) then !avoid numerical errors when temperature is high
+        if ( TemperatureAirRefHeight >= (ConstFreezePoint+8) ) then !avoid numerical errors when temperature is high
             FrozenPrecipFrac = 0.0
         else
             TemperatureDegC = min( 50.0, max(-50.0,(TemperatureAirRefHeight-ConstFreezePoint)) )    ! Kelvin to degree Celsius with limit -50 to +50
diff --git a/src/PhenologyMainMod.F90 b/src/PhenologyMainMod.F90
index e2319350..e455a032 100644
--- a/src/PhenologyMainMod.F90
+++ b/src/PhenologyMainMod.F90
@@ -111,7 +111,7 @@ subroutine PhenologyMain (noahmp)
     ThicknessCanBury      = min(max(SnowDepth-HeightCanopyBot,0.0), (HeightCanopyTop-HeightCanopyBot))
     CanopyFracSnowBury    = ThicknessCanBury / max(1.0e-06, (HeightCanopyTop-HeightCanopyBot))           ! snow buried fraction
     if ( (HeightCanopyTop > 0.0) .and. (HeightCanopyTop <= 1.0) ) then                                   ! MB: change to 1.0 & 0.2 to reflect changes to HeightCanopyTop in MPTABLE
-       SnowDepthVegBury   = HeightCanopyTop * exp(-SnowDepth / 0.2)
+       SnowDepthVegBury   = HeightCanopyTop * exp(-min(SnowDepth,10.0) / 0.2)     ! CH: add limit to avoid numerical floating issue
        CanopyFracSnowBury = min(SnowDepth, SnowDepthVegBury) / SnowDepthVegBury
     endif