Changeset 151 for palm/trunk/SOURCE/init_pegrid.f90

Timestamp:

Mar 7, 2008 1:42:18 PM (16 years ago)

Author:

raasch

Message:

preliminary update for the turbulence recycling method

File:

• palm/trunk/SOURCE/init_pegrid.f90 (modified) (3 diffs)

palm/trunk/SOURCE/init_pegrid.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! Collect on PE0 horizontal index bounds from all other PEs
+! Collect on PE0 horizontal index bounds from all other PEs,
+! broadcast the id of the inflow PE (using the respective communicator)
 ! TEST OUTPUT (TO BE REMOVED) logging mpi2 ierr values
 !
@@ -57 +58 @@
     IMPLICIT NONE
 
-    INTEGER ::  gathered_size, i, ind(5), j, k, maximum_grid_level_l,     &
-                mg_switch_to_pe0_level_l, mg_levels_x, mg_levels_y,       &
-                mg_levels_z, nnx_y, nnx_z, nny_x, nny_z, nnz_x, nnz_y,    &
-                numproc_sqr, nx_total, nxl_l, nxr_l, nyn_l, nys_l, nzb_l, &
-                nzt_l, omp_get_num_threads, subdomain_size
+    INTEGER ::  gathered_size, i, id_inflow_l, ind(5), j, k,                 &
+                maximum_grid_level_l, mg_switch_to_pe0_level_l, mg_levels_x, &
+                mg_levels_y, mg_levels_z, nnx_y, nnx_z, nny_x, nny_z, nnz_x, &
+                nnz_y, numproc_sqr, nx_total, nxl_l, nxr_l, nyn_l, nys_l,    &
+                nzb_l, nzt_l, omp_get_num_threads, subdomain_size
 
     INTEGER, DIMENSION(:), ALLOCATABLE ::  ind_all, nxlf, nxrf, nynf, nysf
@@ -918 +919 @@
     ENDIF
 
+!
+!-- Broadcast the id of the inflow PE
+    IF ( inflow_l )  THEN
+       id_inflow_l = myid
+    ELSE
+       id_inflow_l = 0
+    ENDIF
+    CALL MPI_ALLREDUCE( id_inflow_l, id_inflow, 1, MPI_INTEGER, MPI_SUM, &
+                        comm1dx, ierr )
+
 #else
     IF ( bc_lr == 'dirichlet/radiation' )  THEN