Changeset 4765 for palm/trunk/TUTORIALS/cases/chem_smog/INPUT/chem_smog_p3d
- Timestamp:
- Nov 2, 2020 10:31:45 AM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/TUTORIALS/cases/chem_smog/INPUT/chem_smog_p3d
r4409 r4765 20 20 !-- pressure solver 21 21 !-------------------------------------- 22 23 24 !psolver = 'multigrid',25 !cycle_mg = 'w',26 !mg_cycles = 2,27 !mg_switch_to_pe0_level = -1,22 ! fft_method = 'fftw', 23 24 psolver = 'multigrid', 25 cycle_mg = 'w', 26 mg_cycles = 2, 27 mg_switch_to_pe0_level = -1, 28 28 29 29 ! … … 32 32 initializing_actions = 'set_constant_profiles', ! options: 'set_constant_profiles', or 'inifor', 33 33 ! or both, e.g. 'inifor set_constant_profiles' 34 constant_flux_layer = .T RUE.,34 constant_flux_layer = .T., 35 35 ! 36 36 !-- wind speed … … 48 48 !-- humidity 49 49 !-------------------------------------- 50 humidity = .T RUE.,50 humidity = .T., 51 51 q_surface = 0.001, 52 52 bc_q_b = 'dirichlet', … … 62 62 spinup_pt_amplitude = 20.0, 63 63 dt_spinup = 120.0, ! time step during spinup 64 data_output_during_spinup = .F ALSE., ! set .TRUE. if you want output64 data_output_during_spinup = .F., ! set .T. if you want output 65 65 ! 66 66 !-- building topology … … 87 87 !-- 2d/3d output 88 88 !-------------------------------------- 89 do3d_at_begin = .T RUE.,90 do2d_at_begin = .T RUE.,89 do3d_at_begin = .T., 90 do2d_at_begin = .T., 91 91 dt_data_output = 300.0, 92 92 dt_data_output_av = 600.0, … … 96 96 97 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', '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_RH','kc_RCHO','kc_CO', 'kc_O', 'kc_RO2', 'kc_OH','kc_HO2','kc_HNO3', 98 data_output = 'theta', 'q', 'u', 'v', 'w','e', 99 't_soil', 'm_soil', 100 'tsurf*_xy', 'shf*_xy', 'qsws*_xy', 'ghf*_xy', 101 'rad_net*_xy', 'us*_xy', 't*_xy', 'r_a*_xy', 'm_liq*_xy', 102 'kc_NO2','kc_O3','kc_NO', 'kc_RH', 'kc_RCHO','kc_CO', 103 'kc_OH','kc_HO2','kc_HNO3', 104 'theta_av', 'q_av', 'u_av', 'v_av', 'w_av','e_av', 105 'tsurf*_xy_av', 'shf*_xy_av', 'qsws*_xy_av', 'ghf*_xy_av', 106 'kc_NO2_av','kc_O3_av','kc_NO_av', 99 107 ! 100 108 ! Masked output - here for first layer above ground … … 115 123 dt_averaging_input_pr = 600.0, 116 124 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', 125 data_output_pr = 'e', 'e*', 126 '#u', 'w"u"', 'w*u*', 'wu', 'u*2', 127 '#v', 'w"v"', 'w*v*', 'wv', 'v*2', 128 'w', 'w*2', 129 '#theta', 'w"theta"', 'w*theta*','wtheta', 'theta*2', 130 'kc_NO2','kc_O3','kc_NO', 118 131 / 119 132 … … 127 140 128 141 albedo_type = 5, 129 constant_albedo = .F ALSE.,130 131 surface_reflections = .T RUE.,142 constant_albedo = .F., 143 144 surface_reflections = .T., 132 145 nrefsteps = 5, 133 134 rad_angular_discretization = .FALSE.,135 raytrace_mpi_rma = .FALSE.,136 146 137 147 / ! end of radiation_parameters namelist … … 143 153 &urban_surface_parameters 144 154 145 usm_material_model = .TRUE.,146 155 / ! end of urban_surface_parameters namelist 147 156 … … 160 169 deep_soil_temperature = 285., 161 170 162 aero_resist_kray = .T RUE.,163 constant_roughness = .T RUE.,164 conserve_water_content = .T RUE.,171 aero_resist_kray = .T., 172 constant_roughness = .T., 173 conserve_water_content = .T., 165 174 / ! end of land_surface_parameters namelist 166 175 … … 182 191 183 192 chem_gasphase_on = .TRUE., 184 emissions_anthropogenic = .TRUE., 185 emiss_lod = 0, ! 0 "PARAMETERIZED", 1 "DEFAULT", 2 "PRE-PROCESSED", 186 chem_mechanism = 'smog', 193 chem_mechanism = "smog", 187 194 deposition_dry = .TRUE., 188 195 photolysis_scheme = 'simple', 189 196 190 surface_csflux_name = 'NO', 'NO2', 'RH', 'RCHO', 'CO', !"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, 649.0, 30.0, 30000.0, ! 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" 197 ! initial profiles 198 cs_name = 'O3', 'NO', 'NO2', 'RH','RCHO', 'CO', 199 cs_surface = 0.025, 0.000, 0.010, 0.010, 0.005, 0.100, !units- chem spcs in ppm, and PM10 in ug/m**2 200 201 cs_profile(1,:) = 0.025, 0.040, 0.050, 202 cs_profile(2,:) = 0.000, 0.000, 203 cs_profile(3,:) = 0.010, 0.002, 204 cs_profile(4,:) = 0.010, 0.005, 205 cs_profile(5,:) = 0.005, 0.0025, 206 cs_profile(6,:) = 0.100, 207 208 cs_heights(1,:) = 0.0, 305.0, 495.0, 209 cs_heights(2,:) = 0.0, 295.0, 210 cs_heights(3,:) = 0.0, 295.0, 211 cs_heights(4,:) = 0.0, 295.0, 212 cs_heights(5,:) = 0.0, 295.0, 213 cs_heights(6,:) = 0.0, 214 215 ! Emissions 216 emissions_anthropogenic = .TRUE., 217 emiss_lod = 0, ! 0 "PARAMETERIZED", 1 "DEFAULT", 2 "PRE-PROCESSED", 218 ! Parameters for parameterized traffic emissions 219 surface_csflux_name = 'NO', 'NO2', 'RH', 'RCHO', 'CO', !"cs" stands for chemical species, 220 ! only for parameterized traffic emissions 221 ! NOT to be used when running with emiss_lod>0 222 surface_csflux = 4745.0, 1326.0, 649., 30., 30000., ! surface emission of NO, NO2 in umol m-2 d-1 (PMs in kg m-2 d-1),only for parameterized traffic emissions 194 223 emiss_factor_main = 1.667, 1.667, 1.667, 1.667, 1.667, ! emiss_factor_main*surface_csflux for emission on main streets, one factor per surface_csflux component needed! 195 224 emiss_factor_side = 0.334, 0.334, 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, 225 main_street_id = 11, ! id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90) 197 226 side_street_id = 8, 198 227 max_street_id = 19, 199 228 200 cs_name = 'O3', 'NO', 'NO2', 'RH', 'RCHO', 'CO', 201 cs_surface = 0.020, 0.000, 0.010, 0.010, 0.005, 0.1, !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,:) = 0.010, 0.005, 207 cs_profile(5,:) = 0.005, 0.0025, 208 cs_profile(6,:) = 0.1, 209 210 211 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, 212 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, 213 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, 214 cs_heights(4,:) = 0.0, 295.0, 215 cs_heights(5,:) = 0.0, 295.0, 216 cs_heights(6,:) = 0.0, 217 218 bc_cs_b = 'neumann', !'dirichlet' 229 ! Bottom and top boundary conditions (lateral BCs are cyclic, which is the default) 230 bc_cs_b = 'neumann', 219 231 bc_cs_t = 'neumann', 220 232 221 222 icntrl(3) = 1, ! solver ros2223 icntrl(4) = 500, ! max number of chem-substeps233 ! Parameters chemistry solver (see KPP description) 234 icntrl(3) = 1, ! solver ros2 235 icntrl(4) = 500, ! max number of chem-substeps 224 236 rcntrl(3) = 0.1, ! Hstart, starting value for the integration step size 225 237
Note: See TracChangeset
for help on using the changeset viewer.
