Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

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

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

ganesh_Mk

New member
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.
 
Top