Changeset 4297 for palm/trunk/SOURCE/inflow_turbulence.f90

Timestamp:

Nov 21, 2019 10:37:50 AM (5 years ago)

Author:

oliver.maas

Message:

In case of turbulent inflow: Adjusted recycling_yshift, so that y-shift can be given as multiples of PE (as in y_shift for cyclic BCs).

File:

• palm/trunk/SOURCE/inflow_turbulence.f90 (modified) (2 diffs)

palm/trunk/SOURCE/inflow_turbulence.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! changed recycling_yshift so that the y-shift can be a multiple of PE
+! instead of y-shift of a half domain width
+! 
+! 4183 2019-08-23 07:33:16Z oliver.maas
 ! simplified steering of recycling of absolute values by initialization 
 ! parameter recycling_method_for_thermodynamic_quantities
@@ -240 +244 @@
 !
 !-- y-shift for inflow_dist
-!-- Shift inflow_dist in positive y direction by a distance of INT( npey / 2 )
-    IF ( recycling_yshift .AND. myidx == id_inflow ) THEN
+!-- Shift inflow_dist in positive y direction by a number of
+!-- PEs equal to recycling_yshift
+    IF ( ( recycling_yshift /= 0 ) .AND. myidx == id_inflow ) THEN
+
 !
 !--    Calculate the ID of the PE which sends data to this PE (prev) and of the
 !--    PE which receives data from this PE (next).
-       IF ( myidy >= INT( pdims(2) / 2 ) ) THEN
-          prev = myidy - INT( pdims(2) / 2 )
-       ELSE
-          prev = pdims(2) - ( INT( pdims(2) / 2 ) - myidy )
-       ENDIF
-     
-       IF ( myidy < pdims(2) - INT( pdims(2) / 2 ) ) THEN
-          next = myidy + INT( pdims(2) / 2 )
-       ELSE
-          next = INT( pdims(2) / 2 ) - ( pdims(2) - myidy )
-       ENDIF
-
+       prev = MODULO(myidy - recycling_yshift , pdims(2))
+       next = MODULO(myidy + recycling_yshift , pdims(2))
+       
        local_inflow_dist = 0.0_wp