Commit 56a2e54b authored by Sabine's avatar Sabine
Browse files

BUGFIX: correction gfraction in timemanager

parent 1c0d5e63
......@@ -568,7 +568,8 @@ subroutine timemanager
if ((xscav_frac1(j,ks).lt.0)) then
call get_wetscav(itime,lsynctime,loutnext,j,ks,grfraction,idummy,idummy,wetscav)
if (wetscav(ks).gt.0) then
xscav_frac1(j,ks)=wetscav(ks)*(zpoint2(npoint(j))-zpoint1(npoint(j)))
xscav_frac1(j,ks)=wetscav(ks)* &
(zpoint2(npoint(j))-zpoint1(npoint(j)))*grfraction(1)
else
xmass1(j,ks)=0.
xscav_frac1(j,ks)=0.
......
......@@ -112,8 +112,16 @@ subroutine wetdepo(itime,ltsample,loutnext)
if (wetscav.gt.0.) then
wetdeposit(ks)=xmass1(jpart,ks)* &
(1.-exp(-wetscav*abs(ltsample)))*grfraction(1) ! wet deposition
! if (xscav_frac1(jpart,ks).lt.0) then
! xscav_frac1(jpart,ks)=wetscav*grfraction(1)* &
! (zpoint2(npoint(jpart))-zpoint1(npoint(jpart)))
! endif
else ! if no scavenging
wetdeposit(ks)=0.
! if (xscav_frac1(jpart,ks).lt.0) then
! xscav_frac1(jpart,ks)=0.
! xmass1(jpart,ks)=0.
! endif
endif
restmass = xmass1(jpart,ks)-wetdeposit(ks)
......
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