Changeset f75967d in flexpart.git for src/timemanager.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/timemanager.f90 (modified) (3 diffs)

src/timemanager.f90

@@ -138 +138 @@
 !ZHG 2015
 !CGZ-lifetime: set lifetime to 0
-  checklifetime(:,:)=0
-  species_lifetime(:,:)=0
-  print*, 'Initialized lifetime'
+  ! checklifetime(:,:)=0
+  ! species_lifetime(:,:)=0
+  ! print*, 'Initialized lifetime'
 !CGZ-lifetime: set lifetime to 0
   
@@ -422 +422 @@
         !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
+        ! 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
@@ -594 +594 @@
 !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
+                   ! 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
+                       ! 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