Changeset d4a4250 in flexpart.git for src

Timestamp:

Oct 20, 2016, 2:26:33 PM (8 years ago)

Author:

Sabine <sabine.eckhardt@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie

Children:

0ee2193

Parents:

15fdb5c

Message:

counting how many particles removed at first ts

File:

• src/wetdepo.f90 (modified) (3 diffs)

src/wetdepo.f90

@@ -90 +90 @@
   real :: frac_act, liq_frac, dquer_m
 
-  integer :: blc_count, inc_count
+  integer :: blc_count, inc_count, firsttimerem
   real    :: Si_dummy, wetscav_dummy
   logical :: readclouds_this_nest
 
-
+  firsttimerem=0
 ! Compute interval since radioactive decay of deposited mass was computed
 !************************************************************************
@@ -436 +436 @@
                 xscav_frac1(jpart,ks)=0.
              else
+                 if (xmass1(jpart,ks).eq.0) then
+                     firsttimerem=firsttimerem+1
+                endif
                 xscav_frac1(jpart,ks)=xscav_frac1(jpart,ks)*(-1.)* &
                    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
@@ -464 +467 @@
   tot_inc_count=tot_inc_count+inc_count
 
+!  write (*,*) 'First time removed: ',firsttimerem
+
 end subroutine wetdepo