Changes in src/mpi_mod.f90 [92fab65:f3054ea] in flexpart.git

File:

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

src/mpi_mod.f90

@@ -1 @@
-! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
-! SPDX-License-Identifier: GPL-3.0-or-later
+!**********************************************************************
+! Copyright 1998,1999,2000,2001,2002,2005,2007,2008,2009,2010         *
+! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa,             *
+! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann   *
+!                                                                     *
+! This file is part of FLEXPART.                                      *
+!                                                                     *
+! FLEXPART is free software: you can redistribute it and/or modify    *
+! it under the terms of the GNU General Public License as published by*
+! the Free Software Foundation, either version 3 of the License, or   *
+! (at your option) any later version.                                 *
+!                                                                     *
+! FLEXPART is distributed in the hope that it will be useful,         *
+! but WITHOUT ANY WARRANTY; without even the implied warranty of      *
+! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the       *
+! GNU General Public License for more details.                        *
+!                                                                     *
+! You should have received a copy of the GNU General Public License   *
+! along with FLEXPART.  If not, see <http://www.gnu.org/licenses/>.   *
+!**********************************************************************
 
 module mpi_mod
@@ -2579 +2597 @@
 601 end subroutine mpif_tm_reduce_grid_nest
 
+
+  subroutine mpif_tm_reduce_initcond
+!***********************************************************************
+! Collect init_cond to PID 0, adding from all processes.
+!
+!
+!***********************************************************************
+    use com_mod
+    use unc_mod
+    use par_mod
+
+    implicit none
+
+    integer :: grid_size
+
+!**********************************************************************
+
+    grid_size=numxgrid*numygrid*numzgrid*maxspec* &
+         & maxpointspec_act
+
+! Time for MPI communications
+    if (mp_measure_time) call mpif_mtime('commtime',0)
+
+
+#ifdef USE_MPIINPLACE
+! Using in-place reduction
+    if (lroot) then
+      call MPI_Reduce(MPI_IN_PLACE, init_cond, grid_size, mp_sp, MPI_SUM, id_root, &
+           & mp_comm_used, mp_ierr)
+      if (mp_ierr /= 0) goto 600
+    else
+      call MPI_Reduce(init_cond, 0, grid_size, mp_sp, MPI_SUM, id_root, &
+           & mp_comm_used, mp_ierr)
+      if (mp_ierr /= 0) goto 600
+    end if
+
+#else
+    call MPI_Reduce(init_cond, init_cond0, grid_size, mp_sp, MPI_SUM, id_root, &
+         & mp_comm_used, mp_ierr)
+    if (mp_ierr /= 0) goto 600
+!    if (lroot) init_cond = init_cond0
+#endif
+
+    if (mp_measure_time) call mpif_mtime('commtime',1)
+
+    goto 601
+
+600 write(*,*) "mpi_mod> mp_ierr \= 0", mp_ierr
+    stop
+
+601 end subroutine mpif_tm_reduce_initcond
+    
 
   subroutine mpif_mtime(ident,imode)