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

Last change on this file since 4270 was 4270, checked in by monakurppa, 4 years ago

Updates for the aerosol module salsa

  • Reformat the emission data for the pre-processed mode (lod=2) to limit the size of the input file. Now, the emission is given as total number emission per emission category (ncat) and the number size distribution is given in variable emission_number_fracs(ncat,nbins).
  • Implement offline nesting for salsa. Additionally, offline nesting can be turn off specifically for salsa using the switch nesting_offline_salsa.
  • Remove init_aerosol_type and init_gases_type from salsa_parin and define them based on the initializing_actions
  • parameter definition removed from "season" and "season_z01" is added to parin
  • bugfix in application of index_hh after implementing the new palm_date_time_mod
File size: 10.8 KB
RevLine 
[3924]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,
[4256]12    nz                   = 60,
[3924]13
14    dx                   = 2.0,
15    dy                   = 2.0,
16    dz                   = 2.0,
17
18!
19!-- Numerics
20!------------------------------------------------------------------------------
[4012]21    fft_method           = 'fftw',
[4117]22    monotonic_limiter_z  = .T.,
[3924]23!
24!-- Mode
25!------------------------------------------------------------------------------     
26    humidity             = .T.,
[4012]27    q_surface            = 0.001,
[3924]28
29!
30!-- Initialization
31!------------------------------------------------------------------------------
[4256]32    initializing_actions = 'inifor set_constant_profiles',   !'set_constant_profiles', or 'inifor',
33                                                             ! or both (met. from inifor, chem constant profiles)
[3924]34
[4012]35    pt_surface           = 277.15,
[3924]36    ug_surface           = 1.0,
37    vg_surface           = 0.0,
38
[4227]39    origin_date_time     = '2019-03-06 10:00:00 +00'
[3924]40
41!
[4012]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.,
[3924]49
50!
51!-- Topography
52!------------------------------------------------------------------------------   
53    topography           = 'read_from_file',
54
55!
56!-- Physics
57!------------------------------------------------------------------------------
[4012]58    longitude            = 13.4,
59    latitude             = 52.5,
60
[3924]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!------------------------------------------------------------------------------
[4256]73    end_time              = 31.0,
[3924]74    create_disturbances   = .T.,
75
76!
[4012]77!-- Run-control/timeseries output settings
[3924]78!------------------------------------------------------------------------------
[4012]79    dt_run_control        =  0.0,
[4256]80    dt_dots               =  5.0,
[3924]81
82!
[4012]83!-- Profile output settings
[3924]84!------------------------------------------------------------------------------
[4012]85    skip_time_dopr        =  0.0,
[4256]86    dt_dopr               = 30.0,
87    averaging_interval_pr = 30.0,
[4012]88    dt_averaging_input_pr =  0.0,
[3924]89
[4012]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"',
[4256]94                            '#q',  'q*2',
95                            'e', 'e*', '#km', '#l',
96                            '#t_soil', '#m_soil', 'rad_lw_in', 'rad_lw_out',
[4131]97                            'salsa_Ntot', 'salsa_N_UFP', 'salsa_LDSA',
98                            'salsa_PM0.1', 'salsa_PM2.5', 'salsa_PM10',
[3924]99!
100!-- 2D/3D output settings
101!------------------------------------------------------------------------------
102    do3d_at_begin      = .T.,
[4012]103    do2d_at_begin      = .T.,
[3924]104
[4256]105    dt_data_output     = 60.0,
106    dt_data_output_av  = 60.0,
107    averaging_interval = 60.0,
[4012]108    dt_averaging_input =  0.0,
[3924]109
[4012]110    section_xy         = 0,
[3924]111
[4118]112!     data_output = 'u', 'v', 'w',
113!                   'Ntot', 'PM2.5', 'LDSA',
114!                   's_OC','g_OCNV','g_OCSV',
115!                   'N_bin3', 'm_bin4',
[3924]116
[4012]117
[3924]118 / ! end of d3par namelist
119
120
121!-------------------------------------------------------------------------------
[4012]122!-- RADIATION MODEL PARAMETER NAMELIST
123!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/radpar
124!-------------------------------------------------------------------------------
125  &radiation_parameters
126
[4256]127     radiation_scheme =  'rrtmg',    !'clear-sky' or 'rrtmg',
[4012]128     albedo_type = 5,
129     constant_albedo = .F.,
130
131     dt_radiation = 60.0,
132
133     surface_reflections = .T.,
134     nrefsteps = 3,
135
136     rad_angular_discretization = .F.,
137     raytrace_mpi_rma = .F.,
138
139  /
140
141
142!-------------------------------------------------------------------------------
143!-- LAND SURFACE MODEL PARAMETER NAMELIST
144!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/lsmpar
145!-------------------------------------------------------------------------------
146 &land_surface_parameters
147
148    constant_roughness     = .T.,
149    aero_resist_kray       = .T.,
150
151    vegetation_type        = 2,
152    soil_type              = 3,
153    conserve_water_content = .T.,
154
155    root_fraction          = 0.49, 0.26, 0.25, 0.0, 0.0, 0.0, 0.0, 0.0,
156    soil_temperature       = 277.15, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0,
157    soil_moisture          = 0.30,   0.30,  0.30,  0.30,  0.30, 0.30,  0.30,  0.30,
158    deep_soil_temperature  = 277.0,
159   
160 / ! end of lsm_par namelist
161
162
163!-------------------------------------------------------------------------------
164!-- URBAN SURFACE MODEL PARAMETER NAMELIST
165!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/usmpar
166!-------------------------------------------------------------------------------
167 &urban_surface_parameters
168
169    usm_material_model     = .T.,
170    usm_wall_mod = .T.,
171
172 / ! end of urban_surface_par namelist
173 
174
175!-------------------------------------------------------------------------------
[3924]176!-- PLANT CANOPY MODEL PARAMETER NAMELIST
177!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/canpar
178!-------------------------------------------------------------------------------
179 &plant_canopy_parameters
180 
181    canopy_mode       = 'read_from_file_3d',
[4012]182    canopy_drag_coeff = 0.3,
[3924]183     
184 / ! end of canopy_par namelist
[4256]185 
[3924]186
187!-------------------------------------------------------------------------------
188!-- SALSA PARAMETER NAMELIST
189!-------------------------------------------------------------------------------
190 &salsa_parameters 
191
192!
193!-- Time stepping
[4270]194    dt_salsa             =   2.0, ! time step for calculating aerosol processes (s)
195    skip_time_do_salsa   =  10.0, ! starting time of SALSA (s)
[3924]196
197!
[4270]198!-- If initializing_actions includes 'set_constant_profiles' ,
199!-- define the initial aerosol size distribution by dpg, sigmag and n_lognorm
200    dpg                  = 13.5E-9, 54.0E-9, 864.1E-9, ! mean diameter per mode (in metres)
201    sigmag               =     1.8,    2.16,     2.21, ! standard deviation per mode
202    n_lognorm            =  1.43E9,  4.45E8,    7.0E4, ! number concentration per mode (#/m3)
[3924]203
204!
[4270]205!-- If initializing_actions includes 'set_constant_profiles', apply these initial gas concentrations                   
206    H2SO4_init           = 5.0E12, ! sulphuric acid (#/m3)
207    HNO3_init            = 3.0E12, ! nitric acid (#/m3)
208    NH3_init             = 6.0E12, ! ammonia (#/m3)       
209    OCNV_init            = 1.0E12, ! non-volatile organic gases (#/m3)
210    OCSV_init            = 1.0E12, ! non-volatile organic gases (#/m3)
[3924]211
212!
213!-- List of activated chemical components:
214!-- NOTE! Chemical species have to be activated here even if they are not initially present!
[4270]215    listspec             = 'OC','NO','NH','','','','',       
216    ! listspec             = 'SO4','OC','BC','DU','SS','NO','NH',
[3924]217
218!
219!-- Sectional presentation of the particle size distribution
[4270]220    reglim               = 3.0E-9, 1.0E-8, 2.5E-6, ! limits of the subranges (m)
221    nbin                 = 1, 7,                   ! number of bins per subrange
[3924]222!-- NOTE! Subrange 1 consists only of H2SO4 and/or OC
[4270]223    nf2a                 = 1.0,                    ! Number fraction allocated to subrange 2a
[3924]224
225!
226!-- Aerosol emissions:
[4270]227    salsa_emission_mode  = 'read_from_file', ! 'no_emission','uniform' or 'read_from_file'
[3924]228!-- NOTE! chemical components of the source have to be activated in 'listspec'
229!--       and have to be in the same order
230
231!
232!-- If isdtyp = 0, set the chemical composition of the initial particle size distribution
[4270]233    mass_fracs_a         = 0.8, 0.1, 0.1, 0.0, 0.0, 0.0, 0.0, ! mass fractions of soluble components
234    mass_fracs_b         = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, ! mass fractions of insoluble components
[3924]235!-- NOTE! Set everyhing to zero in mass_fracs_b if you do not want include insoluble species (default)
236
237!
238!-- Boundary conditions
239!-- Decycle at the left & right boundaries: initial concentration 
240!--   is copied to ghost and first three layers at both boundaries at the left boundary
241!--   and a zero gradient is set at the right boundary
242!-- Decycle at the north & south boundaries: initial concentration 
243!--   is copied to ghost and first three layers at both boundaries
[4117]244    decycle_salsa_lr     = .T.,         
245    decycle_salsa_ns     = .F.,
246    decycle_method_salsa = 'dirichlet','dirichlet','dirichlet','dirichlet', 
[4270]247    bc_salsa_b           = 'neumann',    ! surface flux requires 'neumann'
248    bc_salsa_t           = 'dirichlet',  ! top
[3924]249                                                                                           
250!
251!-- Switches for aerosol processes:
[4270]252    nldistupdate         = .T., ! Update size distribution (default .T.)
253    nldepo               = .T., ! Deposition master switch
254    nldepo_pcm           = .T., ! Deposition on vegetation
255    nldepo_surf          = .T., ! Deposition on walls
256    nlcnd                = .F., ! Condensation master switch
257    nlcndgas             = .F., ! Condensation of precursor gases
258    nlcndh2oae           = .F., ! Condensation of H2O on aerosols
259    nlcoag               = .F., ! Coagulation master switch
260    nsnucl               =  0,  ! Nucleation scheme (0 = off)
261    nj3                  =  1,  ! J3 parametrization for nucleation
262
[3924]263!
[4270]264!-- Deposition parametrisations:
265    depo_pcm_par         = 'zhang2001',  ! plant canopy
266    depo_surf_par        = 'zhang2001',  ! ground, roofs and walls
267    season_z01           = 5             ! Season for zhang2001,
268                                         ! 1 = summer, 2 = autumn, 3 = late autumn
269                                         ! 4 = winter, 5 = transitional spring
[3924]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
[4270]274    nest_salsa              = .F.,   ! apply self-nesting for salsa variables
[3924]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
[4270]277    nesting_offline_salsa   = .F.,   ! apply offline nesting for salsa
[3924]278         
279 / ! end of salsa_par namelist
280
281 
Note: See TracBrowser for help on using the repository browser.