diff --git a/ED/src/utils/therm_lib.f90 b/ED/src/utils/therm_lib.f90
index a42cb0849..06197d672 100644
--- a/ED/src/utils/therm_lib.f90
+++ b/ED/src/utils/therm_lib.f90
@@ -2247,13 +2247,14 @@ real(kind=4) function vpdefil(pres,temp,humi,is_shv,useice)
 
 
       !------------------------------------------------------------------------------------!
-      !     Find the vapour pressure and the saturation vapour pressure.                   !
+      !     Find the vapour pressure and the saturation vapour pressure. Use different     !
+      ! functions for ice VS liquid saturation vapour pressure.                            !
       !------------------------------------------------------------------------------------!
       pvap = ( pres * shv ) / ( ep + (1.0 - ep) * shv )
       if (frozen) then
          psat = esif(temp)
       else
-         psat = esif(temp)
+         psat = eslf(temp)
       end if
       !------------------------------------------------------------------------------------!