[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', |
---|
[4478] | 22 | monotonic_limiter_z = .T., |
---|
[3924] | 23 | ! |
---|
| 24 | !-- Mode |
---|
[4478] | 25 | !------------------------------------------------------------------------------ |
---|
[3924] | 26 | humidity = .T., |
---|
[4012] | 27 | q_surface = 0.001, |
---|
[3924] | 28 | |
---|
| 29 | ! |
---|
| 30 | !-- Initialization |
---|
| 31 | !------------------------------------------------------------------------------ |
---|
[4478] | 32 | initializing_actions = 'inifor set_constant_profiles', !'set_constant_profiles', or 'inifor', |
---|
[4256] | 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 |
---|
[4478] | 43 | !------------------------------------------------------------------------------ |
---|
[4012] | 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 |
---|
[4478] | 52 | !------------------------------------------------------------------------------ |
---|
[3924] | 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', |
---|
[4478] | 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 | !------------------------------------------------------------------------------ |
---|
[4478] | 102 | do3d_at_begin = .T., |
---|
[4012] | 103 | do2d_at_begin = .T., |
---|
[3924] | 104 | |
---|
[4478] | 105 | dt_data_output = 60.0, |
---|
[4256] | 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', |
---|
[4508] | 113 | ! 'salsa_Ntot', 'salsa_PM2.5', 'salsa_LDSA', |
---|
| 114 | ! 'salsa_s_OC','salsa_g_OCNV','salsa_g_OCSV', |
---|
| 115 | ! 'salsa_N_bin3', 'salsa_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 |
---|
[4478] | 145 | !------------------------------------------------------------------------------- |
---|
| 146 | &land_surface_parameters |
---|
[4012] | 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, |
---|
[4478] | 156 | soil_temperature = 277.15, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0, |
---|
[4012] | 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, |
---|
[4478] | 159 | |
---|
[4012] | 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 | |
---|
[4478] | 174 | |
---|
[4012] | 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 |
---|
[4478] | 180 | |
---|
[4342] | 181 | canopy_mode = 'read_from_file', |
---|
[4012] | 182 | canopy_drag_coeff = 0.3, |
---|
[4478] | 183 | |
---|
[3924] | 184 | / ! end of canopy_par namelist |
---|
| 185 | |
---|
[4478] | 186 | |
---|
[3924] | 187 | !------------------------------------------------------------------------------- |
---|
| 188 | !-- SALSA PARAMETER NAMELIST |
---|
| 189 | !------------------------------------------------------------------------------- |
---|
[4478] | 190 | &salsa_parameters |
---|
[3924] | 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 | ! |
---|
[4478] | 198 | !-- If initializing_actions includes 'set_constant_profiles' , |
---|
[4270] | 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) |
---|
[4478] | 201 | sigmag = 1.8, 2.16, 2.21, ! standard deviation per mode |
---|
[4270] | 202 | n_lognorm = 1.43E9, 4.45E8, 7.0E4, ! number concentration per mode (#/m3) |
---|
[3924] | 203 | |
---|
| 204 | ! |
---|
[4478] | 205 | !-- If initializing_actions includes 'set_constant_profiles', apply these initial gas concentrations |
---|
[4270] | 206 | H2SO4_init = 5.0E12, ! sulphuric acid (#/m3) |
---|
| 207 | HNO3_init = 3.0E12, ! nitric acid (#/m3) |
---|
[4478] | 208 | NH3_init = 6.0E12, ! ammonia (#/m3) |
---|
[4270] | 209 | OCNV_init = 1.0E12, ! non-volatile organic gases (#/m3) |
---|
| 210 | OCSV_init = 1.0E12, ! non-volatile organic gases (#/m3) |
---|
[3924] | 211 | |
---|
| 212 | ! |
---|
[4478] | 213 | !-- List of activated chemical components: |
---|
[3924] | 214 | !-- NOTE! Chemical species have to be activated here even if they are not initially present! |
---|
[4478] | 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) |
---|
[4478] | 221 | nbin = 1, 7, ! number of bins per subrange |
---|
[3924] | 222 | !-- NOTE! Subrange 1 consists only of H2SO4 and/or OC |
---|
[4478] | 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' |
---|
[4478] | 228 | !-- NOTE! chemical components of the source have to be activated in 'listspec' |
---|
[3924] | 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 |
---|
[4508] | 239 | !-- Use Dirichlet boundary conditions (initial profiles) at left and right boundaries |
---|
| 240 | !-- Use cyclic boundary conditions along y |
---|
| 241 | bc_aer_l = 'dirichlet', |
---|
| 242 | bc_aer_r = 'dirichlet', |
---|
| 243 | bc_aer_s = 'cyclic', |
---|
| 244 | bc_aer_n = 'cyclic', |
---|
| 245 | bc_aer_b = 'neumann', ! surface flux requires 'neumann' |
---|
| 246 | bc_aer_t = 'dirichlet', ! top |
---|
[4478] | 247 | |
---|
[3924] | 248 | ! |
---|
| 249 | !-- Switches for aerosol processes: |
---|
[4270] | 250 | nldistupdate = .T., ! Update size distribution (default .T.) |
---|
| 251 | nldepo = .T., ! Deposition master switch |
---|
| 252 | nldepo_pcm = .T., ! Deposition on vegetation |
---|
| 253 | nldepo_surf = .T., ! Deposition on walls |
---|
| 254 | nlcnd = .F., ! Condensation master switch |
---|
| 255 | nlcndgas = .F., ! Condensation of precursor gases |
---|
| 256 | nlcndh2oae = .F., ! Condensation of H2O on aerosols |
---|
| 257 | nlcoag = .F., ! Coagulation master switch |
---|
| 258 | nsnucl = 0, ! Nucleation scheme (0 = off) |
---|
[4478] | 259 | nj3 = 1, ! J3 parametrization for nucleation |
---|
[4270] | 260 | |
---|
[3924] | 261 | ! |
---|
[4270] | 262 | !-- Deposition parametrisations: |
---|
| 263 | depo_pcm_par = 'zhang2001', ! plant canopy |
---|
| 264 | depo_surf_par = 'zhang2001', ! ground, roofs and walls |
---|
[4478] | 265 | season_z01 = 5 ! Season for zhang2001, |
---|
[4270] | 266 | ! 1 = summer, 2 = autumn, 3 = late autumn |
---|
| 267 | ! 4 = winter, 5 = transitional spring |
---|
[3924] | 268 | ! |
---|
| 269 | !-- Other switches: |
---|
[4478] | 270 | advect_particle_water = .T., ! particle water: advect or calculate at each dt_salsa |
---|
[3924] | 271 | feedback_to_palm = .F., ! feedback to flow due to condensation of water |
---|
[4273] | 272 | nesting_salsa = .F., ! apply self-nesting for salsa variables |
---|
[4270] | 273 | nesting_offline_salsa = .F., ! apply offline nesting for salsa |
---|
[4478] | 274 | |
---|
[3924] | 275 | / ! end of salsa_par namelist |
---|