Ignore:
Timestamp:
Sep 9, 2019 12:55:23 PM (2 years ago)
Author:
gronemeier
Message:

update tutorials flow_around_cube and lsm

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/TUTORIALS/cases/flow_around_cube_noncyclic/INPUT/flow_around_cube_noncyclic_p3d

    r4098 r4222  
     1!-------------------------------------------------------------------------------
     2!-- INITIALIZATION PARAMETER NAMELIST
     3!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar
     4!-------------------------------------------------------------------------------
    15&initialization_parameters
    2    !-- grid setup
    3       nx = 79,
    4       ny = 39,
    5       nz = 40,
     6!
     7!-- grid parameters
     8!-------------------------------------------------------------------------------
     9    nx                         = 79, ! Number of gridboxes in x-direction (nx+1)
     10    ny                         = 39, ! Number of gridboxes in y-direction (ny+1)
     11    nz                         = 40, ! Number of gridboxes in z-direction (nz)
    612
    7       dx = 2.0,
    8       dy = 2.0,
    9       dz = 2.0,
     13    dx                         = 2.0, ! Size of single gridbox in x-direction
     14    dy                         = 2.0, ! Size of single gridbox in y-direction
     15    dz                         = 2.0, ! Size of single gridbox in z-direction
     16!
     17!-- initialization
     18!-------------------------------------------------------------------------------
     19    initializing_actions       = 'cyclic_fill', ! initial conditions
    1020
    11    !-- pressure solver
    12       psolver   = 'multigrid',                ! required for non-cyclic bc
    13       mg_cycles = 2,                          ! needs to be checked
     21    ug_surface                 = 1.0, ! u-comp of geostrophic wind at surface
     22    vg_surface                 = 0.0, ! v-comp of geostrophic wind at surface
    1423
    15    !-- init-action
    16       initializing_actions = 'cyclic_fill',
     24    neutral                    = .T., ! strictly neutral flow
     25!
     26!-- physics
     27!-------------------------------------------------------------------------------
     28    omega                      = 0.0, ! no Coriolis force
     29!
     30!-- boundary conditions
     31!-------------------------------------------------------------------------------
     32    bc_uv_t                    = 'neumann', ! free-slip boundary condition
     33    bc_lr                      = 'dirichlet/radiation',  ! non-cyclic
     34                                                         ! left/right boundary
     35    turbulent_inflow           = .TRUE., ! use turbulent inflow method
     36    recycling_width            = 80.0,   ! Distance between recycling plane and
     37                                         ! inflow boundary
     38    inflow_damping_height      = 80.0, ! height limitation for turbulent inflow
     39!
     40!-- numerics
     41!-------------------------------------------------------------------------------
     42    psolver                    = 'multigrid', ! set pressure solver
     43    mg_cycles                  = 2,  ! set number of cycles for multigrid solver
     44!
     45!-- topography
     46!-------------------------------------------------------------------------------
     47    topography                 = 'read_from_file',  ! use input file
     48                                                    ! to set topography
    1749
    18    !-- advection schemes
    19       momentum_advec = 'ws-scheme',           ! default advection scheme
    20       scalar_advec   = 'ws-scheme',
     50/ ! end of initialization parameter namelist
    2151
    22    !-- wind speed setup
    23       ug_surface = 1.0,
    24       vg_surface = 0.0,
     52!-------------------------------------------------------------------------------
     53!-- RUNTIME PARAMETER NAMELIST
     54!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par
     55!-------------------------------------------------------------------------------
     56&runtime_parameters
     57!
     58!-- run steering
     59!-------------------------------------------------------------------------------
     60    end_time                   = 3600.0,  ! simulation time of the 3D model
     61    create_disturbances        = .FALSE., ! randomly perturbate horiz. velocity
     62!
     63!-- data output
     64!-------------------------------------------------------------------------------
     65    dt_run_control             =   60.0, ! output interval for run control
     66    dt_data_output             = 1800.0, ! output interval for general data
     67    dt_dots                    =   60.0, ! output interval for time-series data
     68    dt_dopr                    = 1800.0, ! output interval for profile data
    2569
    26    !-- temperature setup
    27       neutral = .T.,                          ! strictly neutral flow
     70    data_output                = 'u', 'u_av',
     71                                 'v', 'v_av',
     72                                 'w', 'w_av',  ! 2d and/
     73                                               ! or 3d
     74                                               ! output
    2875
    29    !-- boundary conditions
    30       bc_uv_t = 'neumann',                    ! free-slip boundary condition
    31       bc_lr   = 'dirichlet/radiation',        ! non-cyclic left/right boundary
     76    data_output_pr             = '#u', 'u*2', 'wu', 'w*u*', 'w"u"',
     77                                 '#v', 'v*2', 'wv', 'w*v*', 'w"v"',
     78                                 'w', 'w*2',
     79                                 'e', 'e*',
     80                                 '#km',
     81                                 '#l',      ! Profile
     82                                            ! output
    3283
    33       roughness_length = 0.1,                 ! applies to all surfaces locally
     84    averaging_interval         = 1800.0, ! averaging interval general data
     85    dt_averaging_input         =    6.0, ! averaging general data sampling rate
    3486
    35    !-- special methods
    36       turbulent_inflow      = .TRUE.,
    37       recycling_width       = 80.0,
    38       inflow_damping_height = 80.0,
     87    averaging_interval_pr      = 1800.0, ! averaging interval profile data
     88    dt_averaging_input_pr      =    6.0, ! averaging profile data sampling rate
    3989
    40    !-- special setups
    41       omega = 0.0,                            ! no Coriolis force
    42 
    43       topography = 'read_from_file',
    44       topography_grid_convention = 'cell_center',   ! default
    45    /
    46 
    47 
    48 &runtime_parameters
    49    !-- simulation time
    50       end_time                = 3600.0,
    51       termination_time_needed =  900.0,
    52 
    53    !-- others
    54 
    55    !-- global output settings
    56       skip_time_data_output =    0.0,
    57       dt_data_output        = 1800.0,
    58 
    59       averaging_interval    = 1800.0,
    60       dt_averaging_input    =    6.0,
    61 
    62    !-- time series / run control
    63       dt_run_control        =  60.0,
    64       dt_dots               =  60.0,
    65 
    66    !-- profiles
    67       skip_time_dopr        =    0.0,
    68       dt_dopr               = 1800.0,
    69       averaging_interval_pr = 1800.0,
    70       dt_averaging_input_pr =    6.0,
    71 
    72       data_output_pr = '#u', 'u*2', 'wu', 'w*u*', 'w"u"',
    73                        '#v', 'v*2', 'wv', 'w*v*', 'w"v"',
    74                        'w', 'w*2',
    75                        'e', 'e*',
    76                        '#km',
    77                        '#l',
    78 
    79    !-- output variables
    80       data_output = 'u', 'u_av',
    81                     'v', 'v_av',
    82                     'w', 'w_av',
    83    /
     90/ ! end of runtime parameter namelist
Note: See TracChangeset for help on using the changeset viewer.