Changeset 1228ef7 in flexpart.git for src/releaseparticles_mpi.f90

Timestamp:

Aug 2, 2021, 12:11:23 PM (3 years ago)

Author:

Espen Sollum <eso@…>

Branches:

dev

Parents:

6cb0801

Message:

MPI: fix for mquasilag output

File:

• src/releaseparticles_mpi.f90 (modified) (4 diffs)

src/releaseparticles_mpi.f90

@@ -51 +51 @@
   real(kind=dp) :: juldate,julmonday,jul,jullocal,juldiff
   real,parameter :: eps=nxmax/3.e5,eps2=1.e-6
-  integer :: mind2
+  integer :: mind2, numpartcount_mpi
 ! mind2        eso: pointer to 2nd windfield in memory
-
   integer :: idummy = -7
 !save idummy,xmasssave
@@ -60 +59 @@
   logical :: first_call=.true.
 
-! Use different seed for each process.
-!****************************************************************************
+! Use different random seed for each process
+!*******************************************
   if (first_call) then
     idummy=idummy+mp_seed
@@ -68 +67 @@
 
   mind2=memind(2)
+
+! For mquasilag=1, assign unique particle ID across processes
+  numpartcount_mpi=mp_partid-mp_partgroup_np+1
 
 ! Determine the actual date and time in Greenwich (i.e., UTC + correction for daylight savings time)
@@ -210 +212 @@
                  nclassunc)
             numparticlecount=numparticlecount+1
+! Use a stride equal to number of processes for the MPI version
+            numpartcount_mpi=numpartcount_mpi+mp_partgroup_np
             if (mquasilag.eq.0) then
               npoint(ipart)=i
             else
-              npoint(ipart)=numparticlecount
+              npoint(ipart)=numpartcount_mpi
             endif
             idt(ipart)=mintime               ! first time step