[4409] | 1 | !------------------------------------------------------------------------------ |
---|
| 2 | !-- INITIALIZATION PARAMETER NAMELIST |
---|
| 3 | !-- Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar |
---|
| 4 | !------------------------------------------------------------------------------ |
---|
| 5 | &initialization_parameters |
---|
| 6 | ! |
---|
| 7 | !-- model grid |
---|
| 8 | !-------------------------------------- |
---|
| 9 | nx = 35, |
---|
| 10 | ny = 35, |
---|
| 11 | nz = 80, |
---|
| 12 | |
---|
| 13 | dx = 10.0, |
---|
| 14 | dy = 10.0, |
---|
| 15 | dz = 10.0, |
---|
| 16 | |
---|
| 17 | dz_stretch_level = 1200.0, dz_stretch_factor = 1.02, dz_max = 60.0, |
---|
| 18 | |
---|
| 19 | |
---|
| 20 | !-- pressure solver |
---|
| 21 | !-------------------------------------- |
---|
| 22 | fft_method = 'fftw', |
---|
| 23 | |
---|
| 24 | ! psolver = 'multigrid', |
---|
| 25 | ! cycle_mg = 'w', |
---|
| 26 | ! mg_cycles = 2, |
---|
| 27 | ! mg_switch_to_pe0_level = -1, |
---|
| 28 | |
---|
| 29 | ! |
---|
| 30 | !-- initial and boundary conditions |
---|
| 31 | !-------------------------------------- |
---|
| 32 | initializing_actions = 'set_constant_profiles', ! options: 'set_constant_profiles', or 'inifor', |
---|
| 33 | ! or both, e.g. 'inifor set_constant_profiles' |
---|
| 34 | constant_flux_layer = .TRUE., |
---|
| 35 | ! |
---|
| 36 | !-- wind speed |
---|
| 37 | !-------------------------------------- |
---|
| 38 | ug_surface = 1.0, |
---|
| 39 | vg_surface = 0.0, |
---|
| 40 | ! |
---|
| 41 | !-- temperature |
---|
| 42 | !-------------------------------------- |
---|
| 43 | pt_surface = 290.0, |
---|
| 44 | bc_pt_b = 'dirichlet', |
---|
| 45 | pt_vertical_gradient = 0.1, 1.0, !pt gradient/100m K |
---|
| 46 | pt_vertical_gradient_level= 0.0, 500.0, !height (m) at which pt gradient is effective |
---|
| 47 | ! |
---|
| 48 | !-- humidity |
---|
| 49 | !-------------------------------------- |
---|
| 50 | humidity = .TRUE., |
---|
| 51 | q_surface = 0.001, |
---|
| 52 | bc_q_b = 'dirichlet', |
---|
| 53 | ! |
---|
| 54 | !-- time |
---|
| 55 | !-------------------------------------- |
---|
| 56 | origin_date_time = '2017-07-31 03:00:00 +00' |
---|
| 57 | ! |
---|
| 58 | !-- wall/soil spinup |
---|
| 59 | !-------------------------------------- |
---|
| 60 | spinup_time = 86400.0, ! 24h wall/soil spinup |
---|
| 61 | spinup_pt_mean = 290.0, |
---|
| 62 | spinup_pt_amplitude = 20.0, |
---|
| 63 | dt_spinup = 120.0, ! time step during spinup |
---|
| 64 | data_output_during_spinup = .FALSE., ! set .TRUE. if you want output |
---|
| 65 | ! |
---|
| 66 | !-- building topology |
---|
| 67 | !-------------------------------------- |
---|
| 68 | topography = 'read_from_file', |
---|
| 69 | |
---|
| 70 | / |
---|
| 71 | !------------------------------------------------------------------------------ |
---|
| 72 | !-- RUNTIME PARAMETER NAMELIST |
---|
| 73 | !-- Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par |
---|
| 74 | !------------------------------------------------------------------------------ |
---|
| 75 | &runtime_parameters |
---|
| 76 | ! |
---|
| 77 | !-- run time |
---|
| 78 | !-------------------------------------- |
---|
| 79 | end_time = 7200., |
---|
| 80 | termination_time_needed = 600.0, |
---|
| 81 | ! |
---|
| 82 | !-- processor topology |
---|
| 83 | !-------------------------------------- |
---|
| 84 | ! npex = 6, |
---|
| 85 | ! npey = 6, |
---|
| 86 | ! |
---|
| 87 | !-- 2d/3d output |
---|
| 88 | !-------------------------------------- |
---|
| 89 | do3d_at_begin = .TRUE., |
---|
| 90 | do2d_at_begin = .TRUE., |
---|
| 91 | dt_data_output = 300.0, |
---|
| 92 | dt_data_output_av = 600.0, |
---|
| 93 | dt_do2d_xy = 600.0, |
---|
| 94 | section_xy = 0, |
---|
| 95 | averaging_interval = 600.0, |
---|
| 96 | |
---|
| 97 | |
---|
| 98 | data_output = 'theta', 'q', 'u', 'v', 'w','e', 't_soil', 'm_soil', 'tsurf*_xy', 'shf*_xy', 'qsws*_xy', 'ghf*_xy', 'rad_net*_xy', 'us*_xy', 't*_xy', 'r_a*_xy', 'm_liq*_xy', 'kc_NO2','kc_O3','kc_NO', 'kc_PM10', 'theta_av', 'q_av', 'u_av', 'v_av', 'w_av','e_av', 'tsurf*_xy_av', 'shf*_xy_av', 'qsws*_xy_av', 'ghf*_xy_av', 'kc_NO2_av','kc_O3_av','kc_NO_av', 'kc_PM10_av', |
---|
| 99 | ! |
---|
| 100 | ! Masked output - here for first layer above ground |
---|
| 101 | mask_k_over_surface(1,:) = 0, ! 1. layer above surface (terrain following) |
---|
| 102 | dt_domask = 60., ! time interval for masked output |
---|
| 103 | data_output_masks(1,:) = 'u', 'v', 'kc_NO', 'kc_NO2', 'kc_O3', 'kc_PM10', |
---|
| 104 | ! |
---|
| 105 | !-- timeseries output |
---|
| 106 | !-------------------------------------- |
---|
| 107 | dt_run_control = 0.0, |
---|
| 108 | dt_dots = 60.0, |
---|
| 109 | ! |
---|
| 110 | !-- profile output |
---|
| 111 | !-------------------------------------- |
---|
| 112 | skip_time_dopr = 0.0, |
---|
| 113 | dt_dopr = 600.0, |
---|
| 114 | averaging_interval_pr = 600.0, |
---|
| 115 | dt_averaging_input_pr = 600.0, |
---|
| 116 | |
---|
| 117 | data_output_pr = 'e', 'e*', '#u', 'w"u"', 'w*u*', 'wu', 'u*2', '#v', 'w"v"', 'w*v*', 'wv', 'v*2', 'w', 'w*2', '#theta', 'w"theta"', 'w*theta*','wtheta', 'theta*2', 'kc_NO2','kc_O3','kc_NO', |
---|
| 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 | radiation_scheme = 'clear-sky', !'clear-sky', 'rrtmg', |
---|
| 126 | dt_radiation = 60.0, |
---|
| 127 | |
---|
| 128 | albedo_type = 5, |
---|
| 129 | constant_albedo = .FALSE., |
---|
| 130 | |
---|
| 131 | surface_reflections = .TRUE., |
---|
| 132 | nrefsteps = 5, |
---|
| 133 | |
---|
| 134 | rad_angular_discretization = .FALSE., |
---|
| 135 | raytrace_mpi_rma = .FALSE., |
---|
| 136 | |
---|
| 137 | / ! end of radiation_parameters namelist |
---|
| 138 | |
---|
| 139 | !------------------------------------------------------------------------------ |
---|
| 140 | !-- URBAN SURFACE MODEL PARAMETER NAMELIST |
---|
| 141 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/usmpar |
---|
| 142 | !------------------------------------------------------------------------------ |
---|
| 143 | &urban_surface_parameters |
---|
| 144 | |
---|
| 145 | usm_material_model = .TRUE., |
---|
| 146 | / ! end of urban_surface_parameters namelist |
---|
| 147 | |
---|
| 148 | |
---|
| 149 | !------------------------------------------------------------------------------ |
---|
| 150 | !-- LAND SURFACE MODEL PARAMETER NAMELIST |
---|
| 151 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/lsmpar |
---|
| 152 | !------------------------------------------------------------------------------ |
---|
| 153 | &land_surface_parameters |
---|
| 154 | vegetation_type = 2, |
---|
| 155 | soil_type = 3, |
---|
| 156 | |
---|
| 157 | root_fraction = 0.49, 0.26, 0.25, 0.0, 0.0, 0.0, 0.0, 0.0, |
---|
| 158 | soil_moisture = 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, |
---|
| 159 | soil_temperature = 289., 289., 291., 292., 292., 290., 288., 287., |
---|
| 160 | deep_soil_temperature = 285., |
---|
| 161 | |
---|
| 162 | aero_resist_kray = .TRUE., |
---|
| 163 | constant_roughness = .TRUE., |
---|
| 164 | conserve_water_content = .TRUE., |
---|
| 165 | / ! end of land_surface_parameters namelist |
---|
| 166 | |
---|
| 167 | |
---|
| 168 | !------------------------------------------------------------------------------ |
---|
| 169 | !-- PLANT CANOPY MODEL PARAMETER NAMELIST |
---|
| 170 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/canpar |
---|
| 171 | !------------------------------------------------------------------------------ |
---|
| 172 | &plant_canopy_parameters |
---|
| 173 | canopy_mode = 'read_from_file', |
---|
| 174 | canopy_drag_coeff = 0.3, |
---|
| 175 | / ! end of plant_canopy_parameters namelist |
---|
| 176 | |
---|
| 177 | !------------------------------------------------------------------------------ |
---|
| 178 | !-- CHEMISTRY MODEL PARAMETER NAMELIST |
---|
| 179 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/chempar |
---|
| 180 | !------------------------------------------------------------------------------ |
---|
| 181 | &chemistry_parameters |
---|
| 182 | |
---|
| 183 | chem_gasphase_on = .TRUE., |
---|
| 184 | emissions_anthropogenic = .TRUE., |
---|
| 185 | emiss_lod = 0, ! 0 "PARAMETERIZED", 1 "DEFAULT", 2 "PRE-PROCESSED", |
---|
| 186 | chem_mechanism = "phstatp", |
---|
| 187 | deposition_dry = .TRUE., |
---|
| 188 | photolysis_scheme = 'simple', |
---|
| 189 | |
---|
| 190 | surface_csflux_name = 'NO', 'NO2', 'PM10', !"cs" stands for chemical species, NOT to be used when running |
---|
| 191 | ! with mode_emis = "PRE-PROCESSED" or mode_emis = "DEFAULT" |
---|
| 192 | surface_csflux = 4745.0, 1326.0, 2700e-09, ! surface emission of NO, NO2 in umol m-2 d-1 (PMs in kg m-2 d-1), |
---|
| 193 | ! NOT to be used when running with mode_emis = "PRE-PROCESSED" or mode_emis = "DEFAULT" |
---|
| 194 | emiss_factor_main = 1.667, 1.667, 1.667, ! emiss_factor_main*surface_csflux for emission on main streets, one factor per surface_csflux component needed! |
---|
| 195 | emiss_factor_side = 0.334, 0.334, 0.334, ! emiss_factor_side*surface_csflux for emission on side streets, one factor per surface_csflux component needed! |
---|
| 196 | main_street_id = 11, ! id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90) |
---|
| 197 | side_street_id = 8, |
---|
| 198 | max_street_id = 19, |
---|
| 199 | |
---|
| 200 | cs_name = 'O3', 'NO', 'NO2', 'PM10', |
---|
| 201 | cs_surface = 0.020, 0.000, 0.010, 2.0e-09, !units- chem spcs in ppm, and PM10 in ug/m**2 |
---|
| 202 | |
---|
| 203 | cs_profile(1,:) = 0.015, 0.018, 0.020, 0.021, 0.022, 0.023, 0.024, 0.025, 0.027, 0.030, 0.050, !units for gasphas ppm |
---|
| 204 | cs_profile(2,:) = 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, !units for gasphas ppm |
---|
| 205 | cs_profile(3,:) = 0.010, 0.008, 0.006, 0.005, 0.004, 0.003, 0.002, 0.002, 0.002, 0.002, 0.002, |
---|
| 206 | cs_profile(4,:) = 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 1.0e-09, |
---|
| 207 | |
---|
| 208 | |
---|
| 209 | cs_heights(1,:) = 0.0, 15.0, 25.0, 45.0, 65.0, 85.0, 105.0, 125.0, 155.0, 185.0, 295.0, |
---|
| 210 | cs_heights(2,:) = 0.0, 15.0, 25.0, 45.0, 65.0, 85.0, 105.0, 125.0, 155.0, 185.0, 295.0, |
---|
| 211 | cs_heights(3,:) = 0.0, 15.0, 25.0, 45.0, 65.0, 85.0, 105.0, 125.0, 155.0, 185.0, 295.0, |
---|
| 212 | cs_heights(4,:) = 0.0, 15.0, 15.0, 45.0, 65.0, 85.0, 105.0, 125.0, 155.0, 185.0, 295.0, |
---|
| 213 | |
---|
| 214 | bc_cs_b = 'neumann', !'dirichlet' |
---|
| 215 | bc_cs_t = 'neumann', |
---|
| 216 | |
---|
| 217 | |
---|
| 218 | icntrl(3) = 1, !solver ros2 |
---|
| 219 | icntrl(4) = 500, !max number of chem-substeps |
---|
| 220 | rcntrl(3) = 0.1, ! Hstart, starting value for the integration step size |
---|
| 221 | |
---|
| 222 | / ! end of chemistry_parameters namelist |
---|