Ignore:
Timestamp:
Jun 13, 2016 7:12:51 AM (5 years ago)
Author:
hellstea
Message:

last commit documented

File:
1 edited

Legend:

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

    r1917 r1933  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! 
    2222!
    2323! Former revisions:
     
    2525! $Id$
    2626!
     27! 1917 2016-05-27 14:28:12Z witha
     28! Initial version of purely vertical nesting introduced.
     29!
    2730! 1914 2016-05-26 14:44:07Z witha
    2831! Added call to wind turbine model for reading of &wind_turbine_par
     
    245248
    246249    USE pmc_interface,                                                         &
    247         ONLY:  nested_run
     250        ONLY:  nested_run, nesting_mode
    248251
    249252    USE profil_parameter,                                                      &
     
    429432
    430433!
    431 !--       In case of nested runs, explicitly set nesting boundary conditions
    432 !--       except for the root domain. This will overwrite the user settings.
    433           IF ( nest_domain )  THEN
    434              bc_lr   = 'nested'
    435              bc_ns   = 'nested'
    436              bc_uv_t = 'nested'
    437              bc_pt_t = 'nested'
    438              bc_q_t  = 'nested'
    439              bc_p_t  = 'neumann'
     434!--       In case of nested runs, explicitly set nesting boundary conditions.
     435!--       This will overwrite the user settings. bc_lr and bc_ns always need
     436!--       to be cyclic for vertical nesting.
     437          IF ( nesting_mode == 'vertical' )  THEN
     438             IF (bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic' )  THEN
     439                WRITE ( message_string, *)  'bc_lr and bc_ns were set to ,',   &
     440                                      'cyclic for vertical nesting'
     441                CALL message( 'parin', 'PA0428', 0, 0, 0, 6, 0 )
     442                bc_lr   = 'cyclic'
     443                bc_ns   = 'cyclic'
     444             ENDIF
     445             IF ( nest_domain )  THEN
     446                bc_uv_t = 'nested'
     447                bc_pt_t = 'nested'
     448                bc_q_t  = 'nested'
     449                bc_p_t  = 'neumann'
     450             ENDIF
     451          ELSE
     452         
     453!
     454!--       For other nesting modes only set boundary conditions for
     455!--       nested domains.
     456             IF ( nest_domain )  THEN
     457                bc_lr   = 'nested'
     458                bc_ns   = 'nested'
     459                bc_uv_t = 'nested'
     460                bc_pt_t = 'nested'
     461                bc_q_t  = 'nested'
     462                bc_p_t  = 'neumann'
     463             ENDIF
    440464          ENDIF
     465           
    441466!
    442467!--       Check validity of lateral boundary conditions. This has to be done
Note: See TracChangeset for help on using the changeset viewer.