From bde771c14ab9f9308a6a3f3e196d41699e0c83ee Mon Sep 17 00:00:00 2001 From: Yixin Ma Date: Thu, 11 Dec 2025 12:40:00 -0500 Subject: [PATCH] Updates to the function vpdefil in therm_lib.f90: use different functions to calculate saturation vapor pressure for ice VS liquid. The old scripts always use ice saturation vapor pressure. --- ED/src/utils/therm_lib.f90 | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) 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 !------------------------------------------------------------------------------------!