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

src/timemanager_mpi.f90

@@ -151 +151 @@
 
 !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
 
@@ -564 +564 @@
         
         !CGZ-lifetime: output species lifetime
-        ! if (lroot) then
+        if (lroot) then
         !   write(*,*) 'Overview species lifetime in days', &
         !        real((species_lifetime(:,1)/species_lifetime(:,2))/real(3600.0*24.0))
@@ -573 +573 @@
         !     write (*,*) 'timemanager> starting simulation'
         !   end if
-        ! end if
+        end if
 
 45      format(i13,' SECONDS SIMULATED: ',i13, ' PARTICLES:    Uncertainty: ',3f7.3)
@@ -757 +757 @@
                    
                    !CGZ-lifetime: Check mass fraction left/save lifetime
-                   if(lroot.and.real(npart(npoint(j)))*xmass1(j,ks)/xmass(npoint(j),ks).lt.0.01.and.checklifetime(j,ks).eq.0.)then
+                   ! if(lroot.and.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