Changeset 56a2e54 in flexpart.git

Timestamp:

Dec 6, 2016, 2:22:08 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:

6481010

Parents:

1c0d5e6

Message:

BUGFIX: correction gfraction in timemanager

Location:

src

Files:

• timemanager.f90 (modified) (1 diff)
• wetdepo.f90 (modified) (1 diff)

src/timemanager.f90

@@ -569 +569 @@
             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.

src/wetdepo.f90

@@ -113 +113 @@
         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