Changeset 33279f7 in flexpart.git

Timestamp:

Feb 10, 2017, 2:47:25 PM (7 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:

9287c01

Parents:

6473ad3

Message:

counting how often xscav is less than 1

File:

• src/conccalc.f90 (modified) (4 diffs)

src/conccalc.f90

@@ -59 +59 @@
   real :: xl,yl,wx,wy,w
   real,parameter :: factor=.596831, hxmax=6.0, hymax=4.0, hzmax=150.
-
+  integer xscav_count
 
   ! For forward simulations, make a loop over the number of species;
@@ -65 +65 @@
   ! releasepoints
   !***************************************************************************
-
+  xscav_count=0
   do i=1,numpart
     if (itra1(i).ne.itime) goto 20
@@ -76 +76 @@
 33   continue
 
-
+  if (xscav_frac1(i,1).lt.0) xscav_count=xscav_count+1
+           
   ! For special runs, interpolate the air density to the particle position
   !************************************************************************
@@ -438 +439 @@
 20  continue
   end do
+  write(*,*) 'xscav count:',xscav_count
 
   !***********************************************************************