Ignore:
Timestamp:
Mar 27, 2018 3:52:42 PM (3 years ago)
Author:
suehring
Message:

Nesting in RANS-LES and RANS-RANS mode enabled; synthetic turbulence generator at all lateral boundaries in nesting or non-cyclic forcing mode; revised Inifor initialization in nesting mode

File:
1 edited

Legend:

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

    r2934 r2938  
    2525! -----------------
    2626! $Id$
     27! - Revise start and end indices for imposing random disturbances in case of
     28!   nesting or large-scale forcing.
     29! - Remove check for inifor initialization in case of nested runs.
     30! - Adapt call to check for synthetic turbulence geneartor settings.
     31!
     32! 2936 2018-03-27 14:49:27Z suehring
    2733! Check spinup in case of nested runs, spinup time and timestep must be
    2834! identical to assure synchronuous simulations.
     
    13821388    ENDIF
    13831389!
    1384 !-- In case of nested run assure that all domains are initialized the same
    1385 !-- way, i.e. if at least at one domain is initialized with soil and
    1386 !-- atmospheric data provided by COSMO, all domains must be initialized the
    1387 !-- same way, to assure that soil and atmospheric quantities are
    1388 !-- consistent.
    1389     IF ( nested_run )  THEN
    1390        check_nest = .TRUE.
    1391 #if defined( __parallel )
    1392        CALL MPI_ALLREDUCE( TRIM( initializing_actions ) == 'inifor',           &
    1393                            check_nest, 1, MPI_LOGICAL,                         &
    1394                            MPI_LAND, MPI_COMM_WORLD, ierr )
    1395 
    1396        IF ( TRIM( initializing_actions ) == 'inifor'  .AND.                    &
    1397             .NOT.  check_nest )  THEN
    1398           message_string = 'In case of nesting, if at least in one ' //        &
    1399                            'domain initializing_actions = inifor, '  //        &
    1400                            'all domains need to be initialized that way.'
    1401           CALL message( 'netcdf_data_input_mod', 'PA0430', 3, 2, 0, 6, 0 )
    1402        ENDIF
    1403 #endif
    1404     ENDIF
    1405 !
    14061390!-- In case of spinup and nested run, spinup end time must be identical
    14071391!-- in order to have synchronously running simulations.
     
    14361420
    14371421!
    1438 !-- When synthetic turbulence generator is used, peform addtional checks
    1439     IF ( syn_turb_gen )  CALL stg_check_parameters
     1422!-- Check for synthetic turbulence generator settings
     1423    CALL stg_check_parameters
    14401424!
    14411425!-- When plant canopy model is used, peform addtional checks
     
    38983882          dist_nxr(1) = MIN( inflow_disturbance_end, nxr )
    38993883       ELSEIF ( bc_lr == 'nested'  .OR.  bc_lr == 'forcing' )  THEN
    3900           dist_nxl = MAX( inflow_disturbance_begin, nxl )
     3884          dist_nxl    = MAX( inflow_disturbance_begin, nxl )
     3885          dist_nxr    = MIN( nx - inflow_disturbance_begin, nxr )
    39013886       ENDIF
    39023887       IF ( bc_ns == 'dirichlet/radiation' )  THEN
     
    39073892          dist_nyn(1) = MIN( inflow_disturbance_end, nyn )
    39083893       ELSEIF ( bc_ns == 'nested'  .OR.  bc_ns == 'forcing' )  THEN
    3909           dist_nys = MAX( inflow_disturbance_begin, nys )
     3894          dist_nys    = MAX( inflow_disturbance_begin, nys )
     3895          dist_nyn    = MIN( ny - inflow_disturbance_begin, nyn )
    39103896       ENDIF
    39113897    ELSE
     
    39183904          dist_nxl    = inflow_disturbance_begin
    39193905          dist_nxr(1) = inflow_disturbance_end
     3906       ELSEIF ( bc_lr == 'nested'  .OR.  bc_lr == 'forcing' )  THEN
     3907          dist_nxr    = nx - inflow_disturbance_begin
     3908          dist_nxl    = inflow_disturbance_begin
    39203909       ENDIF
    39213910       IF ( bc_ns == 'dirichlet/radiation' )  THEN
     
    39253914          dist_nys    = inflow_disturbance_begin
    39263915          dist_nyn(1) = inflow_disturbance_end
     3916       ELSEIF ( bc_ns == 'nested'  .OR.  bc_ns == 'forcing' )  THEN
     3917          dist_nyn    = ny - inflow_disturbance_begin
     3918          dist_nys    = inflow_disturbance_begin
    39273919       ENDIF
    39283920    ENDIF
Note: See TracChangeset for help on using the changeset viewer.