Ignore:
Timestamp:
Aug 23, 2019 7:33:16 AM (5 years ago)
Author:
oliver.maas
Message:

simplified steering of recycling of absolute values. Replaced initialization parameter recycle_absolute_quantities by recycling_method_for_thermodynamic_quantities.

File:
1 edited

Legend:

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

    r4182 r4183  
    2525! -----------------
    2626! 4172 2019-08-20 11:55:33Z oliver.maas
     27! removed conversion from recycle_absolute_quantities to raq, added check and
     28! error message for correct input of recycling_method_for_thermodynamic_quantities
     29!
     30! 11:55:33Z oliver.maas
    2731! Corrected "Former revisions" section
    2832!
     
    164168    INTEGER(iwp) ::  k                               !< loop index
    165169    INTEGER(iwp) ::  kk                              !< loop index
    166     INTEGER(iwp) ::  r                               !< loop index
    167170    INTEGER(iwp) ::  mid                             !< masked output running index
    168171    INTEGER(iwp) ::  netcdf_data_format_save         !< initial value of netcdf_data_format
     
    29282931
    29292932!
    2930 !-- In case of turbulent inflow calculate the index of the recycling plane
     2933!-- In case of turbulent inflow
    29312934    IF ( turbulent_inflow )  THEN
     2935
     2936!
     2937!--    Calculate the index of the recycling plane
    29322938       IF ( recycling_width <= dx  .OR.  recycling_width >= nx * dx )  THEN
    29332939          WRITE( message_string, * )  'illegal value for recycling_width: ',   &
     
    29482954
    29492955!
    2950 !--    Convert recycle_absolute_quantities (list of strings that define the quantities for
    2951 !--    absolute recycling) to raq (list of logicals with length 7 corresponding to u,v,w,pt,e,q,s).
    2952 !--    Output error message for not implemented quantities.
    2953        DO r = LBOUND( recycle_absolute_quantities, 1 ), UBOUND( recycle_absolute_quantities, 1 )
    2954           SELECT CASE ( TRIM( recycle_absolute_quantities(r) ) )
    2955              CASE ( 'theta' )
    2956                 raq(4) = .TRUE.
    2957              CASE ( 'q' )
    2958                 raq(6) = .TRUE.
    2959              CASE ( '' )
    2960                 CONTINUE
    2961              CASE DEFAULT
    2962                 message_string = 'absolute recycling not implemented for variable ' // &
    2963                 TRIM( recycle_absolute_quantities(r) )
    2964                 CALL message( 'inflow_turbulence', 'PA0184', 1, 2, 0, 6, 0 )
    2965           END SELECT
    2966        ENDDO
     2956!--   Check for correct input of recycling method for thermodynamic quantities
     2957       IF ( TRIM( recycling_method_for_thermodynamic_quantities ) /= 'turbulent_fluctuation' .AND. &
     2958            TRIM( recycling_method_for_thermodynamic_quantities ) /= 'absolute_value' )  THEN
     2959          WRITE( message_string, * )  'unknown recycling method for thermodynamic quantities: ',   &
     2960               TRIM( recycling_method_for_thermodynamic_quantities )
     2961          CALL message( 'check_parameters', 'PA0184', 1, 2, 0, 6, 0 )
     2962       ENDIF
     2963
    29672964    ENDIF
    29682965
Note: See TracChangeset for help on using the changeset viewer.