[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 | |
---|
[4012] | 39 | day_of_year_init = 65, ! march 06, to avoid extreme temperatures |
---|
| 40 | time_utc_init = 36000.0, ! (0.0 s = 00 UTC midnight), |
---|
[3924] | 41 | |
---|
| 42 | ! |
---|
[4012] | 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., |
---|
[3924] | 50 | |
---|
| 51 | ! |
---|
| 52 | !-- Topography |
---|
| 53 | !------------------------------------------------------------------------------ |
---|
| 54 | topography = 'read_from_file', |
---|
| 55 | |
---|
| 56 | ! |
---|
| 57 | !-- Physics |
---|
| 58 | !------------------------------------------------------------------------------ |
---|
[4012] | 59 | longitude = 13.4, |
---|
| 60 | latitude = 52.5, |
---|
| 61 | |
---|
[3924] | 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 | ! |
---|
[4012] | 78 | !-- Run-control/timeseries output settings |
---|
[3924] | 79 | !------------------------------------------------------------------------------ |
---|
[4012] | 80 | dt_run_control = 0.0, |
---|
| 81 | dt_dots = 10.0, |
---|
[3924] | 82 | |
---|
| 83 | ! |
---|
[4012] | 84 | !-- Profile output settings |
---|
[3924] | 85 | !------------------------------------------------------------------------------ |
---|
[4012] | 86 | skip_time_dopr = 0.0, |
---|
| 87 | dt_dopr = 60.0, |
---|
| 88 | averaging_interval_pr = 60.0, |
---|
| 89 | dt_averaging_input_pr = 0.0, |
---|
[3924] | 90 | |
---|
[4012] | 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', |
---|
[4131] | 96 | 'salsa_Ntot', 'salsa_N_UFP', 'salsa_LDSA', |
---|
| 97 | 'salsa_PM0.1', 'salsa_PM2.5', 'salsa_PM10', |
---|
[3924] | 98 | ! |
---|
| 99 | !-- 2D/3D output settings |
---|
| 100 | !------------------------------------------------------------------------------ |
---|
| 101 | do3d_at_begin = .T., |
---|
[4012] | 102 | do2d_at_begin = .T., |
---|
[3924] | 103 | |
---|
[4012] | 104 | dt_data_output = 10.0, |
---|
| 105 | dt_data_output_av = 10.0, |
---|
| 106 | averaging_interval = 10.0, |
---|
| 107 | dt_averaging_input = 0.0, |
---|
[3924] | 108 | |
---|
[4012] | 109 | section_xy = 0, |
---|
[3924] | 110 | |
---|
[4118] | 111 | ! data_output = 'u', 'v', 'w', |
---|
| 112 | ! 'Ntot', 'PM2.5', 'LDSA', |
---|
| 113 | ! 's_OC','g_OCNV','g_OCSV', |
---|
| 114 | ! 'N_bin3', 'm_bin4', |
---|
[3924] | 115 | |
---|
[4012] | 116 | |
---|
[3924] | 117 | / ! end of d3par namelist |
---|
| 118 | |
---|
| 119 | |
---|
| 120 | !------------------------------------------------------------------------------- |
---|
[4012] | 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 | !------------------------------------------------------------------------------- |
---|
[3924] | 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', |
---|
[4012] | 181 | canopy_drag_coeff = 0.3, |
---|
[3924] | 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 |
---|
[4117] | 247 | decycle_salsa_lr = .T., |
---|
| 248 | decycle_salsa_ns = .F., |
---|
| 249 | decycle_method_salsa = 'dirichlet','dirichlet','dirichlet','dirichlet', |
---|
[3924] | 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 | |
---|