Changeset 3884

Timestamp:

Apr 10, 2019 1:31:55 PM (6 years ago)

Author:

Giersch

Message:

id_recycling is only calculated in case of tubulent inflow

File:

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

palm/trunk/SOURCE/init_pegrid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! id_recycling is only calculated in case of tubulent inflow
+! 
+! 3761 2019-02-25 15:31:42Z raasch
 ! unused variable removed
 ! 
@@ -263 +266 @@
                message_string, mg_switch_to_pe0_level, momentum_advec,         &
                psolver, outflow_source_plane, recycling_width, scalar_advec,   &
-               subdomain_size, turbulent_outflow, y_shift
+               subdomain_size, turbulent_inflow, turbulent_outflow, y_shift
 
     USE grid_variables,                                                        &
@@ -1352 +1355 @@
 !-- Broadcast the id of the recycling plane
 !-- WARNING: needs to be adjusted in case of inflows other than from left side!
-    IF ( NINT( recycling_width / dx ) >= nxl  .AND. &
-         NINT( recycling_width / dx ) <= nxr )  THEN
-       id_recycling_l = myidx
-    ELSE
-       id_recycling_l = 0
-    ENDIF
-    IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
-    CALL MPI_ALLREDUCE( id_recycling_l, id_recycling, 1, MPI_INTEGER, MPI_SUM, &
-                        comm1dx, ierr )
+    IF ( turbulent_inflow ) THEN
+    
+       IF ( NINT( recycling_width / dx ) >= nxl  .AND. &
+            NINT( recycling_width / dx ) <= nxr )  THEN
+          id_recycling_l = myidx
+       ELSE
+          id_recycling_l = 0
+       ENDIF
+       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
+       CALL MPI_ALLREDUCE( id_recycling_l, id_recycling, 1, MPI_INTEGER, MPI_SUM, &
+                           comm1dx, ierr )
+                           
+    ENDIF
 
 !