Commit d4a42507 authored by Sabine's avatar Sabine
Browse files

counting how many particles removed at first ts

parent 15fdb5c1
...@@ -89,11 +89,11 @@ subroutine wetdepo(itime,ltsample,loutnext) ...@@ -89,11 +89,11 @@ subroutine wetdepo(itime,ltsample,loutnext)
real, parameter :: bcls(6) = (/22.7, 0.0, 0.0, 1321.0, 381.0, 0.0/) !now (Kyro et al 2009) 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, dquer_m
integer :: blc_count, inc_count integer :: blc_count, inc_count, firsttimerem
real :: Si_dummy, wetscav_dummy real :: Si_dummy, wetscav_dummy
logical :: readclouds_this_nest logical :: readclouds_this_nest
firsttimerem=0
! Compute interval since radioactive decay of deposited mass was computed ! Compute interval since radioactive decay of deposited mass was computed
!************************************************************************ !************************************************************************
...@@ -435,9 +435,12 @@ subroutine wetdepo(itime,ltsample,loutnext) ...@@ -435,9 +435,12 @@ subroutine wetdepo(itime,ltsample,loutnext)
xmass1(jpart,ks)=0. xmass1(jpart,ks)=0.
xscav_frac1(jpart,ks)=0. xscav_frac1(jpart,ks)=0.
else else
if (xmass1(jpart,ks).eq.0) then
firsttimerem=firsttimerem+1
endif
xscav_frac1(jpart,ks)=xscav_frac1(jpart,ks)*(-1.)* & xscav_frac1(jpart,ks)=xscav_frac1(jpart,ks)*(-1.)* &
wetdeposit(ks)/xmass1(jpart,ks) wetdeposit(ks)/xmass1(jpart,ks)
! write (*,*) 'paricle kept: ',jpart,ks,wetdeposit(ks),xscav_frac1(jpart,ks) ! write (*,*) 'paricle kept: ',jpart,ks,wetdeposit(ks),xscav_frac1(jpart,ks),xmass(jpart,ks)
endif endif
endif endif
endif endif
...@@ -463,4 +466,6 @@ subroutine wetdepo(itime,ltsample,loutnext) ...@@ -463,4 +466,6 @@ subroutine wetdepo(itime,ltsample,loutnext)
tot_blc_count=tot_blc_count+blc_count tot_blc_count=tot_blc_count+blc_count
tot_inc_count=tot_inc_count+inc_count tot_inc_count=tot_inc_count+inc_count
! write (*,*) 'First time removed: ',firsttimerem
end subroutine wetdepo end subroutine wetdepo
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