Changeset 18adf60 in flexpart.git for src/timemanager_mpi.f90

Timestamp:

May 9, 2016, 9:35:01 AM (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:

3b80e98, 93786a1

Parents:

62e65c7

Message:

Bug found in the MQUASILAG=1 option, xmass array accessed with out-of-bounds index. Quick fix implemented (see changes in timemanager and advance).

File:

• src/timemanager_mpi.f90 (modified) (2 diffs)

src/timemanager_mpi.f90

@@ -546 +546 @@
 ! Decide whether to write an estimate of the number of particles released, 
 ! or exact number (require MPI reduce operation)
-        numpart_tot_mpi = numpart*mp_partgroup_np
-
-        if (mp_exact_numpart.and..not.(lmpreader.and.lmp_use_reader)) then
+        if (mp_dev_mode) then
+          numpart_tot_mpi = numpart
+        else
+          numpart_tot_mpi = numpart*mp_partgroup_np
+        end if
+
+        if (mp_exact_numpart.and..not.(lmpreader.and.lmp_use_reader).and.&
+             &.not.mp_dev_mode) then
           call MPI_Reduce(numpart, numpart_tot_mpi, 1, MPI_INTEGER, MPI_SUM, id_root, &
                & mp_comm_used, mp_ierr)
@@ -738 +743 @@
             endif
 
-
-            if (mdomainfill.eq.0) then
+! Skip check on mass fraction when npoint represents particle number
+            if (mdomainfill.eq.0.and.mquasilag.eq.0) then
               if (xmass(npoint(j),ks).gt.0.)then 
                    xmassfract=max(xmassfract,real(npart(npoint(j)))* &