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

Last change on this file since 4131 was 4131, checked in by monakurppa, 2 years ago

Several changes in the salsa aerosol module:

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