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/dispersion_eulerian_and_lpm_extended/INPUT/dispersion_eulerian_and_lpm_extended_p3d

    r4098 r4224  
     1!-------------------------------------------------------------------------------
     2!-- INITIALIZATION PARAMETER NAMELIST
     3!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar
     4!-------------------------------------------------------------------------------
    15&initialization_parameters
     6!
     7!-- grid parameters
     8!-------------------------------------------------------------------------------
     9    nx                      = 39, ! 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)
    212
    3 !-- grid setup
    4     nx = 39,
    5     ny = 39,
    6     nz = 40,
     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    = 'set_constant_profiles', ! initial conditions
    720
    8     dx = 2.0,
    9     dy = 2.0,
    10     dz = 2.0,
     21    ug_surface              = 1.0, ! u-comp of geostrophic wind at surface
     22    vg_surface              = 0.0, ! v-comp of geostrophic wind at surface
     23   
     24    s_surface               = 0.0008, ! equivalent to 420ppm CO2 as background
     25                                   ! concentration
     26!
     27!-- boundary conditions
     28!-------------------------------------------------------------------------------
     29    bc_uv_t                 = 'neumann', ! free-slip boundary condition
     30   
     31    surface_scalarflux      = 0.02, ! kg emission per square meter and second
     32    bc_s_b                  = 'neumann', ! bottom boundary condition for scalar
     33                                         ! concentration. Required with scalar
     34                                         !flux
     35                                   
     36    roughness_length        = 0.1, ! roughness length (in m) for momentum
     37!
     38!-- numerics
     39!-------------------------------------------------------------------------------
     40    fft_method              = 'temperton-algorithm', ! build-in fft method
     41     
     42    momentum_advec          = 'ws-scheme', ! default advection scheme
     43    scalar_advec            = 'ws-scheme', ! default advection scheme
     44!
     45!-- mode
     46!-------------------------------------------------------------------------------
     47    neutral                 = .TRUE., ! strictly neutral flow
    1148
    12 !-- init-action
    13     initializing_actions = 'set_constant_profiles',
     49    passive_scalar          = .TRUE., ! switch on prognostic equation for pas.
     50                                      ! scalar
    1451
    15 !-- advection schemes
    16     momentum_advec = 'ws-scheme',           ! default advection scheme
    17     scalar_advec   = 'ws-scheme',
     52    dp_external             = .TRUE., ! use ext. pressure gradient as driving force
     53    dpdxy                   = -0.0002, 0.0, ! external pressure gradient applied
     54                                            ! in x- and y-direction (Pa/m)
     55!
     56!-- physics
     57!-------------------------------------------------------------------------------
     58    omega                   = 0.0, ! angular velocity of the rotating system,
     59                                   ! no Coriolis force
     60!
     61!-- topography
     62!-------------------------------------------------------------------------------
     63    topography              = 'read_from_file', ! ASCII Input file contains the 
     64                                                ! arbitrary topography height
     65                                                ! information in m.
     66                                   
     67/ ! end of initialization parameter namelist
    1868
    19 !-- wind speed setup
    20     ug_surface = 1.0,
    21     vg_surface = 0.0,
     69!-------------------------------------------------------------------------------
     70!-- RUNTIME PARAMETER NAMELIST
     71!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par
     72!-------------------------------------------------------------------------------
     73&runtime_parameters
     74!
     75!-- run steering
     76!-------------------------------------------------------------------------------
     77    end_time                = 10800.0, ! simulation time of the 3D model
    2278
    23 !-- temperature setup
    24     neutral = .T.,                          ! strictly neutral flow
     79    create_disturbances     = .TRUE., ! randomly perturbate horiz. velocity
     80!
     81!-- data output
     82!-------------------------------------------------------------------------------
     83    skip_time_data_output   =    0.0, ! no data output before this time interval
     84    dt_data_output          = 3600.0, ! output interval for general data
    2585
    26 !-- tracer setup
    27     passive_scalar     = .T.,
    28     s_surface          = 0.0008,            !equivalent to 420ppm CO2 as background concentration
    29     surface_scalarflux = 0.02,              !kg emission per square meter and second
    30     bc_s_b             = 'neumann',
     86    averaging_interval      = 3600.0, ! averaging interval general data
     87    dt_averaging_input      =    6.0, ! general data sampling rate for averaging
    3188
    32 !-- boundary conditions
    33     bc_uv_t = 'neumann',                    ! free-slip boundary condition
    34     roughness_length = 0.1,                 ! applies to all surfaces locally
     89    dt_run_control          =  6.0, ! output interval for run control
     90    dt_dots                 =  6.0, ! output interval for time series
    3591
    36 !-- special methods
    37     dp_external = .T.,                      ! use horizontal pressure gradient
    38     dpdxy = -0.0002, 0.0,                   ! for forcing
     92    skip_time_dopr          =    0.0, ! no profile data output before this time
     93                                      ! interval
     94    dt_dopr                 = 3600.0, ! output interval for profile data
     95    averaging_interval_pr   = 3600.0, ! averaging interval profile data
     96    dt_averaging_input_pr   =    6.0, ! profile data sampling rate for averaging
    3997
    40 !-- special setups
    41     omega = 0.0,                            ! no Coriolis force
     98    do3d_at_begin           = .TRUE., ! output of 3d volume data at the beginning
     99    do2d_at_begin           = .TRUE., ! output of 2d data at the beginning
    42100
    43     topography = 'read_from_file',
    44     topography_grid_convention = 'cell_center',   ! default
    45 /
     101    dt_do2d_xy              = 3600.0, ! output interval for xy cross sections
     102    dt_do2d_xz              = 3600.0, ! output interval for xz cross sections
     103
     104    data_output_pr          = '#u', 'u*2', 'wu', 'w*u*', 'w"u"',
     105                              '#v', 'v*2', 'wv', 'w*v*', 'w"v"',
     106                              'w', 'w*2',
     107                              'e', 'e*',
     108                              '#km',
     109                              '#l',
     110                              '#s', 'ws', 'w*s*', 'w"s"', ! profile output
     111
     112    data_output             = 'u', 'u_av',
     113                              'u_xy', 'u_xz', 'u_xz_av',
     114                              'v', 'v_av',
     115                              'w', 'w_av',
     116                              's', 's_av',
     117                              'ssws*_xy',
     118                              'pc', 'pc_av',
     119                              'pc_xy', 'pc_xz', 'pc_xz_av', ! 2d and/or 3d
     120                                                            ! output
     121                 
     122    section_xy              = -1, 21, 22, 23, ! grid index for 2D XY cross
     123                                              ! sections
     124    section_xz              = -1, 17, 18, 19, ! grid index for 2D XZ cross
     125                                              ! sections
     126
     127/ ! end of runtime parameter namelist
     128
     129!-------------------------------------------------------------------------------
     130!-- PARTICLE PARAMETER NAMELIST
     131!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/parpar
     132!-------------------------------------------------------------------------------
     133&particle_parameters
     134!
     135!-- boundary conditions for particle transport
     136!-------------------------------------------------------------------------------
     137    bc_par_lr               = 'absorb', ! lateral boundary condition
     138    bc_par_b                = 'reflect', ! bottom boundary condition
     139    bc_par_t                = 'reflect', ! top boundary condition
     140!
     141!-- define point source
     142!-------------------------------------------------------------------------------
     143    psb                     = 45.0, ! bottom edge of a particle source 
     144    pst                     = 45.0, ! top edge of a particle source
     145    pss                     = 40.0, ! front ("south") edge of a particle source
     146    psn                     = 40.0, ! rear ("north") edge of a particle source
     147    psl                     = 25.0, ! left edge of a particle source
     148    psr                     = 25.0, ! right edge of a particle source
     149   
     150    particles_per_point     = 100, ! particles per point (default = 1)
     151   
     152    random_start_position   = .FALSE., ! add random fluctuations for initial
     153                                       ! particle position
     154!
     155!-- define time variables of particle relaese
     156!-------------------------------------------------------------------------------
     157    dt_prel                 = 1.0, ! time after new particle set is released
     158   
     159    step_dealloc            = 50, ! time steps after arrays are deallocated
     160   
     161    particle_maximum_age    = 200.0, ! maximum allowed age of particles in s
     162!
     163!-- physics of particles
     164!-------------------------------------------------------------------------------
     165    use_sgs_for_particles   = .TRUE., ! SGS-velocities for particles
     166!
     167!-- particle output
     168!-------------------------------------------------------------------------------
     169    dt_dopts                = 10.0, ! output interval for particle quantities
     170 
     171/ ! end of particle parameter namelist
     172
     173!-------------------------------------------------------------------------------
     174!-- USER PARAMETER NAMELIST
     175!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/userpar
     176!-------------------------------------------------------------------------------
     177&user_parameters
     178!
     179!-- definition of emission stripe
     180!-------------------------------------------------------------------------------
     181    emission_stripe_width_y = 6.0, ! width of emission stripe in meters
     182    emission_stripe_orig_y  = 12.0, 62.0, ! origin of emission stripe in meters
     183!
     184!-- defition of time series output
     185!-------------------------------------------------------------------------------
     186    s_ts_pos_x              = 17.0, 61.0, 75.0, ! x position of scalar time
     187                                                ! series output
     188    s_ts_pos_y              = 25.0, 15.0, 41.0, ! y-position of scalar time
     189                                                ! series output
     190           
     191/ ! end of user parameter namelist
    46192
    47193
    48 &runtime_parameters
    49 
    50 !-- simulation time
    51     end_time                = 10800.0,
    52     termination_time_needed =   900.0,
    53 
    54 !-- others
    55     create_disturbances = .T.,
    56 
    57 !-- global output settings
    58     skip_time_data_output =    0.0,
    59     dt_data_output        = 3600.0,
    60 
    61     averaging_interval    = 3600.0,
    62     dt_averaging_input    =    6.0,
    63 
    64 !-- time series / run control
    65     dt_run_control        =  6.0,
    66     dt_dots               =  6.0,
    67 
    68 !-- profiles
    69     skip_time_dopr        =    0.0,
    70     dt_dopr               = 3600.0,
    71     averaging_interval_pr = 3600.0,
    72     dt_averaging_input_pr =    6.0,
    73 
    74     do3d_at_begin         = .TRUE.,
    75     do2d_at_begin         = .TRUE.,
    76 
    77     dt_do2d_xy            = 3600.0,
    78     dt_do2d_xz            = 3600.0,
    79 
    80     data_output_pr = '#u', 'u*2', 'wu', 'w*u*', 'w"u"',
    81                      '#v', 'v*2', 'wv', 'w*v*', 'w"v"',
    82                      'w', 'w*2',
    83                      'e', 'e*',
    84                      '#km',
    85                      '#l',
    86                      '#s', 'ws', 'w*s*', 'w"s"',
    87 
    88 !-- output variables
    89     data_output = 'u', 'u_av',
    90                   'u_xy', 'u_xz', 'u_xz_av',
    91                   'v', 'v_av',
    92                   'w', 'w_av',
    93                   's', 's_av',
    94                   'ssws*_xy',
    95                   'pc', 'pc_av', 'pc_xy', 'pc_xz', 'pc_xz_av',
    96                  
    97 
    98     section_xy = -1,21,22,23,
    99     section_xz = -1,17,18,19,
    100 
    101 /
    102 
    103 &particle_parameters
    104 
    105 !-- Boundary conditions for particles
    106     bc_par_lr = 'absorb',
    107     bc_par_b  = 'reflect',
    108     bc_par_t  = 'reflect',
    109 
    110 !-- Define point source
    111     psb = 45.0, 
    112     pst = 45.0,
    113     pss = 40.0,
    114     psn = 40.0,
    115     psl = 25.0,
    116     psr = 25.0,
    117 
    118 !-- Time after new particle set is released
    119     dt_prel = 1.0,
    120 
    121 !-- Time steps after array are deallocated
    122     step_dealloc = 50,
    123 
    124 !-- Particles per point (default = 1)
    125     particles_per_point = 100,
    126 
    127 !-- Maximum time for particles
    128     particle_maximum_age = 200.0,
    129 
    130 !-- Add random fluctuations for initial particle position
    131     random_start_position = .FALSE.,
    132 
    133 !-- SGS-velocities for particles
    134     use_sgs_for_particles = .TRUE.,
    135 
    136 !-- Particle time series interval
    137     dt_dopts = 10.0,
    138 
    139 /
    140 
    141 &user_parameters
    142     emission_stripe_width_y = 6.0,  !in meters
    143     emission_stripe_orig_y  = 12.0, 62.0,  !in meters
    144     s_ts_pos_x = 17.0, 61.0, 75.0,
    145     s_ts_pos_y = 25.0, 15.0, 41.0,
    146            
    147 /
    148 
    149 
Note: See TracChangeset for help on using the changeset viewer.