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