Changeset 2599 for palm/trunk/SOURCE/parin.f90

Timestamp:

Nov 1, 2017 1:18:45 PM (7 years ago)

Author:

hellstea

Message:

i/o grouping update for nested runs

File:

• palm/trunk/SOURCE/parin.f90 (modified) (5 diffs)

palm/trunk/SOURCE/parin.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! The i/o grouping is updated to work correctly also in nested runs.
+! 
+! 2575 2017-10-24 09:57:58Z maronga
 ! Renamed phi -> latitude, added longitude
 ! 
@@ -390 +393 @@
     INTEGER(iwp) ::  i      !<
     INTEGER(iwp) ::  ioerr  !< error flag for open/read/write
-
+    INTEGER(iwp) ::  myworldid       !<
+    INTEGER(iwp) ::  numworldprocs   !<
 
     NAMELIST /inipar/  aerosol_bulk, alpha_surface, approximation, bc_e_b,     &
@@ -513 +517 @@
 
     CALL location_message( 'finished', .TRUE. )
-
 !
 !-- Calculate the number of groups into which parallel I/O is split.
@@ -522 +525 @@
 !-- system.
 !-- First, set the default:
+    CALL MPI_COMM_RANK( MPI_COMM_WORLD, myworldid, ierr ) 
+    CALL MPI_COMM_SIZE( MPI_COMM_WORLD, numworldprocs, ierr )
     IF ( maximum_parallel_io_streams == -1  .OR.                               &
-         maximum_parallel_io_streams > numprocs )  THEN
-       maximum_parallel_io_streams = numprocs
+         maximum_parallel_io_streams > numworldprocs )  THEN
+       maximum_parallel_io_streams = numworldprocs
     ENDIF
 !
@@ -532 +537 @@
 !-- These settings are repeated in init_pegrid for the communicator comm2d,
 !-- which is not available here
-    io_blocks = numprocs / maximum_parallel_io_streams
-    io_group  = MOD( myid+1, io_blocks )
-
+    !io_blocks = numprocs / maximum_parallel_io_streams
+    io_blocks = numworldprocs / maximum_parallel_io_streams
+    !io_group  = MOD( myid+1, io_blocks )
+    io_group  = MOD( myworldid+1, io_blocks )
+    
     CALL location_message( 'reading NAMELIST parameters from PARIN', .FALSE. )
 !