source: palm/trunk/TUTORIALS/cases/dispersion_eulerian_and_lpm_extended/INPUT/dispersion_eulerian_and_lpm_extended_p3d

Last change on this file was 4409, checked in by Giersch, 3 years ago

Tutorials updated, all tested with revision 4402

File size: 8.8 KB
Line 
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                      = 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)
12
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
20
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
48
49    passive_scalar          = .TRUE., ! switch on prognostic equation for pas.
50                                      ! scalar
51
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
68
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
78
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
85
86    averaging_interval      = 3600.0, ! averaging interval general data
87    dt_averaging_input      =    6.0, ! general data sampling rate for averaging
88
89    dt_run_control          =  6.0, ! output interval for run control
90    dt_dots                 =  6.0, ! output interval for time series
91
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
97
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
100
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"', '#v', 'v*2', 'wv', 'w*v*', 'w"v"', 'w', 'w*2', 'e', 'e*', '#km', '#l', '#s', 'ws', 'w*s*', 'w"s"',
105
106    data_output             = 'u', 'u_av', 'u_xy', 'u_xz', 'u_xz_av', 'v', 'v_av', 'w', 'w_av', 's', 's_av', 'ssws*_xy', 'pc', 'pc_av', 'pc_xy', 'pc_xz', 'pc_xz_av',
107                 
108    section_xy              = -1, 21, 22, 23, ! grid index for 2D XY cross
109                                              ! sections
110    section_xz              = -1, 17, 18, 19, ! grid index for 2D XZ cross
111                                              ! sections
112
113/ ! end of runtime parameter namelist
114
115!-------------------------------------------------------------------------------
116!-- PARTICLE PARAMETER NAMELIST
117!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/parpar
118!-------------------------------------------------------------------------------
119&particle_parameters
120!
121!-- boundary conditions for particle transport
122!-------------------------------------------------------------------------------
123    bc_par_lr               = 'absorb', ! lateral boundary condition
124    bc_par_b                = 'reflect', ! bottom boundary condition
125    bc_par_t                = 'reflect', ! top boundary condition
126!
127!-- define point source
128!-------------------------------------------------------------------------------
129    psb                     = 45.0, ! bottom edge of a particle source 
130    pst                     = 45.0, ! top edge of a particle source
131    pss                     = 40.0, ! front ("south") edge of a particle source
132    psn                     = 40.0, ! rear ("north") edge of a particle source
133    psl                     = 25.0, ! left edge of a particle source
134    psr                     = 25.0, ! right edge of a particle source
135   
136    particles_per_point     = 100, ! particles per point (default = 1)
137   
138    random_start_position   = .FALSE., ! add random fluctuations for initial
139                                       ! particle position
140!
141!-- define time variables of particle relaese
142!-------------------------------------------------------------------------------
143    dt_prel                 = 1.0, ! time after new particle set is released
144   
145    step_dealloc            = 50, ! time steps after arrays are deallocated
146   
147    particle_maximum_age    = 200.0, ! maximum allowed age of particles in s
148!
149!-- physics of particles
150!-------------------------------------------------------------------------------
151    use_sgs_for_particles   = .TRUE., ! SGS-velocities for particles
152!
153!-- particle output
154!-------------------------------------------------------------------------------
155    dt_dopts                = 10.0, ! output interval for particle quantities
156 
157/ ! end of particle parameter namelist
158
159!-------------------------------------------------------------------------------
160!-- USER PARAMETER NAMELIST
161!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/userpar
162!-------------------------------------------------------------------------------
163&user_parameters
164!
165!-- definition of emission stripe
166!-------------------------------------------------------------------------------
167    emission_stripe_width_y = 6.0, ! width of emission stripe in meters
168    emission_stripe_orig_y  = 12.0, 62.0, ! origin of emission stripe in meters
169!
170!-- defition of time series output
171!-------------------------------------------------------------------------------
172    s_ts_pos_x              = 17.0, 61.0, 75.0, ! x position of scalar time
173                                                ! series output
174    s_ts_pos_y              = 25.0, 15.0, 41.0, ! y-position of scalar time
175                                                ! series output
176           
177/ ! end of user parameter namelist
178
179
Note: See TracBrowser for help on using the repository browser.