Ignore:
Timestamp:
Sep 10, 2019 12:49:41 PM (3 years ago)
Author:
Giersch
Message:

Three Tutorial cases updated for newest version

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/TUTORIALS/cases/example_neutral/INPUT/example_neutral_p3d

    r4098 r4224  
    1  &initialization_parameters
     1!-------------------------------------------------------------------------------
     2!-- INITIALIZATION PARAMETER NAMELIST
     3!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar
     4!-------------------------------------------------------------------------------
     5&initialization_parameters
     6!
     7!-- grid parameters
     8!-------------------------------------------------------------------------------
     9    nx                         = 49, ! Number of gridboxes in x-direction (nx+1)
     10    ny                         = 49, ! Number of gridboxes in y-direction (ny+1)
     11    nz                         = 100,! Number of gridboxes in z-direction (nz)
    212
    3      !-- Numerical grid   
    4          nx = 49,
    5          ny = 49,
    6          nz = 100,
     13    dx                         = 20.0, ! Size of single gridbox in x-direction
     14    dy                         = 20.0, ! Size of single gridbox in y-direction
     15    dz                         = 20.0, ! Size of single gridbox in z-direction
     16!
     17!-- initialization
     18!-------------------------------------------------------------------------------
     19    initializing_actions       = 'set_constant_profiles' , ! initial conditions
    720
    8          dx = 20.0,
    9          dy = 20.0,
    10          dz = 20.0,
     21    ug_surface                 = 5.0 , ! u-comp of geostrophic wind at surface
     22    vg_surface                 = 0.0 , ! v-comp of geostrophic wind at surface
    1123
    12      !-- Pressure solver 
    13          fft_method = 'temperton-algorithm',
     24    pt_surface                 = 300.0,     ! initial surface potential temp
     25    pt_vertical_gradient       = 0.0,       ! piecewise temp gradients
     26    pt_vertical_gradient_level = 0.0, ! height level of temp gradients
     27!
     28!-- boundary conditions
     29!-------------------------------------------------------------------------------
     30    bc_uv_b                    =  'dirichlet' ,
     31    bc_uv_t                    =   'neumann',
    1432
    15      !-- Initialization
    16 !         initializing_actions = 'set_constant_profiles',
     33    bc_p_b                     = 'neumann',
     34    bc_p_t                     = 'dirichlet',
    1735
    18          initializing_actions = 'set_1d-model_profiles',
    19          end_time_1d = 864000.0,
    20          damp_level_1d = 1250.0,
    21          mixing_length_1d = 'blackadar',
    22          dissipation_1d = 'detering',
     36    constant_flux_layer        = .T. ,
     37   
     38    roughness_length           = 0.1 ,
     39!
     40!-- numerics
     41!-------------------------------------------------------------------------------
     42    fft_method                 = 'temperton-algorithm',   ! build-in fft method     
    2343
    24          dt_run_control_1d = 3600.0,
    25          dt_pr_1d = 86400.0,
     44    rayleigh_damping_height    = 1250.0,
     45    rayleigh_damping_factor    = 0.01,
     46     
     47    galilei_transformation     = .T. ,
     48   
     49/ ! end of initialization parameter namelist
    2650
    27      !-- Rayleigh damping
    28          rayleigh_damping_height = 1250.0,
    29          rayleigh_damping_factor = 0.01,
    30          galilei_transformation = .T.,
     51!-------------------------------------------------------------------------------
     52!-- RUNTIME PARAMETER NAMELIST
     53!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par
     54!-------------------------------------------------------------------------------
     55&runtime_parameters 
     56!
     57!-- run steering
     58!-------------------------------------------------------------------------------
     59    end_time                   = 172800.0, ! simulation time of the 3D model
    3160
    32      !-- Roughness length         
    33          roughness_length = 0.1,
     61    create_disturbances        = .T.,      ! randomly perturbate horiz. velocity
     62    dt_disturb                 = 150.0,             ! interval for random perturbations
     63    disturbance_energy_limit   = 0.01,! upper limit for perturbation energy
     64!
     65!-- data output
     66!-------------------------------------------------------------------------------
     67    netcdf_data_format         = 2, ! use NetCDF3
     68    dt_run_control             = 60.0,
    3469
    35      !-- Geostrophic wind (large-scale pressure gradient)
    36          ug_surface = 5.0,
    37          vg_surface = 0.0,
     70    dt_dopr                    = 10800.0,
     71         
     72    averaging_interval_pr      = 10800.0,
     73    dt_averaging_input_pr      = 50.0,
     74    data_output_pr             = '#u', '#v', 'w', '#theta',
     75                                 'w"u"', 'w*u*', 'wu', 'w"v"', 'w*v*', 'wv',
     76                                 'w*2', 'theta*2',
     77                       
     78/! end of runtime parameter namelist
    3879
    39      !-- Neutral stratification
    40          pt_surface = 300.0,
    41          pt_vertical_gradient = 0.0,
    42          pt_vertical_gradient_level = 0.0,
    43          neutral = .T.,
    44 
    45      !-- Boundary conditions
    46          bc_uv_b =  'dirichlet',
    47          bc_uv_t = 'neumann',
    48 
    49          bc_p_b = 'neumann',
    50          bc_p_t = 'dirichlet',
    51 
    52          constant_flux_layer = .T.,
    53          /
    54 
    55  &runtime_parameters 
    56 
    57      !-- Simulation time 2 days
    58          end_time = 172800.0,
    59 
    60      !-- Disturbances
    61          create_disturbances = .T.,
    62          dt_disturb = 150.0,
    63          disturbance_energy_limit = 0.01,
     80!-------------------------------------------------------------------------------
     81!-- SPECTRA PARAMETER NAMELIST
     82!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/sppar
     83!-------------------------------------------------------------------------------         
     84&spectra_parameters
    6485 
    65      !-- Run control
    66          dt_run_control = 60.0,
    67 
    68      !-- Time series
    69          dt_dopr = 10800.0,
     86    data_output_sp            =  'u', 'v', 'w',
     87    dt_dosp                   = 10800.0,
     88    comp_spectra_level        =  2, 5, 10, 15, 20,
     89    spectra_direction         =  'x' , 'y' , 'x',
     90    averaging_interval_sp    = 10800.0,
    7091         
    71      !-- Vertical profiles
    72          averaging_interval_pr = 10800.0,
    73          dt_averaging_input_pr = 50.0,
    74          data_output_pr = '#u', '#v', 'w', '#theta',
    75                           'w"u"', 'w*u*', 'wu', 'w"v"', 'w*v*', 'wv','w*2', 'theta*2',
    76                          
    77          /
    78          
    79  &spectra_parameters
     92/! end of spectra parameter namelist
    8093 
    81      !-- Output settings for spectra
    82          data_output_sp =  'u', 'v', 'w',
    83          dt_dosp = 10800.0,
    84          comp_spectra_level =  2, 5, 10, 15, 20,
    85          spectra_direction =  'x' , 'y' , 'x',
    86          averaging_interval_sp = 10800.0,
    87          
    88          /
    89  
Note: See TracChangeset for help on using the changeset viewer.