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