Commit 964f28f5 authored by Espen Sollum's avatar Espen Sollum
Browse files

Updated wetdepo from Sabine

parent 4f5fa008
......@@ -86,7 +86,7 @@ subroutine wetdepo(itime,ltsample,loutnext)
!ZHG aerosol below-cloud scavenging removal polynomial constants for rain and snow
real, parameter :: bclr(6) = (/274.35758, 332839.59273, 226656.57259, 58005.91340, 6588.38582, 0.244984/) !rain (Laakso et al 2003)
real, parameter :: bcls(6) = (/22.7, 0.0, 0.0, 1321.0, 381.0, 0.0/) !now (Kyro et al 2009)
real :: frac_act, liq_frac, dquer_m
real :: frac_act, liq_frac, ice_frac, dquer_m
integer :: blc_count, inc_count
real :: Si_dummy, wetscav_dummy
......@@ -321,15 +321,19 @@ subroutine wetdepo(itime,ltsample,loutnext)
if (in_aero(ks).lt.0.) in_aero(ks)=0.
!ZHG updated parameterization of cloud water to better reproduce the values coming from ECMWF
cl=1.6E-6*prec(1)**0.36
! cl=1.6E-6*prec(1)**0.36
cl=1E6*2E-7*prec(1)**0.36
!ZHG: Calculate the partition between liquid and water phase water.
if (act_temp .le. 253.) then
liq_frac=0
ice_frac=1
else if (act_temp .ge. 273.) then
liq_frac=1
ice_frac=0
else
liq_frac =((act_temp-273.)/(273.-253.))**2.
ice_frac =((act_temp-273.)/(273.-253.))**2.
liq_frac=max(0.,1.-ice_frac)
endif
! ZHG: Calculate the aerosol partition based on cloud phase and Ai and Bi
frac_act = liq_frac*ccn_aero(ks) +(1-liq_frac)*in_aero(ks)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment