Ignore:
Timestamp:
Feb 14, 2020 1:39:31 PM (5 years ago)
Author:
Giersch
Message:

Tutorials updated, all tested with revision 4402

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/TUTORIALS/cases/lsm_long/INPUT/lsm_long_p3d

    r4397 r4409  
    77!-- grid parameters
    88!-------------------------------------------------------------------------------
    9     nx                         = 19, ! Number of gridboxes in x-direction (nx+1)
    10     ny                         = 19, ! Number of gridboxes in y-direction (ny+1)
    11     nz                         = 50, ! Number of gridboxes in z-direction (nz)
     9    nx                         = 19,
     10    ny                         = 19,
     11    nz                         = 50,
    1212
    13     dx                         = 100.0, ! Size of single gridbox in x-direction
    14     dy                         = 100.0, ! Size of single gridbox in y-direction
    15     dz                         = 25.0, ! Size of single gridbox in z-direction
     13    dx                         = 100.0,
     14    dy                         = 100.0,
     15    dz                         = 25.0,
    1616
    17     dz_stretch_level           = 100.0, ! Height (in m) where stretching starts
    18     dz_max                     = 50.0,  ! maximum vertical grid size
     17    dz_stretch_level           = 100.0,
     18    dz_max                     = 50.0,
    1919!
    2020!-- initialization
    2121!------------------------------------------------------------------------------
    22     initializing_actions       = 'set_constant_profiles', ! initial conditions
     22    initializing_actions       = 'set_constant_profiles',
    2323
    24     ug_surface                 = 2.5, ! u-comp of geostrophic wind at surface
    25     vg_surface                 = 0.0, ! v-comp of geostrophic wind at surface
     24    ug_surface                 = 2.5,
     25    vg_surface                 = 0.0,
    2626
    27     pt_surface                 = 293.0, ! initial surface potential temp
     27    pt_surface                 = 293.0,
    2828
    29     pt_vertical_gradient       = 0.0,
    30                                  3.0,
    31                                  0.0,    ! piecewise temp gradients
    32     pt_vertical_gradient_level = 0.0,
    33                                  800.0,
    34                                  1600.0, ! height level of temp gradients
     29    pt_vertical_gradient       = 0.0, 3.0, 0.0,
     30    pt_vertical_gradient_level = 0.0, 800.0, 1600.0,
    3531 
    3632    q_surface                  = 0.01,   ! mixing ratio --> rel. hum. = ?
     
    3834    surface_pressure           = 1000.0, ! in hPa
    3935
    40     origin_date_time           = '2019-06-21 12:00:00 +00'
     36    origin_date_time           = '2019-06-21 12:00:00 +00', ! does the year matter?
    4137!
    4238!-- physics
    4339!------------------------------------------------------------------------------
    44     latitude                   = 53.37, ! geographical latitude
    45     longitude                  = 9.73,  ! geographical longitude
     40    latitude                   = 53.37,
     41    longitude                  = 9.73,
    4642!
    4743!-- boundary conditions
     
    6460    fft_method                 = 'temperton-algorithm', ! build-in fft method
    6561
    66     rayleigh_damping_height    = 1600.0, ! starting height for Rayleigh damping
    67     rayleigh_damping_factor    = 0.01,   ! strength of Rayleigh damping
     62    rayleigh_damping_height    = 1600.0,
     63    rayleigh_damping_factor    = 0.01,
    6864
    6965/
     
    7773!-- run steering
    7874!-------------------------------------------------------------------------------
    79     end_time                   = 216000.0, ! simulation time of the 3D model
    80                                            ! (2.5 days)
     75    end_time                   = 216000.0, ! (2.5 days)
    8176
    82     create_disturbances        = .T.,    ! randomly perturbate horiz. velocity
    83     dt_disturb                 = 120.0,  ! interval for random perturbations
    84     disturbance_energy_limit   = 0.0001, ! upper limit for perturbation energy
     77    create_disturbances        = .TRUE.,
     78    dt_disturb                 = 120.0,
     79    disturbance_energy_limit   = 0.0001,
    8580!
    8681!-- data output
     
    8883    netcdf_data_format         = 2, ! use NetCDF3
    8984
    90     dt_run_control             = 900.0,   ! output interval for run control
    91     dt_data_output             = 21600.0, ! output interval for general data
    92     dt_data_output_av          = 3600.0,  ! output interval for averaged data
     85    dt_run_control             = 900.0,
     86    dt_data_output             = 21600.0,
     87    dt_data_output_av          = 3600.0,
    9388
    94     averaging_interval         = 1800.0,  ! averaging interval general data
     89    averaging_interval         = 1800.0,
    9590
    96     data_output_pr             = '#u', '#v', '#theta', '#q', '#thetav',
    97                                  '#km', '#kh', '#l',
    98                                  'w', 'e', 'e*', 'p',
    99                                  'w"u"', 'w*u*', 'w"v"', 'w*v*',
    100                                  'w"theta"', 'w*theta*', 'w"thetav"',
    101                                  'w*thetav*',
    102                                  'w"q"', 'w*q*', 'w*e*',
    103                                  'theta*2', 'q*2', 'u*2', 'v*2', 'w*2',
    104                                  '#t_soil', '#m_soil',  ! Profile output
     91    data_output_pr             = '#u', '#v', '#theta', '#q', '#thetav', '#km', '#kh', '#l', 'w', 'e', 'e*', 'p', 'w"u"', 'w*u*', 'w"v"', 'w*v*', 'w"theta"', 'w*theta*', 'w"thetav"', 'w*thetav*', 'w"q"', 'w*q*', 'w*e*', 'theta*2', 'q*2', 'u*2', 'v*2', 'w*2', '#t_soil', '#m_soil',
    10592
    106     section_xy                 = 0, ! surface variables only
     93    section_xy                 = 0,
    10794
    108     data_output                = 'theta',         'theta_av',
    109                                  'thetav',        'thetav_av',
    110                                  'q',             'q_av',
    111                                  'u',             'u_av',
    112                                  'v',             'v_av',
     95    data_output                = 'theta', 'theta_av', 'thetav', 'thetav_av', 'q', 'q_av', 'u', 'u_av', 'v', 'v_av', 't_soil', 't_soil_av', 'm_soil', 'm_soil_av', 'us*_xy', 'us*_xy_av', 't*_xy', 't*_xy_av', 'r_a*_xy', 'r_a*_xy_av', 'r_s*_xy', 'r_s*_xy_av', 'tsurf*_xy', 'tsurf*_xy_av', 'ghf*_xy', 'ghf*_xy_av', 'shf*_xy', 'shf*_xy_av', 'qsws_liq*_xy', 'qsws_liq*_xy_av', 'qsws_soil*_xy', 'qsws_soil*_xy_av', 'qsws_veg*_xy', 'qsws_veg*_xy_av', 'c_liq*_xy', 'c_liq*_xy_av', 'c_soil*_xy', 'c_soil*_xy_av', 'c_veg*_xy', 'c_veg*_xy_av',
    11396
    114                                  't_soil',        't_soil_av',
    115                                  'm_soil',        'm_soil_av',
    116 
    117                                  'us*_xy',        'us*_xy_av',
    118                                  't*_xy',         't*_xy_av',
    119                                  'r_a*_xy',       'r_a*_xy_av',
    120                                  'r_s*_xy',       'r_s*_xy_av',
    121                                  'tsurf*_xy',     'tsurf*_xy_av'
    122 
    123                                  'ghf*_xy',       'ghf*_xy_av',
    124                                  'shf*_xy',       'shf*_xy_av',
    125 
    126                                  'qsws_liq*_xy',  'qsws_liq*_xy_av',
    127                                  'qsws_soil*_xy', 'qsws_soil*_xy_av',
    128                                  'qsws_veg*_xy',  'qsws_veg*_xy_av',
    129 
    130                                  'c_liq*_xy',     'c_liq*_xy_av',
    131                                  'c_soil*_xy',    'c_soil*_xy_av',
    132                                  'c_veg*_xy',     'c_veg*_xy_av', ! 2d and/or
    133                                                                   ! 3d output
    13497/
    13598
     
    139102!-------------------------------------------------------------------------------
    140103&land_surface_parameters
    141     soil_type                  = 3,  ! type to be used in soil model
    142     soil_moisture              = 0.35,  0.35,  0.35,  0.35, ! moisture at every
    143                                  0.35,  0.35,  0.35,  0.35, ! soil level
    144     soil_temperature           = 293.0, 293.0, 293.0, 293.0, ! temperature at
    145                                  293.0, 293.0, 293.0, 293.0, ! every soil level
    146     deep_soil_temperature      = 293.0,  ! temperature of deep soil
     104    soil_type                  = 3,  ! medium-fine
     105    soil_moisture              = 0.35,  0.35,  0.35,  0.35, 0.35,  0.35,  0.35,  0.35,
     106    soil_temperature           = 293.0, 293.0, 293.0, 293.0, 293.0, 293.0, 293.0, 293.0,
     107    deep_soil_temperature      = 293.0,
    147108
    148     surface_type               = 'vegetation',  ! surface classification
    149     vegetation_type            = 3, ! vegetation to be used in surface model
    150 
    151     vegetation_coverage        = 0.95, ! part of surface covered by vegetation
    152                                        ! (0 - 1)
    153     c_surface                  = 0.0,  ! surface heat capacity
     109    surface_type               = 'vegetation',
     110    vegetation_type            = 3, ! short grass
     111    vegetation_coverage        = 0.95,
    154112    conserve_water_content     = .TRUE.,  ! closed (T) or open (F) lower boundary
    155113
     
    162120&radiation_parameters
    163121    radiation_scheme           = 'clear-sky', ! Radiation scheme to be used
    164     dt_radiation               = 60.0, ! Time step of the radiation model
     122    dt_radiation               = 60.0,
    165123
    166124/
Note: See TracChangeset for help on using the changeset viewer.