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

Timestamp:

Nov 1, 2017 2:11:20 PM (7 years ago)

Author:

raasch

Message:

small changes concerning r2599, cycle number are now three digits wide

File:

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

palm/trunk/SOURCE/parin.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! some comments added and variables renamed concerning r2599
+! 
+! 2599 2017-11-01 13:18:45Z hellstea
 ! The i/o grouping is updated to work correctly also in nested runs.
 ! 
@@ -391 +394 @@
     IMPLICIT NONE
 
-    INTEGER(iwp) ::  i      !<
-    INTEGER(iwp) ::  ioerr  !< error flag for open/read/write
-    INTEGER(iwp) ::  myworldid       !<
-    INTEGER(iwp) ::  numworldprocs   !<
+    INTEGER(iwp) ::  global_id      !< process id with respect to MPI_COMM_WORLD
+    INTEGER(iwp) ::  global_procs   !< # of procs with respect to MPI_COMM_WORLD
+    INTEGER(iwp) ::  i              !<
+    INTEGER(iwp) ::  ioerr          !< error flag for open/read/write
 
     NAMELIST /inipar/  aerosol_bulk, alpha_surface, approximation, bc_e_b,     &
@@ -524 +527 @@
 !-- severe problems depending on the configuration of the underlying file
 !-- system.
+!-- Calculation of the number of blocks and the I/O group must be based on all
+!-- PEs involved in this run. Since myid and numprocs are related to the
+!-- comm2d communicator, which gives only a subset of all PEs in case of
+!-- nested runs, that information must be inquired again from the global
+!-- communicator.
 !-- First, set the default:
-    CALL MPI_COMM_RANK( MPI_COMM_WORLD, myworldid, ierr ) 
-    CALL MPI_COMM_SIZE( MPI_COMM_WORLD, numworldprocs, ierr )
+    CALL MPI_COMM_RANK( MPI_COMM_WORLD, global_id, ierr )
+    CALL MPI_COMM_SIZE( MPI_COMM_WORLD, global_procs, ierr )
     IF ( maximum_parallel_io_streams == -1  .OR.                               &
-         maximum_parallel_io_streams > numworldprocs )  THEN
-       maximum_parallel_io_streams = numworldprocs
+         maximum_parallel_io_streams > global_procs )  THEN
+       maximum_parallel_io_streams = global_procs
     ENDIF
 !
@@ -535 +543 @@
 !-- respective PE belongs. I/O of the groups is done in serial, but in parallel
 !-- for all PEs belonging to the same group.
-!-- These settings are repeated in init_pegrid for the communicator comm2d,
-!-- which is not available here
-    !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 )
+    io_blocks = global_procs / maximum_parallel_io_streams
+    io_group  = MOD( global_id+1, io_blocks )
     
     CALL location_message( 'reading NAMELIST parameters from PARIN', .FALSE. )