Changeset d6a0977 in flexpart.git for src/timemanager.f90

Timestamp:

Dec 14, 2015, 3:10:04 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:

f75967d

Parents:

88d8c3d

Message:

Updates to Henrik's wet depo scheme

File:

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

src/timemanager.f90

@@ -136 +136 @@
   !**********************************************************************
 
+!ZHG 2015
+!CGZ-lifetime: set lifetime to 0
+  checklifetime(:,:)=0
+  species_lifetime(:,:)=0
+  print*, 'Initialized lifetime'
+!CGZ-lifetime: set lifetime to 0
+  
+
 
   if (verbosity.gt.0) then
@@ -411 +419 @@
         if (iflux.eq.1) call fluxoutput(itime)
         write(*,45) itime,numpart,gridtotalunc,wetgridtotalunc,drygridtotalunc
+ 
+        !CGZ-lifetime: output species lifetime
+!ZHG
+        write(*,*) 'Overview species lifetime in days', &
+             real((species_lifetime(:,1)/species_lifetime(:,2))/real(3600.0*24.0))
+        write(*,*) 'all info:',species_lifetime
+!ZHG
+        !CGZ-lifetime: output species lifetime
+
         !write(*,46) float(itime)/3600,itime,numpart
 45      format(i9,' SECONDS SIMULATED: ',i8, ' PARTICLES:    Uncertainty: ',3f7.3)
@@ -575 +592 @@
                    xmassfract=max(xmassfract,real(npart(npoint(j)))* &
                    xmass1(j,ks)/xmass(npoint(j),ks))
+!ZHG 2015
+                  !CGZ-lifetime: Check mass fraction left/save lifetime
+                   if(real(npart(npoint(j)))*xmass1(j,ks)/xmass(npoint(j),ks).lt.0.01.and.checklifetime(j,ks).eq.0.)then
+                       !Mass below 1% of initial >register lifetime
+                       checklifetime(j,ks)=abs(itra1(j)-itramem(j))
+                       species_lifetime(ks,1)=species_lifetime(ks,1)+abs(itra1(j)-itramem(j))
+                       species_lifetime(ks,2)= species_lifetime(ks,2)+1
+                   endif
+                   !CGZ-lifetime: Check mass fraction left/save lifetime
+!ZHG 2015
             else
               xmassfract=1.