Changeset 2826 for palm/trunk


Ignore:
Timestamp:
Feb 21, 2018 12:39:28 PM (4 years ago)
Author:
hellstea
Message:

Bugfix in setting the default boundary conditions for nest domains

File:
1 edited

Legend:

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

    r2817 r2826  
    2525! -----------------
    2626! $Id$
     27! Bugfix in setting the default boundary conditions for nest domains
     28!
     29! 2817 2018-02-19 16:32:21Z knoop
    2730! Preliminary gust module interface implemented
    2831!
     
    615618!
    616619!--       In case of nested runs, explicitly set nesting boundary conditions.
    617 !--       This will overwrite the user settings. bc_lr and bc_ns always need
    618 !--       to be cyclic for vertical nesting.
    619           IF ( nesting_mode == 'vertical' )  THEN
    620              IF (bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic' )  THEN
    621                 WRITE ( message_string, *)  'bc_lr and bc_ns were set to ,',   &
    622                                       'cyclic for vertical nesting'
    623                 CALL message( 'parin', 'PA0428', 0, 0, 0, 6, 0 )
    624                 bc_lr   = 'cyclic'
    625                 bc_ns   = 'cyclic'
    626              ENDIF
    627              IF ( nest_domain )  THEN
    628                 bc_uv_t  = 'nested'
    629                 bc_pt_t  = 'nested'
    630                 bc_q_t   = 'nested'
    631                 bc_s_t   = 'nested'
    632                 bc_cs_t  = 'nested'
    633                 bc_p_t   = 'neumann' 
    634              ENDIF       
    635 !
    636 !--       For other nesting modes only set boundary conditions for
    637 !--       nested domains.
    638              IF ( nest_domain )  THEN
    639                 bc_lr    = 'nested'
    640                 bc_ns    = 'nested'
    641                 bc_uv_t  = 'nested'
    642                 bc_pt_t  = 'nested'
    643                 bc_q_t   = 'nested'
    644                 bc_s_t   = 'nested'
    645                 bc_cs_t = 'nested'
    646                 bc_p_t   = 'neumann'
     620!--       This will overwrite the user settings and basic defaults.
     621!--       bc_lr and bc_ns always need to be cyclic for vertical nesting.
     622          IF ( nested_run )  THEN
     623             IF ( nesting_mode == 'vertical' )  THEN
     624                IF (bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic' )  THEN
     625                   WRITE ( message_string, *) 'bc_lr and bc_ns were set to ,', &
     626                        'cyclic for vertical nesting'
     627                   CALL message( 'parin', 'PA0428', 0, 0, 0, 6, 0 )
     628                   bc_lr   = 'cyclic'
     629                   bc_ns   = 'cyclic'
     630                ENDIF
     631                IF ( nest_domain )  THEN
     632                   bc_uv_t  = 'nested'
     633                   bc_pt_t  = 'nested'
     634                   bc_q_t   = 'nested'
     635                   bc_s_t   = 'nested'
     636                   bc_cs_t  = 'nested'
     637                   bc_p_t   = 'neumann' 
     638                ENDIF
     639!
     640!--          For other nesting modes only set boundary conditions for
     641!--          nested domains.
     642             ELSE
     643                IF ( nest_domain )  THEN
     644                   bc_lr    = 'nested'
     645                   bc_ns    = 'nested'
     646                   bc_uv_t  = 'nested'
     647                   bc_pt_t  = 'nested'
     648                   bc_q_t   = 'nested'
     649                   bc_s_t   = 'nested'
     650                   bc_cs_t  = 'nested'
     651                   bc_p_t   = 'neumann'
     652                ENDIF
    647653             ENDIF
    648654          ENDIF
Note: See TracChangeset for help on using the changeset viewer.