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

Timestamp:

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

Author:

raasch

Message:

preliminary update for the turbulence recycling method

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! Case of reading data for recycling included in initializing_actions 
+! Case of reading data for recycling included in initializing_actions,
+! check of turbulent_inflow and calculation of recycling_plane
 !
 ! Former revisions:
@@ -2981 +2982 @@
 
 !
+!-- A turbulent inflow requires Dirichlet conditions at the respective inflow
+!-- boundary (so far, a turbulent inflow is realized from the left side only)
+    IF ( turbulent_inflow  .AND.  bc_lr /= 'dirichlet/radiation' )  THEN
+       IF ( myid == 0 )  THEN
+          PRINT*, '+++ check_parameters:'
+          PRINT*, '    turbulent_inflow = .T. requires a Dirichlet condition', &
+                  ' at the inflow boundary'
+       ENDIF
+       CALL local_stop
+    ENDIF
+
+!
+!-- In case of turbulent inflow calculate the index of the recycling plane
+    IF ( turbulent_inflow )  THEN
+       IF ( recycling_width == 9999999.9 )  THEN
+!
+!--       Set the default value for the width of the recycling domain
+          recycling_width = 0.1 * nx * dx
+       ELSE
+          IF ( recycling_width < dx  .OR.  recycling_width > nx * dx )  THEN
+             IF ( myid == 0 )  THEN
+                PRINT*, '+++ check_parameters:'
+                PRINT*, '    illegal value for recycling_width: ', &
+                        recycling_width
+             ENDIF
+             CALL local_stop
+          ENDIF
+       ENDIF
+!
+!--    Calculate the index
+       recycling_plane = recycling_width / dx
+    ENDIF
+
+!
 !-- Check random generator
     IF ( random_generator /= 'system-specific'  .AND. &