source: palm/trunk/TUTORIALS/cases/chem_smog/INPUT/chem_smog_p3d

Last change on this file was 4765, checked in by Giersch, 4 years ago

Tutorials updated for r4761

  • Property svn:executable set to *
File size: 9.5 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!-- model grid
8!--------------------------------------
9    nx = 35,
10    ny = 35,
11    nz = 80,
12
13    dx = 10.0,
14    dy = 10.0,
15    dz = 10.0,
16
17    dz_stretch_level = 1200.0, dz_stretch_factor = 1.02, dz_max = 60.0,
18
19
20!-- pressure solver
21!--------------------------------------
22!   fft_method =  'fftw',
23
24      psolver = 'multigrid',
25      cycle_mg  = 'w',
26      mg_cycles = 2,
27      mg_switch_to_pe0_level = -1,
28
29!
30!-- initial and boundary conditions
31!--------------------------------------
32    initializing_actions = 'set_constant_profiles',  ! options: 'set_constant_profiles', or 'inifor',
33                                                     ! or both, e.g. 'inifor set_constant_profiles'
34    constant_flux_layer = .T.,
35!
36!-- wind speed
37!--------------------------------------
38     ug_surface = 1.0,
39     vg_surface = 0.0,
40!
41!-- temperature
42!--------------------------------------
43    pt_surface = 290.0,
44    bc_pt_b    = 'dirichlet',
45    pt_vertical_gradient      = 0.1, 1.0,     !pt gradient/100m K
46    pt_vertical_gradient_level= 0.0, 500.0,   !height (m) at which pt gradient is effective
47!
48!-- humidity
49!--------------------------------------
50    humidity   = .T.,
51    q_surface  = 0.001,
52    bc_q_b     = 'dirichlet',
53!
54!-- time
55!--------------------------------------
56    origin_date_time     = '2017-07-31 06:00:00 +00'
57!
58!-- wall/soil spinup
59!--------------------------------------
60     spinup_time         = 86400.0,         ! 24h wall/soil spinup
61     spinup_pt_mean      = 290.0,
62     spinup_pt_amplitude = 20.0,
63     dt_spinup           = 120.0,           ! time step during spinup
64     data_output_during_spinup = .F.,       ! set .T. if you want output
65!
66!-- building topology
67!--------------------------------------
68    topography = 'read_from_file',
69
70/
71!------------------------------------------------------------------------------
72!-- RUNTIME PARAMETER NAMELIST
73!-- Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par
74!------------------------------------------------------------------------------
75 &runtime_parameters
76!
77!-- run time
78!--------------------------------------
79    end_time                =  7200.,
80    termination_time_needed =  600.0,   
81!
82!-- processor topology
83!--------------------------------------
84!    npex = 6,
85!    npey = 6,
86!
87!-- 2d/3d output
88!--------------------------------------
89    do3d_at_begin         = .T.,
90    do2d_at_begin         = .T.,
91    dt_data_output        = 300.0,
92    dt_data_output_av     = 600.0,
93    dt_do2d_xy            = 600.0,
94    section_xy            = 0,
95    averaging_interval    = 600.0,
96
97     
98    data_output = 'theta', 'q', 'u', 'v', 'w','e',
99                   't_soil', 'm_soil',
100                   'tsurf*_xy', 'shf*_xy', 'qsws*_xy', 'ghf*_xy',
101                   'rad_net*_xy', 'us*_xy', 't*_xy', 'r_a*_xy', 'm_liq*_xy',
102                   'kc_NO2','kc_O3','kc_NO', 'kc_RH', 'kc_RCHO','kc_CO',
103                   'kc_OH','kc_HO2','kc_HNO3',
104                   'theta_av', 'q_av', 'u_av', 'v_av', 'w_av','e_av',
105                   'tsurf*_xy_av', 'shf*_xy_av', 'qsws*_xy_av', 'ghf*_xy_av',
106                   'kc_NO2_av','kc_O3_av','kc_NO_av',
107!
108! Masked output - here for first layer above ground
109   mask_k_over_surface(1,:) = 0,          ! 1. layer above surface (terrain following)
110   dt_domask        = 60.,                ! time interval for masked output
111   data_output_masks(1,:) = 'u', 'v', 'kc_NO', 'kc_NO2', 'kc_O3',
112!
113!-- timeseries output
114!--------------------------------------
115    dt_run_control =   0.0,
116    dt_dots        =   60.0,
117!
118!-- profile output
119!--------------------------------------
120    skip_time_dopr        =   0.0,
121    dt_dopr               =  600.0,
122    averaging_interval_pr =  600.0,
123    dt_averaging_input_pr =  600.0,
124
125    data_output_pr = 'e', 'e*',
126                     '#u', 'w"u"', 'w*u*', 'wu', 'u*2',
127                     '#v', 'w"v"', 'w*v*', 'wv', 'v*2',
128                     'w',  'w*2',
129                     '#theta', 'w"theta"', 'w*theta*','wtheta', 'theta*2',
130                     'kc_NO2','kc_O3','kc_NO',
131 /
132
133!-------------------------------------------------------------------------------
134!-- RADIATION MODEL PARAMETER NAMELIST
135!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/radpar
136!-------------------------------------------------------------------------------
137 &radiation_parameters
138    radiation_scheme =  'clear-sky',    !'clear-sky', 'rrtmg',
139    dt_radiation     = 60.0,
140
141    albedo_type      = 5,
142    constant_albedo  = .F.,
143
144    surface_reflections = .T.,
145    nrefsteps = 5,
146
147 / ! end of radiation_parameters namelist
148
149!------------------------------------------------------------------------------
150!-- URBAN SURFACE MODEL PARAMETER NAMELIST
151!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/usmpar
152!------------------------------------------------------------------------------
153 &urban_surface_parameters
154
155 / ! end of urban_surface_parameters namelist
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    vegetation_type = 2,
164    soil_type       = 3,
165
166    root_fraction = 0.49, 0.26, 0.25, 0.0, 0.0, 0.0, 0.0, 0.0,
167    soil_moisture = 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3,
168    soil_temperature = 289., 289., 291., 292., 292., 290., 288., 287.,
169    deep_soil_temperature = 285.,
170
171    aero_resist_kray       = .T.,
172    constant_roughness     = .T.,
173    conserve_water_content = .T.,
174 / ! end of land_surface_parameters namelist
175
176
177!------------------------------------------------------------------------------
178!-- PLANT CANOPY MODEL PARAMETER NAMELIST
179!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/canpar
180!------------------------------------------------------------------------------
181 &plant_canopy_parameters
182    canopy_mode = 'read_from_file',
183    canopy_drag_coeff = 0.3,
184 / ! end of plant_canopy_parameters namelist
185
186!------------------------------------------------------------------------------
187!-- CHEMISTRY MODEL PARAMETER NAMELIST
188!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/chempar
189!------------------------------------------------------------------------------
190 &chemistry_parameters
191
192    chem_gasphase_on           = .TRUE.,
193    chem_mechanism             = "smog",
194    deposition_dry             = .TRUE.,
195    photolysis_scheme          = 'simple',
196
197! initial profiles
198    cs_name                    = 'O3',  'NO',  'NO2',  'RH','RCHO',  'CO',
199    cs_surface                 = 0.025, 0.000, 0.010, 0.010, 0.005, 0.100,  !units- chem spcs in ppm, and PM10 in ug/m**2
200
201    cs_profile(1,:)            = 0.025, 0.040, 0.050,
202    cs_profile(2,:)            = 0.000, 0.000,
203    cs_profile(3,:)            = 0.010, 0.002,
204    cs_profile(4,:)            = 0.010, 0.005,
205    cs_profile(5,:)            = 0.005, 0.0025,
206    cs_profile(6,:)            = 0.100,
207
208    cs_heights(1,:)            = 0.0,  305.0,  495.0,
209    cs_heights(2,:)            = 0.0,   295.0,
210    cs_heights(3,:)            = 0.0,   295.0,
211    cs_heights(4,:)            = 0.0,   295.0,
212    cs_heights(5,:)            = 0.0,   295.0,
213    cs_heights(6,:)            = 0.0, 
214
215! Emissions
216    emissions_anthropogenic    = .TRUE.,
217    emiss_lod                  = 0,                 ! 0 "PARAMETERIZED", 1 "DEFAULT", 2 "PRE-PROCESSED",
218! Parameters for parameterized traffic emissions
219    surface_csflux_name        = 'NO',   'NO2',   'RH', 'RCHO', 'CO',   !"cs" stands for chemical species,
220                                                            ! only for parameterized traffic emissions
221                                                            ! NOT to be used when running with emiss_lod>0
222    surface_csflux             = 4745.0, 1326.0, 649., 30., 30000., ! surface emission of NO, NO2 in umol m-2 d-1 (PMs in kg m-2 d-1),only for parameterized traffic emissions
223    emiss_factor_main          = 1.667, 1.667, 1.667, 1.667, 1.667,      ! emiss_factor_main*surface_csflux for emission on main streets, one factor per surface_csflux component needed!
224    emiss_factor_side          = 0.334, 0.334, 0.334, 0.334, 0.334,      ! emiss_factor_side*surface_csflux for emission on side streets, one factor per surface_csflux component needed!
225    main_street_id             = 11,                        ! id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90)
226    side_street_id             =  8,
227    max_street_id              = 19,
228
229! Bottom and top boundary conditions (lateral BCs are cyclic, which is the default)
230    bc_cs_b                    = 'neumann',
231    bc_cs_t                    = 'neumann',
232
233! Parameters chemistry solver (see KPP description)
234    icntrl(3)                  = 1,    ! solver ros2
235    icntrl(4)                  = 500,  ! max number of chem-substeps
236    rcntrl(3)                  = 0.1,  ! Hstart, starting value for the integration step size
237
238 / ! end of chemistry_parameters namelist
Note: See TracBrowser for help on using the repository browser.