Ignore:
Timestamp:
Aug 20, 2019 11:55:33 AM (2 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/check_parameters.f90

    r4142 r4172  
    2525! -----------------
    2626! $Id$
     27! added conversion from recycle_absolute_quantities to raq for recycling of
     28! absolute quantities and added error message PA184 for not implemented quantities
     29!
     30! 4142 2019-08-05 12:38:31Z suehring
    2731! Consider spinup in number of output timesteps for averaged 2D output (merge
    2832! from branch resler).
     
    866870    INTEGER(iwp) ::  k                               !< loop index
    867871    INTEGER(iwp) ::  kk                              !< loop index
     872    INTEGER(iwp) ::  r                               !< loop index
    868873    INTEGER(iwp) ::  mid                             !< masked output running index
    869874    INTEGER(iwp) ::  netcdf_data_format_save         !< initial value of netcdf_data_format
     
    36473652          CALL message( 'check_parameters', 'PA0421', 1, 2, 0, 6, 0 )
    36483653       ENDIF
     3654
     3655!
     3656!--    Convert recycle_absolute_quantities (list of strings that define the quantities for
     3657!--    absolute recycling) to raq (list of logicals with length 7 corresponding to u,v,w,pt,e,q,s).
     3658!--    Output error message for not implemented quantities.
     3659       DO r = LBOUND( recycle_absolute_quantities, 1 ), UBOUND( recycle_absolute_quantities, 1 )
     3660          SELECT CASE ( TRIM( recycle_absolute_quantities(r) ) )
     3661             CASE ( 'theta' )
     3662                raq(4) = .TRUE.
     3663             CASE ( 'q' )
     3664                raq(6) = .TRUE.
     3665             CASE ( '' )
     3666                CONTINUE
     3667             CASE DEFAULT
     3668                message_string = 'absolute recycling not implemented for variable ' // &
     3669                TRIM( recycle_absolute_quantities(r) )
     3670                CALL message( 'inflow_turbulence', 'PA184', 1, 2, 0, 6, 0 )
     3671          END SELECT
     3672       ENDDO
    36493673    ENDIF
    36503674
Note: See TracChangeset for help on using the changeset viewer.