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

Timestamp:

Aug 20, 2019 11:55:33 AM (5 years ago)

Author:

oliver.maas

Message:

Added optional method for recycling of absolute values of pt and q instead of recycling only the deviations from the mean profile at the recycling plane. With the new method two problems are solved: 1. A horizontally homogeneous temperature and humidity field is achieved, even for non-neutral boundary layers and thus no thermal circulation is triggered. 2. No gravity waves build up at the inflow due to cyclic boundary conditions for pt and q.

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! added conversion from recycle_absolute_quantities to raq for recycling of 
+! absolute quantities and added error message PA184 for not implemented quantities
+! 
+! 4142 2019-08-05 12:38:31Z suehring
 ! Consider spinup in number of output timesteps for averaged 2D output (merge
 ! from branch resler).
@@ -866 +870 @@
     INTEGER(iwp) ::  k                               !< loop index
     INTEGER(iwp) ::  kk                              !< loop index
+    INTEGER(iwp) ::  r                               !< loop index
     INTEGER(iwp) ::  mid                             !< masked output running index
     INTEGER(iwp) ::  netcdf_data_format_save         !< initial value of netcdf_data_format
@@ -3647 +3652 @@
           CALL message( 'check_parameters', 'PA0421', 1, 2, 0, 6, 0 )
        ENDIF
+
+!
+!--    Convert recycle_absolute_quantities (list of strings that define the quantities for 
+!--    absolute recycling) to raq (list of logicals with length 7 corresponding to u,v,w,pt,e,q,s). 
+!--    Output error message for not implemented quantities.
+       DO r = LBOUND( recycle_absolute_quantities, 1 ), UBOUND( recycle_absolute_quantities, 1 )
+          SELECT CASE ( TRIM( recycle_absolute_quantities(r) ) )
+             CASE ( 'theta' )
+                raq(4) = .TRUE.
+             CASE ( 'q' )
+                raq(6) = .TRUE.
+             CASE ( '' )
+                CONTINUE
+             CASE DEFAULT
+                message_string = 'absolute recycling not implemented for variable ' // &
+                TRIM( recycle_absolute_quantities(r) )
+                CALL message( 'inflow_turbulence', 'PA184', 1, 2, 0, 6, 0 )
+          END SELECT
+       ENDDO
     ENDIF