Merge pull request #424 from yixin98/vpdefil_update

Updates to VPD calculation

Author: mpaiao

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
       !------------------------------------------------------------------------------------!