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

Timestamp:

Jan 9, 2008 8:17:38 AM (17 years ago)

Author:

raasch

Message:

second preliminary update for turbulent inflow

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
 ! TEST OUTPUT (TO BE REMOVED) logging mpi2 ierr values
 !
@@ -458 +458 @@
     DEALLOCATE( nxlf , nxrf , nynf , nysf )
 
+!
+!-- Collect index bounds from other PEs (to be written to restart file later)
+    ALLOCATE( hor_index_bounds(4,0:numprocs-1) )
+
+    IF ( myid == 0 )  THEN
+
+       hor_index_bounds(1,0) = nxl
+       hor_index_bounds(2,0) = nxr
+       hor_index_bounds(3,0) = nys
+       hor_index_bounds(4,0) = nyn
+
+!
+!--    Receive data from all other PEs
+       DO  i = 1, numprocs-1
+          CALL MPI_RECV( ibuf, 4, MPI_INTEGER, i, MPI_ANY_TAG, comm2d, status, &
+                         ierr )
+          hor_index_bounds(:,i) = ibuf(1:4)
+       ENDDO
+
+    ELSE
+!
+!--    Send index bounds to PE0
+       ibuf(1) = nxl
+       ibuf(2) = nxr
+       ibuf(3) = nys
+       ibuf(4) = nyn
+       CALL MPI_SEND( ibuf, 4, MPI_INTEGER, 0, myid, comm2d, ierr )
+
+    ENDIF
+
 #if defined( __print )
 !
@@ -606 +636 @@
     nzta = nz
     nnz  = nz
+
+    ALLOCATE( hor_index_bounds(4,0:0) )
+    hor_index_bounds(1,0) = nxl
+    hor_index_bounds(2,0) = nxr
+    hor_index_bounds(3,0) = nys
+    hor_index_bounds(4,0) = nyn
 
 !