- Timestamp:
- Nov 2, 2020 10:31:45 AM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/TUTORIALS/cases/chem_phstatp/INPUT/chem_phstatp_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', … … 54 54 !-- time 55 55 !-------------------------------------- 56 origin_date_time = '2017-07-31 0 3:00:00 +00'56 origin_date_time = '2017-07-31 06:00:00 +00' 57 57 ! 58 58 !-- wall/soil spinup … … 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', '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', 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_PM10', 103 'theta_av', 'q_av', 'u_av', 'v_av', 'w_av','e_av', 104 'tsurf*_xy_av', 'shf*_xy_av', 'qsws*_xy_av', 'ghf*_xy_av', 105 'kc_NO2_av','kc_O3_av','kc_NO_av', 'kc_PM10_av', 99 106 ! 100 107 ! Masked output - here for first layer above ground … … 115 122 dt_averaging_input_pr = 600.0, 116 123 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', 124 data_output_pr = 'e', 'e*', 125 '#u', 'w"u"', 'w*u*', 'wu', 'u*2', 126 '#v', 'w"v"', 'w*v*', 'wv', 'v*2', 127 'w', 'w*2', 128 '#theta', 'w"theta"', 'w*theta*','wtheta', 'theta*2', 129 'kc_NO2','kc_O3','kc_NO', 118 130 / 119 131 … … 127 139 128 140 albedo_type = 5, 129 constant_albedo = .F ALSE.,130 131 surface_reflections = .T RUE.,141 constant_albedo = .F., 142 143 surface_reflections = .T., 132 144 nrefsteps = 5, 133 134 rad_angular_discretization = .FALSE.,135 raytrace_mpi_rma = .FALSE.,136 145 137 146 / ! end of radiation_parameters namelist … … 143 152 &urban_surface_parameters 144 153 145 usm_material_model = .TRUE.,146 154 / ! end of urban_surface_parameters namelist 147 155 … … 160 168 deep_soil_temperature = 285., 161 169 162 aero_resist_kray = .T RUE.,163 constant_roughness = .T RUE.,164 conserve_water_content = .T RUE.,170 aero_resist_kray = .T., 171 constant_roughness = .T., 172 conserve_water_content = .T., 165 173 / ! end of land_surface_parameters namelist 166 174 … … 182 190 183 191 chem_gasphase_on = .TRUE., 184 emissions_anthropogenic = .TRUE.,185 emiss_lod = 0, ! 0 "PARAMETERIZED", 1 "DEFAULT", 2 "PRE-PROCESSED",186 192 chem_mechanism = "phstatp", 187 193 deposition_dry = .TRUE., 188 194 photolysis_scheme = 'simple', 189 195 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) 196 ! initial profiles 197 cs_name = 'O3', 'NO', 'NO2', 'PM10', 198 cs_surface = 0.025, 0.000, 0.010, 2.0e-09, !units- chem spcs in ppm, and PM10 in ug/m**2 199 200 cs_profile(1,:) = 0.025, 0.040, 0.050, 201 cs_profile(2,:) = 0.000, 0.000, 202 cs_profile(3,:) = 0.010, 0.002, 203 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, 204 205 cs_heights(1,:) = 0.0, 305.0, 495.0, 206 cs_heights(2,:) = 0.0, 295.0, 207 cs_heights(3,:) = 0.0, 295.0, 208 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, 209 210 ! Emissions 211 emissions_anthropogenic = .TRUE., 212 emiss_lod = 0, ! 0 "PARAMETERIZED", 1 "DEFAULT", 2 "PRE-PROCESSED", 213 ! Parameters for parameterized traffic emissions 214 surface_csflux_name = 'NO', 'NO2', 'PM10', !"cs" stands for chemical species, 215 ! only for parameterized traffic emissions 216 ! NOT to be used when running with emiss_lod>0 217 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),only for parameterized traffic emissions 218 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! 219 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! 220 main_street_id = 11, ! id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90) 197 221 side_street_id = 8, 198 222 max_street_id = 19, 199 223 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' 224 ! Bottom and top boundary conditions (lateral BCs are cyclic, which is the default) 225 bc_cs_b = 'neumann', 215 226 bc_cs_t = 'neumann', 216 227 217 218 icntrl(3) = 1, ! solver ros2219 icntrl(4) = 500, ! max number of chem-substeps228 ! Parameters chemistry solver (see KPP description) 229 icntrl(3) = 1, ! solver ros2 230 icntrl(4) = 500, ! max number of chem-substeps 220 231 rcntrl(3) = 0.1, ! Hstart, starting value for the integration step size 221 232
Note: See TracChangeset
for help on using the changeset viewer.
