Commit a5579720 authored by Sabine's avatar Sabine

BUGFIX for wetdeposition, where icefraction and liquidfraction where interchanged

parent bf7ca87f
......@@ -335,16 +335,25 @@ subroutine wetdepo(itime,ltsample,loutnext)
cl=1.6E-6*prec(1)**0.36
endif
!ZHG: Calculate the partition between liquid and water phase water.
!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.
! sec bugfix after FLEXPART paper review, liq_frac was 1-liq_frac
! IP bugfix v10.4, calculate ice_frac and liq_frac
ice_frac= ((act_temp-273.)/(273.-253.))**2.
!liq_frac = 1-ice_frac !((act_temp-253.)/(273.-253.))**2.
liq_frac=max(0.,1.-ice_frac)
end if
! 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)
! frac_act = liq_frac*ccn_aero(ks) +(1-liq_frac)*in_aero(ks)
! IP, use ice_frac and liq_frac
frac_act = liq_frac*ccn_aero(ks) + ice_frac*in_aero(ks)
!ZHG Use the activated fraction and the liqid water to calculate the washout
......
Markdown is supported
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