Changeset 1228ef7 in flexpart.git for src/mpi_mod.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/mpi_mod.f90 (modified) (4 diffs)

src/mpi_mod.f90

@@ -178 +178 @@
     implicit none
 
-    integer :: i,j,s,addmaxpart=0
+    integer :: i,j,s
 
 ! Each process gets an ID (mp_pid) in the range 0,..,mp_np-1
@@ -486 +486 @@
 
     integer,intent(in) :: num_part
-    integer :: i,jj, addone
+    integer :: i
 
 ! Exit if too many particles
@@ -589 +589 @@
     integer, intent(in) :: itime
     real :: pmin,z
-    integer :: i,jj,nn, num_part=1,m,imin, num_trans
-    logical :: first_iter
+    integer :: i,jj,nn,m,imin, num_trans
     integer,allocatable,dimension(:) :: idx_arr
     real,allocatable,dimension(:) :: sorted ! TODO: we don't really need this
@@ -613 +612 @@
     end do
 
-! Do not rebalance particles for ipout=3    
-    if (ipout.eq.3) return
+! Do not rebalance particles for ipout=3 or mquasilag=1
+    if ((ipout.eq.3).or.(mquasilag.eq.1)) return
 
 ! For each successive element in index array, see if a lower value exists