Changeset f75967d in flexpart.git for src/wetdepo.f90

Timestamp:

Jan 5, 2016, 12:20:01 PM (8 years ago)

Author:

Espen Sollum ATMOS <eso@…>

Branches:

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

Children:

d8107c2

Parents:

d6a0977

Message:

Fixed a bug with number of occurances of wet scavenging written to stdout

File:

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

src/wetdepo.f90

@@ -73 +73 @@
   integer :: ngrid,itage,nage,hz,il,interp_time, n, clouds_v
   integer :: ks, kp
-  integer :: n1,n2, icbot,ictop, indcloud !TEST
+!  integer :: n1,n2, icbot,ictop, indcloud !TEST
   real :: S_i, act_temp, cl, cle ! in cloud scavenging
   real :: clouds_h ! cloud height for the specific grid point
@@ -105 +105 @@
   ! Loop over all particles
   !************************
-
 
   blc_count=0
@@ -185 +184 @@
     n=memind(2)
     if (abs(memtime(1)-interp_time).lt.abs(memtime(2)-interp_time)) &
-    n=memind(1)
+         n=memind(1)
 
     if (ngrid.eq.0) then
@@ -199 +198 @@
   ! scavenging is done
 
-     if (clouds_v.le.1) goto 20
+    if (clouds_v.le.1) goto 20
   
   ! 1) Parameterization of the the area fraction of the grid cell where the
@@ -311 +310 @@
 ! NIK 13 may 2015: only do incloud if positive in-cloud scavenging parameters are given in species file
         if (weta_in(ks).gt.0. .or. wetb_in(ks).gt.0.) then 
+          inc_count=inc_count+1
 ! if negative coefficients (turned off) set to zero for use in equation
           if (weta_in(ks).lt.0.) weta_in(ks)=0.
@@ -355 +355 @@
            !OLD 
           if (readclouds) then
-           wetscav=S_i*(prec(1)/3.6E6)
+            wetscav=S_i*(prec(1)/3.6E6)
           else
-          wetscav=S_i*(prec(1)/3.6E6)/clouds_h
+            wetscav=S_i*(prec(1)/3.6E6)/clouds_h
           endif