source: palm/trunk/TESTS/cases/urban_environment/INPUT/urban_environment_p3d @ 3653

Last change on this file since 3653 was 3653, checked in by forkel, 6 years ago

parameter chem_mechanism added to &chemistry_parameters

  • Property svn:executable set to *
File size: 9.9 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!
8!-- Grid
9!------------------------------------------------------------------------------
10    nx                   = 19,
11    ny                   = 19,
12    nz                   = 60,
13
14    dx                   = 2.0,
15    dy                   = 2.0,
16    dz                   = 2.0,
17
18!
19!-- Numerics
20!------------------------------------------------------------------------------
21    fft_method           = 'fftw',
22
23!
24!-- Mode
25!------------------------------------------------------------------------------     
26    humidity             = .T.,
27    q_surface            = 0.001,
28
29!
30!-- Initialization
31!------------------------------------------------------------------------------
32    initializing_actions = 'inifor set_constant_profiles',   !'set_constant_profiles', or 'inifor',
33                                                             ! or both (met. from inifor, chem constant profiles)
34
35    pt_surface           = 277.15,
36    ug_surface           = 1.0,
37    vg_surface           = 0.0,
38
39    day_of_year_init     = 65,             ! march 06, to avoid extreme temperatures
40    time_utc_init        = 36000.0,        ! (0.0 s = 00 UTC midnight),
41
42!
43!-- Wall/soil spinup
44!------------------------------------------------------------------------------   
45    spinup_time         = 1200.0,
46    spinup_pt_mean      = 284.15,
47    spinup_pt_amplitude = 10.0,             ! give a dirunal cycle from 6 - 16 °C
48    dt_spinup           = 120.0,
49    data_output_during_spinup = .F.,
50
51!
52!-- Topography
53!------------------------------------------------------------------------------   
54    topography           = 'read_from_file',
55
56!
57!-- Physics
58!------------------------------------------------------------------------------
59    longitude            = 13.4,
60    latitude             = 52.5,
61
62 / ! end of inipar namelist
63
64
65!-------------------------------------------------------------------------------
66!-- RUNTIME PARAMETER NAMELIST
67!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par
68!-------------------------------------------------------------------------------
69 &runtime_parameters
70
71!
72!-- Run steering
73!------------------------------------------------------------------------------
74    end_time              = 70.0,
75    create_disturbances   = .T.,
76
77!
78!-- Run-control/timeseries output settings
79!------------------------------------------------------------------------------
80    dt_run_control        =  0.0,
81    dt_dots               =  10.0,
82
83!
84!-- Profile output settings
85!------------------------------------------------------------------------------
86    skip_time_dopr        =  0.0,
87    dt_dopr               = 60.0,
88    averaging_interval_pr = 60.0,
89    dt_averaging_input_pr =  0.0,
90
91    data_output_pr        = '#u',  'u*2',  'wu',  'w*u*',  'w"u"',
92                            '#v',  'v*2',  'wv',  'w*v*',  'w"v"',
93!                            'w',   'w*2',
94                            '#theta', 'theta*2', 'wtheta', 'w*theta*', 'w"theta"',
95                            '#q',  'q*2',
96                            'e', 'e*', '#km', '#l',
97                            '#t_soil', '#m_soil', 'rad_lw_in', 'rad_lw_out',
98                            'kc_NO', 'kc_NO2', 'kc_O3', 'kc_PM10',
99!
100!-- 2D/3D output settings
101!------------------------------------------------------------------------------
102    do3d_at_begin      = .T.,
103    do2d_at_begin      = .T.,
104
105    dt_data_output     = 60.0,
106    dt_data_output_av  = 60.0,
107    averaging_interval = 60.0,
108    dt_averaging_input =  0.0,
109
110    section_xy         = 0,
111
112    data_output = 'u',  'u_av',
113                  'v',  'v_av',
114                  'w',  'w_av',
115                  'theta', 'theta_av',
116                  'q',  'q_av',
117                  'e',  'e_av',
118                  'p',  'p_av',
119                  'm_soil', 'm_soil_av',
120                  't_soil', 't_soil_av',
121                  'shf*_xy', 'shf*_xy_av',
122                  'ghf*_xy', 'ghf*_xy_av',
123                  'us*_xy', 'z0*_xy',
124                  't*_xy', 'ol*_xy',
125                  'm_liq*_xy',
126                  'r_s*_xy', 'r_a*_xy',
127                  'rad_net*_xy', 'rad_lw_in*_xy', 'rad_lw_out*_xy', 'rad_sw_in*_xy', 'rad_sw_out*_xy',
128                  'qsws*_xy', 'qsws*_xy_av',
129                  'tsurf*_xy', 'tsurf*_xy_av',
130                  'pcm_lad', 'pcm_heatrate',
131                  'kc_NO', 'kc_NO2', 'kc_O3', 'kc_PM10',
132                  'kc_NO_av', 'kc_NO2_av', 'kc_O3_av', 'kc_PM10_av',
133
134 / ! end of d3par namelist
135
136
137!-------------------------------------------------------------------------------
138!-- RADIATION MODEL PARAMETER NAMELIST
139!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/radpar
140!-------------------------------------------------------------------------------
141  &radiation_parameters
142
143     radiation_scheme =  'rrtmg',    !'clear-sky' or 'rrtmg',
144     albedo_type = 5,
145     constant_albedo = .F.,
146
147     dt_radiation = 60.0,
148
149     surface_reflections = .T.,
150     nrefsteps = 3,
151
152     rad_angular_discretization = .F.,
153     raytrace_mpi_rma = .F.,
154
155  /
156
157
158!-------------------------------------------------------------------------------
159!-- LAND SURFACE MODEL PARAMETER NAMELIST
160!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/lsmpar
161!-------------------------------------------------------------------------------
162 &land_surface_parameters
163
164    constant_roughness     = .T.,
165    aero_resist_kray       = .T.,
166
167    vegetation_type        = 2,
168    soil_type              = 3,
169    conserve_water_content = .T.,
170
171    root_fraction          = 0.49, 0.26, 0.25, 0.0, 0.0, 0.0, 0.0, 0.0,
172    soil_temperature       = 277.15, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0,
173    soil_moisture          = 0.30,   0.30,  0.30,  0.30,  0.30, 0.30,  0.30,  0.30,
174    deep_soil_temperature  = 277.0,
175   
176 / ! end of lsm_par namelist
177
178
179!-------------------------------------------------------------------------------
180!-- URBAN SURFACE MODEL PARAMETER NAMELIST
181!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/usmpar
182!-------------------------------------------------------------------------------
183 &urban_surface_parameters
184
185    usm_material_model     = .T.,
186    usm_wall_mod = .T.,
187
188 / ! end of urban_surface_par namelist
189 
190
191!-------------------------------------------------------------------------------
192!-- PLANT CANOPY MODEL PARAMETER NAMELIST
193!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/canpar
194!-------------------------------------------------------------------------------
195 &plant_canopy_parameters
196 
197    canopy_mode       = 'read_from_file_3d',
198    canopy_drag_coeff = 0.3,
199     
200 / ! end of canopy_par namelist
201 
202
203!-------------------------------------------------------------------------------
204!-- CHEMISTRY MODEL PARAMETER NAMELIST
205!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/chempar
206!-------------------------------------------------------------------------------
207 &chemistry_parameters
208
209    chem_mechanism             = 'phstatp',
210    chem_gasphase_on           = .T.,
211    call_chem_at_all_substeps  = .FALSE.,
212   
213    photolysis_scheme          = 'simple',
214
215    cs_name                    = 'NO', 'NO2', 'O3', 'PM10',
216    cs_surface                 = 0.02,  0.03, 0.02, 25.0e-9,   ! surface concentration of PM10, PM25 in kg m-3,
217                                                               ! units for gases in ppm   
218!-- initial profiles (also for lateral boundaries if 'set_constant_profiles' is chosen)
219!-- units for PM: kg m-3, units for gases: ppm
220    cs_profile(1,:)            =  0.020, 0.020, 0.015, 0.015, 0.015, 0.010, 0.010, 0.010, 0.010, 0.010, 0.007,
221    cs_profile(2,:)            =  0.030, 0.030, 0.020, 0.020, 0.020, 0.015, 0.015, 0.015, 0.015, 0.015, 0.010,
222    cs_profile(3,:)            =  0.020, 0.020, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050,
223    cs_profile(1,:)            =  25.e-09, 25.e-09, 23.0e-09, 23.e-09, 23.e-09, 23.e-09, 23.e-09, 23.e-09, 22.e-09, 22.e-09, 22.0e-09,
224
225    cs_heights(1,:)            =  0.0,   5.0,   15.0,   25.0,  35.0,  45.0,  55.0,  65.0,  75.0,  85.0, 95.0,
226    cs_heights(2,:)            =  0.0,   5.0,   15.0,   25.0,  35.0,  45.0,  55.0,  65.0,  75.0,  85.0, 95.0,
227    cs_heights(3,:)            =  0.0,   5.0,   15.0,   25.0,  35.0,  45.0,  55.0,  65.0,  75.0,  85.0, 95.0,
228    cs_heights(4,:)            =  0.0,   5.0,   15.0,   25.0,  35.0,  45.0,  55.0,  65.0,  75.0,  85.0, 95.0,
229
230    bc_cs_b                    = 'neumann',
231    bc_cs_t                    = 'neumann',
232   
233    do_emis                    = .T.,
234    mode_emis                  = 'PARAMETERIZED',  ! emission depending on street type from static driver
235    surface_csflux_name        = 'NO',   'NO2',  'PM10',   !"cs" stands for chemical species
236    surface_csflux             = 4745.0, 1326.0, 2700.e-09,! surface emission of PM10, PM25 in kg m-2 d-1,
237                                                           ! for gases in umol m-2 d-1 for PARAMETERIZED
238    emiss_factor_main          = 1.667,  1.667,  1.667,    ! emiss_factor_main*surface_csflux for emission on main streets
239    emiss_factor_side          = 0.334,  0.334,  0.334,    !emiss_factor_side*surface_csflux for emission on side streets
240    main_street_id = 12,       ! id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90)
241    side_street_id =  8,
242    max_street_id  = 19,
243
244! -- KPP solver steering. For details  see comments in chem_gasphase_mod.f90, SUBROUTINE rosenbrock
245!                         (or other solver subroutine, if another solver is chosen)
246    icntrl(3)                  = 1,   ! solver ros2
247    icntrl(4)                  = 500, ! max. number of chem-substeps
248    rcntrl(3)                  = 0.1, ! Hstart,starting value for the integration step size
249
250 / ! end of chemistry_par namelist
Note: See TracBrowser for help on using the repository browser.