source: palm/trunk/TESTS/cases/urban_environment_salsa/INPUT/urban_environment_salsa_p3d @ 4117

Last change on this file since 4117 was 4117, checked in by monakurppa, 23 months ago

Apply decycling flags for salsa and set decycling boundary conditions only at the ghost points not at the prognostic grid points. See changeset 4109 for more details.

File size: 10.3 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                   = 22,
13
14    dx                   = 2.0,
15    dy                   = 2.0,
16    dz                   = 2.0,
17
18!
19!-- Numerics
20!------------------------------------------------------------------------------
21    fft_method           = 'fftw',
22    monotonic_limiter_z  = .T.,
23!
24!-- Mode
25!------------------------------------------------------------------------------     
26    humidity             = .T.,
27    q_surface            = 0.001,
28
29!
30!-- Initialization
31!------------------------------------------------------------------------------
32    initializing_actions = '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!
97!-- 2D/3D output settings
98!------------------------------------------------------------------------------
99    do3d_at_begin      = .T.,
100    do2d_at_begin      = .T.,
101
102    dt_data_output     = 10.0,
103    dt_data_output_av  = 10.0,
104    averaging_interval = 10.0,
105    dt_averaging_input =  0.0,
106
107    section_xy         = 0,
108
109    data_output = 'u', 'v', 'w',
110                  'Ntot', 'PM2.5', 'LDSA',
111                  's_OC','g_OCNV','g_OCSV',
112                  'N_bin3', 'm_bin4',
113
114
115 / ! end of d3par namelist
116
117
118!-------------------------------------------------------------------------------
119!-- RADIATION MODEL PARAMETER NAMELIST
120!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/radpar
121!-------------------------------------------------------------------------------
122  &radiation_parameters
123
124     radiation_scheme =  'clear-sky',    !'clear-sky' or 'rrtmg',
125     albedo_type = 5,
126     constant_albedo = .F.,
127
128     dt_radiation = 60.0,
129
130     surface_reflections = .T.,
131     nrefsteps = 3,
132
133     rad_angular_discretization = .F.,
134     raytrace_mpi_rma = .F.,
135
136  /
137
138
139!-------------------------------------------------------------------------------
140!-- LAND SURFACE MODEL PARAMETER NAMELIST
141!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/lsmpar
142!-------------------------------------------------------------------------------
143 &land_surface_parameters
144
145    constant_roughness     = .T.,
146    aero_resist_kray       = .T.,
147
148    vegetation_type        = 2,
149    soil_type              = 3,
150    conserve_water_content = .T.,
151
152    root_fraction          = 0.49, 0.26, 0.25, 0.0, 0.0, 0.0, 0.0, 0.0,
153    soil_temperature       = 277.15, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0,
154    soil_moisture          = 0.30,   0.30,  0.30,  0.30,  0.30, 0.30,  0.30,  0.30,
155    deep_soil_temperature  = 277.0,
156   
157 / ! end of lsm_par namelist
158
159
160!-------------------------------------------------------------------------------
161!-- URBAN SURFACE MODEL PARAMETER NAMELIST
162!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/usmpar
163!-------------------------------------------------------------------------------
164 &urban_surface_parameters
165
166    usm_material_model     = .T.,
167    usm_wall_mod = .T.,
168
169 / ! end of urban_surface_par namelist
170 
171
172!-------------------------------------------------------------------------------
173!-- PLANT CANOPY MODEL PARAMETER NAMELIST
174!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/canpar
175!-------------------------------------------------------------------------------
176 &plant_canopy_parameters
177 
178    canopy_mode       = 'read_from_file_3d',
179    canopy_drag_coeff = 0.3,
180     
181 / ! end of canopy_par namelist
182
183
184!-------------------------------------------------------------------------------
185!-- SALSA PARAMETER NAMELIST
186!-------------------------------------------------------------------------------
187 &salsa_parameters 
188
189!
190!-- Time stepping
191    dt_salsa           =   2.0, ! time step for calculating aerosol processes (s)
192    skip_time_do_salsa =  10.0, ! starting time of SALSA (s)
193
194!
195!-- Concentration initial types: 0 = based on the parameter file, 1 = read from a NetCDF file
196    init_aerosol_type  = 0, ! size distribution
197    init_gases_type    = 0, ! gases
198
199!
200!-- If isdtyp = 0, define the initial aerosol size distribution by dpg, sigmag and n_lognorm
201    dpg       = 13.5E-9, 54.0E-9, 864.1E-9, ! mean diameter per mode (in metres)
202    sigmag    =     1.8,    2.16,     2.21, ! standard deviation per mode
203    n_lognorm =  1.43E9,  4.45E8,    7.0E4, ! number concentration per mode (#/m3)
204
205!
206!-- If igctyp = 0, apply these initial gas concentrations                       
207    H2SO4_init = 5.0E12, ! sulphuric acid (#/m3)
208    HNO3_init  = 3.0E12, ! nitric acid (#/m3)
209    NH3_init   = 6.0E12, ! ammonia (#/m3)       
210    OCNV_init  = 1.0E12, ! non-volatile organic gases (#/m3)
211    OCSV_init  = 1.0E12, ! non-volatile organic gases (#/m3)
212
213!
214!-- List of activated chemical components:
215!-- NOTE! Chemical species have to be activated here even if they are not initially present!
216    listspec = 'OC','NO','NH','','','','',  ! List of actived aerosols       
217    ! listspec = 'SO4','OC','BC','DU','SS','NO','NH',
218
219!
220!-- Sectional presentation of the particle size distribution
221    reglim = 3.0E-9, 1.0E-8, 2.5E-6, ! limits of the subranges (m)
222    nbin   = 1, 7,                   ! number of bins per subrange
223!-- NOTE! Subrange 1 consists only of H2SO4 and/or OC
224    nf2a   = 1.0,                    ! Number fraction allocated to subrange 2a (b-bins will get 1-nf2a)
225
226!
227!-- Aerosol emissions:
228    salsa_emission_mode = 'read_from_file', ! 'no_emission','uniform' or 'read_from_file'
229!-- NOTE! chemical components of the source have to be activated in 'listspec'
230!--       and have to be in the same order
231
232!
233!-- If isdtyp = 0, set the chemical composition of the initial particle size distribution
234    mass_fracs_a = 0.8, 0.1, 0.1, 0.0, 0.0, 0.0, 0.0, ! mass fractions of chemical components: soluble
235    mass_fracs_b = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, ! mass fractions of chemical components: insoluble
236!-- NOTE! Set everyhing to zero in mass_fracs_b if you do not want include insoluble species (default)
237
238!
239!-- Boundary conditions
240!-- Decycle at the left & right boundaries: initial concentration 
241!--   is copied to ghost and first three layers at both boundaries at the left boundary
242!--   and a zero gradient is set at the right boundary
243!-- Decycle at the north & south boundaries: initial concentration 
244!--   is copied to ghost and first three layers at both boundaries
245    decycle_salsa_lr     = .T.,         
246    decycle_salsa_ns     = .F.,
247    decycle_method_salsa = 'dirichlet','dirichlet','dirichlet','dirichlet', 
248    bc_salsa_b     = 'neumann', !'dirichlet',  ! surface flux requires 'neumann'
249    bc_salsa_t     = 'dirichlet',  ! top
250                                                                                           
251!
252!-- Switches for aerosol processes:
253    nldistupdate = .T., ! update aerosol size distribution
254    nldepo       = .T., ! Deposition master switch
255    nldepo_pcm   = .T., ! Deposition on vegetation
256    nldepo_surf  = .T., ! Deposition on walls
257    nlcnd        = .T., ! Condensation master switch
258    nlcndgas     = .T., ! Condensation of precursor gases
259    nlcndh2oae   = .F., ! Condensation of H2O on aerosols                         
260    nlcoag       = .T., ! Coagulation master switch
261    nsnucl       =  0,  ! Nucleation scheme (0 = off)
262    nj3          =  1,  ! J3 parametrization for nucleation
263!
264!-- Deposition on vegetation (pcm = plant canopy model)
265    depo_pcm_par = 'zhang2001',  ! or 'petroff2010'
266!
267!-- Deposition on on ground, walls and roofs
268    depo_surf_par = 'zhang2001', ! or 'petroff2010'
269
270!
271!-- Other switches:
272    advect_particle_water   = .T.,   ! particle water: advect or calculate at each dt_salsa
273    feedback_to_palm        = .F.,   ! feedback to flow due to condensation of water
274    nest_salsa              = .F.,   ! apply nesting for salsa
275    read_restart_data_salsa = .F.,   ! skip reading restart data even if it's done for the flow
276    write_binary_salsa      = .F.,   ! skip writing restart data even if it's done for the flow
277         
278 / ! end of salsa_par namelist
279
280 
Note: See TracBrowser for help on using the repository browser.