Low PM2.5 concentration using WRF-Chem V3.9.1

Post Reply
ganesh_Mk
Posts: 8
Joined: Fri Mar 22, 2019 4:37 am

Low PM2.5 concentration using WRF-Chem V3.9.1

Post by ganesh_Mk » Tue Oct 22, 2019 5:22 am

Dear Users,

I am conducting experimental simulations to study the PM 2.5 concentration.I have tried several chemistry options ranging from RADM2 to MOZCART aerosol chemistry. The simulations were successful but I am getting very low values of PM 2.5 of the order of 0.1 to 0.9 ugm-3 which should be of the order 100 to 900 ugm-3. I am confused as to where I am making mistake. has anyone come across similar kind of issue with PM 2.5 concentration?. what should be my next step in resolving this issue?. My namelist file details are below here for your reference.

&time_control
run_days = 7,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2018, 2018, 2018,
start_month = 10, 10, 10,
start_day = 20, 20, 20,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2018, 2018, 2018
end_month = 10, 10, 10,
end_day = 27, 27, 27,
end_hour = 00, 00, 00,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600,
input_from_file = .true.,.true.,.true.,
history_interval = 60, 60, 60,
frames_per_outfile = 1000, 1000, 1000,
restart = .false.,
restart_interval = 5000,
auxinput5_inname = 'wrfchemi_d<domain>_<date>',
auxinput6_inname = 'wrfbiochemi_d<domain>',
auxinput7_inname = 'wrffirechemi_d<domain>_<date>',
auxinput5_interval_m = 360, 360, 360,
auxinput7_interval_m = 360, 360, 360,
frames_per_auxinput5 = 1, 1, 1,
frames_per_auxinput6 = 1, 1, 1,
frames_per_auxinput7 = 1, 1, 1,
io_form_auxinput5 = 2,
io_form_auxinput6 = 2,
io_form_auxinput7 = 2,
io_form_auxinput12 = 0,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,
/

&domains
time_step = 54,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 100, 73, 163,
e_sn = 79, 55, 145,
e_vert = 30, 30, 30,
dx = 9000, 3000, 3000,
dy = 9000, 3000, 3000,
p_top_requested = 5000,
num_metgrid_levels = 32,
num_metgrid_soil_levels = 4,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 57, 99,
j_parent_start = 1, 15, 71,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
/

&physics
mp_physics = 2, 2, 2,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 2, 2, 2,
radt = 30, 30, 30,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 3, 3, 3,
isfflx = 1,
ifsnow = 1,
mp_zero_out = 2,
progn = 1, 1, 1,
sf_urban_physics = 0, 0, 0,
cu_diag = 1, 1, 1,
cu_rad_feedback = .true.,
/

&fdda
/

&dynamics
w_damping = 1,
diff_opt = 1, 1, 1,
km_opt = 4, 4, 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
tke_adv_opt = 2, 2, 2,
moist_adv_opt = 2, 2, 2,
scalar_adv_opt = 2, 2, 2,
chem_adv_opt = 2, 2, 2,
/

&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
nested = .false., .true., .true.,
/

&grib2
/

&chem
kemit = 19,
chem_opt = 202, 202, 202,
bioemdt = 30, 30, 30,
photdt = 30, 30, 30,
chemdt = 0, 0, 0,
io_style_emissions = 2,
emiss_inpt_opt = 1, 1, 1,
emiss_opt = 3, 3, 3,
chem_in_opt = 1, 1, 1,
phot_opt = 2, 2, 2,
gas_drydep_opt = 1, 1, 1,
aer_drydep_opt = 1, 1, 1,
bio_emiss_opt = 3, 3, 3,
dustwd_onoff = 1, 1, 1,
ne_area = 400,
dust_opt = 4,
dust_schme = 3,
dmsemis_opt = 0,
seas_opt = 2,
gas_bc_opt = 1, 1, 1,
gas_ic_opt = 1, 1, 1,
aer_bc_opt = 1, 1, 1,
aer_ic_opt = 1, 1, 1,
gaschem_onoff = 1, 1, 1,
aerchem_onoff = 1, 1, 1,
wetscav_onoff = 0, 0, 0,
cldchem_onoff = 0, 0, 0,
vertmix_onoff = 1, 1, 1,
chem_conv_tr = 1, 1, 1,
conv_tr_wetscav = 1, 1, 1,
conv_tr_aqchem = 1, 1, 1,
biomass_burn_opt = 2, 2, 3,
plumerisefire_frq = 60, 60, 60,
aer_ra_feedback = 1, 1, 1,
chemdiag = 1, 1, 1,
opt_pars_out = 1,
have_bcs_chem = .false., .false., .false.,
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

Kindly help me in this regard. Thanks for your time and support.

Post Reply

Return to “Aerosols”