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

update tutorials flow_around_cube and lsm

File:
1 edited

Legend:

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

    r4098 r4222  
    22!
    33!-- grid parameters
    4 !------------------------------------------------------------------------------
    5    nx                         = 19,
    6    ny                         = 19,
    7    nz                         = 50,
     4!-------------------------------------------------------------------------------
     5    nx                         = 19, ! Number of gridboxes in x-direction (nx+1)
     6    ny                         = 19, ! Number of gridboxes in y-direction (ny+1)
     7    nz                         = 50, ! Number of gridboxes in z-direction (nz)
    88
    9    dx                         = 100.0,
    10    dy                         = 100.0,
    11    dz                         = 25.0,
     9    dx                         = 100.0, ! Size of single gridbox in x-direction
     10    dy                         = 100.0, ! Size of single gridbox in y-direction
     11    dz                         = 25.0, ! Size of single gridbox in z-direction
    1212
    13    dz_stretch_level           = 100.0,
    14    dz_max                     = 50
     13    dz_stretch_level           = 100.0, ! Height (in m) where stretching starts
     14    dz_max                     = 50.0,  ! maximum vertical grid size
    1515!
    1616!-- initialization
    1717!------------------------------------------------------------------------------
    18    initializing_actions       = 'set_constant_profiles',
     18    initializing_actions       = 'set_constant_profiles', ! initial conditions
    1919
    20    ug_surface                 = 2.5,
    21    vg_surface                 = 0.0,
     20    ug_surface                 = 2.5, ! u-comp of geostrophic wind at surface
     21    vg_surface                 = 0.0, ! v-comp of geostrophic wind at surface
    2222
    23    pt_surface                 = 293.0,
    24    pt_vertical_gradient       =    0.0,  3.0,   0,
    25    pt_vertical_gradient_level = 0.0,  800,  1600,
     23    pt_surface                 = 293.0, ! initial surface potential temp
    2624
    27    q_surface                  = 0.013, ! mixing ratio --> rel. hum. = ?
     25    pt_vertical_gradient       = 0.0,
     26                                 3.0,
     27                                 0.0,     ! piecewise temp gradients
     28    pt_vertical_gradient_level = 0.0,
     29                                 800.0,
     30                                 1600.0,  ! height level of temp gradients
     31 
     32    q_surface                  = 0.013,   ! mixing ratio --> rel. hum. = ?
    2833
    29    surface_pressure           = 1000.0,
     34    surface_pressure           = 1000.0,  ! in hPa
    3035
    31    day_of_year_init           = 172, ! June 21
    32    time_utc_init              = 43200.0, ! 12:00 UTC
     36    day_of_year_init           = 172,    ! June 21
     37    time_utc_init              = 43200.0, ! 12:00 UTC
    3338!
    3439!-- physics
    3540!------------------------------------------------------------------------------
    36    latitude                   = 53.37,
    37    longitude                  = 9.73, ! Hannover
     41    latitude                   = 53.37, ! geographical latitude
     42    longitude                  = 9.73,  ! geographical longitude
    3843!
    3944!-- boundary conditions
    4045!------------------------------------------------------------------------------
    41    constant_flux_layer        = .TRUE.,
     46    constant_flux_layer        = .TRUE., ! use a constant flux layer at
     47                                         ! the bottom boundary
    4248
    43 
    44    bc_pt_b                    = 'dirichlet', ! required when using LSM
    45    bc_q_b                     = 'dirichlet', ! required when using LSM
    46 
    47    rayleigh_damping_height    = 1600.0,
    48    rayleigh_damping_factor    = 0.01,
     49    bc_pt_b                    = 'dirichlet', ! required when using LSM
     50    bc_q_b                     = 'dirichlet', ! required when using LSM
    4951!
    5052!-- mode
    5153!------------------------------------------------------------------------------
    52    humidity                   = .TRUE.,
    53    !precipitation              = .FALSE.,
     54    humidity                   = .TRUE.,  ! use prognostic equation for humidity
    5455
    55    reference_state            = 'horizontal_average',
     56    reference_state            = 'horizontal_average', ! set reference state for
     57                                                       ! buoyancy term
    5658!
    5759!-- numerics
    5860!------------------------------------------------------------------------------
    59    fft_method                 = 'temperton-algorithm', !'fftw', !
    60 /
     61    fft_method                 = 'temperton-algorithm', ! build-in fft method
    6162
     63    rayleigh_damping_height    = 1600.0, ! starting height for Rayleigh damping
     64    rayleigh_damping_factor    = 0.01,   ! strength of Rayleigh damping
     65
     66/ ! end of initialization parameter namelist
     67
     68!-------------------------------------------------------------------------------
     69!-- RUNTIME PARAMETER NAMELIST
     70!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par
     71!-------------------------------------------------------------------------------
    6272&runtime_parameters
    6373!
    6474!-- run steering
     75!-------------------------------------------------------------------------------
     76    end_time                   = 216000.0, ! simulation time of the 3D model
     77                                           ! (2.5 days)
     78
     79    create_disturbances        = .T.,    ! randomly perturbate horiz. velocity
     80    dt_disturb                 = 120.0,  ! interval for random perturbations
     81    disturbance_energy_limit   = 0.0001, ! upper limit for perturbation energy
     82!
     83!-- data output
    6584!------------------------------------------------------------------------------
    66    end_time                   = 216000.0, ! 2.5 days
     85    netcdf_data_format         = 2, ! use NetCDF3
    6786
    68    create_disturbances        = .T.,
    69    dt_disturb                 = 120.0,
    70    disturbance_energy_limit   = 0.0001,
    71 !
    72 !-- general output setting
    73 !------------------------------------------------------------------------------
    74    netcdf_data_format         = 2, ! NetCDF3
    75    dt_run_control             = 900.0,
    76    dt_data_output             = 21600.0, ! 6 hours
    77    dt_data_output_av          = 3600.0,
    78    averaging_interval         = 1800.0,
    79 !
    80 !-- profile output setting
    81 !------------------------------------------------------------------------------
    82    data_output_pr             = '#u', '#v', '#theta', '#q', '#thetav',
    83                                 '#km', '#kh', '#l',
    84                                 'w', 'e', 'e*', 'p',
    85                                 'w"u"', 'w*u*', 'w"v"', 'w*v*',
    86                                 'w"theta"', 'w*theta*', 'w"thetav"', 'w*thetav*',
    87                                 'w"q"', 'w*q*', 'w*e*',
    88                                 'theta*2', 'q*2', 'u*2', 'v*2', 'w*2',
    89                                 '#t_soil', '#m_soil',
    90 !
    91 !-- 2D and 3D output setting
    92 !------------------------------------------------------------------------------
    93    section_xy                 = 0, ! surface variables only
     87    dt_run_control             = 900.0,   ! output interval for run control
     88    dt_data_output             = 21600.0, ! output interval for general data
     89    dt_data_output_av          = 3600.0,  ! output interval for averaged data
    9490
    95    data_output                = 'theta', 'thetav', 'q', 'u', 'v',
    96                                 'theta_av', 'thetav_av', 'q_av', 'u_av', 'v_av',
     91    averaging_interval         = 1800.0,  ! averaging interval general data
    9792
    98                                 't_soil',  't_soil_av',
    99                                 'm_soil',  'm_soil_av',
     93    data_output_pr             = '#u', '#v', '#theta', '#q', '#thetav',
     94                                 '#km', '#kh', '#l',
     95                                 'w', 'e', 'e*', 'p',
     96                                 'w"u"', 'w*u*', 'w"v"', 'w*v*',
     97                                 'w"theta"', 'w*theta*', 'w"thetav"',
     98                                 'w*thetav*',
     99                                 'w"q"', 'w*q*', 'w*e*',
     100                                 'theta*2', 'q*2', 'u*2', 'v*2', 'w*2',
     101                                 '#t_soil', '#m_soil',  ! Profile
     102                                                        ! output
    100103
    101                                 'us*_xy',   'us*_xy_av',
    102                                 't*_xy',   't*_xy_av',
    103                                 'r_a*_xy', 'r_a*_xy_av',
    104                                 'r_s*_xy', 'r_s*_xy_av',
    105                                 'tsurf*_xy', 'tsurf*_xy_av'
     104    section_xy                 = 0, ! surface variables only
    106105
    107                                 'ghf*_xy', 'ghf*_xy_av',
    108                                 'shf*_xy', 'shf*_xy_av',
     106    data_output                = 'theta',         'theta_av',
     107                                 'thetav',        'thetav_av',
     108                                 'q',             'q_av',
     109                                 'u',             'u_av',
     110                                 'v',             'v_av',
    109111
    110                                 'qsws_liq*_xy',  'qsws_liq*_xy_av',
    111                                 'qsws_soil*_xy', 'qsws_soil*_xy_av',
    112                                 'qsws_veg*_xy',  'qsws_veg*_xy_av',
     112                                 't_soil',        't_soil_av',
     113                                 'm_soil',        'm_soil_av',
    113114
    114                                 'c_liq*_xy',  'c_liq*_xy_av',
    115                                 'c_soil*_xy', 'c_soil*_xy_av',
    116                                 'c_veg*_xy',  'c_veg*_xy_av',
    117 /
     115                                 'us*_xy',        'us*_xy_av',
     116                                 't*_xy',         't*_xy_av',
     117                                 'r_a*_xy',       'r_a*_xy_av',
     118                                 'r_s*_xy',       'r_s*_xy_av',
     119                                 'tsurf*_xy',     'tsurf*_xy_av'
    118120
     121                                 'ghf*_xy',       'ghf*_xy_av',
     122                                 'shf*_xy',       'shf*_xy_av',
     123
     124                                 'qsws_liq*_xy',  'qsws_liq*_xy_av',
     125                                 'qsws_soil*_xy', 'qsws_soil*_xy_av',
     126                                 'qsws_veg*_xy',  'qsws_veg*_xy_av',
     127
     128                                 'c_liq*_xy',     'c_liq*_xy_av',
     129                                 'c_soil*_xy',    'c_soil*_xy_av',
     130                                 'c_veg*_xy',     'c_veg*_xy_av', ! 2d and/
     131                                                                  ! or 3d
     132                                                                  ! output
     133/ ! end of runtime parameter namelist
     134
     135!-------------------------------------------------------------------------------
     136!-- LAND-SURFACE PARAMETER NAMELIST
     137!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/lsmpar
     138!-------------------------------------------------------------------------------
    119139&land_surface_parameters
    120 !
    121 !-- soil setup
    122 !------------------------------------------------------------------------------
    123    soil_type                  = 3,
    124    soil_temperature           = 293.0, 293.0, 293.0, 293.0,
    125                                 293.0, 293.0, 293.0, 293.0,
    126    deep_soil_temperature      = 293.0,
    127    soil_moisture              = 0.35,  0.35,  0.35,  0.35, 
    128                                 0.35,  0.35,  0.35,  0.35,
    129 !
    130 !-- boundary conditions
    131 !------------------------------------------------------------------------------
    132    surface_type               = 'vegetation',
    133    vegetation_type            = 3,
    134    vegetation_coverage        = 0.95,
    135    c_surface                  = 0.0,
    136    conserve_water_content     = .T.,
    137 /
     140    soil_type                  = 3,  ! type to be used in soil model
     141    soil_moisture              = 0.35,  0.35,  0.35,  0.35, ! moisture at every
     142                                 0.35,  0.35,  0.35,  0.35, ! soil level
     143    soil_temperature           = 293.0, 293.0, 293.0, 293.0, ! temperature at
     144                                 293.0, 293.0, 293.0, 293.0, ! every soil level
     145    deep_soil_temperature      = 293.0,  ! temperature of deep soil
    138146
     147    surface_type               = 'vegetation',  ! surface classification
     148    vegetation_type            = 3, ! vegetation to be used in surface model
     149
     150    vegetation_coverage        = 0.95, ! part of surface covered by vegetation
     151                                       ! (0 - 1)
     152    c_surface                  = 0.0,  ! surface heat capacity
     153    conserve_water_content     = .T.,  ! closed (T) or open (F) lower boundary
     154
     155/ ! end of land-surface parameter namelist
     156
     157!-------------------------------------------------------------------------------
     158!-- RADIATION PARAMETER NAMELIST
     159!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/radpar
     160!-------------------------------------------------------------------------------
    139161&radiation_parameters
    140 !
    141 !-- general setup
    142 !------------------------------------------------------------------------------
    143    radiation_scheme           = 'clear-sky',
    144    dt_radiation               = 60.0,
    145 /
     162    radiation_scheme           = 'clear-sky', ! Radiation scheme to be used
     163    dt_radiation               = 60.0, ! Time step of the radiation model
     164
     165/ ! end of radiation parameter namelist
Note: See TracChangeset for help on using the changeset viewer.