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

Timestamp:

May 5, 2008 2:09:05 PM (16 years ago)

Author:

raasch

Message:

bugfixes for turbulent inflow in init_pegrid, inflow_turbulence, and init_3d_model

File:

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

palm/trunk/SOURCE/init_pegrid.f90

@@ -48 +48 @@
     USE control_parameters
     USE fft_xy
+    USE grid_variables
     USE indices
     USE pegrid
@@ -58 +59 @@
     IMPLICIT NONE
 
-    INTEGER ::  gathered_size, i, id_inflow_l, ind(5), j, k,                 &
+    INTEGER ::  gathered_size, i, id_inflow_l, id_recycling_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, &
@@ -922 +923 @@
 !-- Broadcast the id of the inflow PE
     IF ( inflow_l )  THEN
-       id_inflow_l = myid
+       id_inflow_l = myidx
     ELSE
        id_inflow_l = 0
     ENDIF
     CALL MPI_ALLREDUCE( id_inflow_l, id_inflow, 1, MPI_INTEGER, MPI_SUM, &
+                        comm1dx, ierr )
+
+!
+!-- Broadcast the id of the recycling plane
+!-- WARNING: needs to be adjusted in case of inflows other than from left side!
+    IF ( ( recycling_width / dx ) >= nxl  .AND.  ( recycling_width / dx ) <= nxr ) &
+    THEN
+       id_recycling_l = myidx
+    ELSE
+       id_recycling_l = 0
+    ENDIF
+    CALL MPI_ALLREDUCE( id_recycling_l, id_recycling, 1, MPI_INTEGER, MPI_SUM, &
                         comm1dx, ierr )