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

Timestamp:

Nov 8, 2016 3:00:55 PM (8 years ago)

Author:

gronemeier

Message:

Implement turbulent outflow condition

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Implement turbulent outflow condition
 ! 
 ! Former revisions:
@@ -3703 +3703 @@
     ENDIF
 
+
+    IF ( turbulent_outflow )  THEN
+!
+!--    Turbulent outflow requires Dirichlet conditions at the respective inflow
+!--    boundary (so far, a turbulent outflow is realized at the right side only)
+       IF ( bc_lr /= 'dirichlet/radiation' )  THEN
+          message_string = 'turbulent_outflow = .T. requires ' //              &
+                           'bc_lr = "dirichlet/radiation"'
+          CALL message( 'check_parameters', 'PA0038', 1, 2, 0, 6, 0 )
+       ENDIF
+!
+!--    The ouflow-source plane must lay inside the model domain
+       IF ( outflow_source_plane < dx  .OR.  &
+            outflow_source_plane > nx * dx )  THEN
+          WRITE( message_string, * )  'illegal value for outflow_source'//     &
+                                      '_plane: ', outflow_source_plane
+          CALL message( 'check_parameters', 'PA0145', 1, 2, 0, 6, 0 )
+       ENDIF
+    ENDIF
+
 !
 !-- Determine damping level index for 1D model