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

Abnormal nitrate aerosols concentrations

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.

jingliang0811

New member
Deal all,

I am running WRF-Chem v4.2.1 with chem_opt = 11 (RADM2-MADE/SORGAM-SQ) for one domain. The default initial and boundary condition data are generated by ./real.exe. There isn't any emission file for use. However, the model always hangs up after a few hours of simulation. After my check, this error is caused by abnormal simulation of nitrate aerosols and no2 concentrations in the last output, as shown below. These abnormal values just occur suddenly in the last output, like from emission sources. My settings in the namelist.input are also shown below. Thank you in advance for any help and suggestion.


no3ai.png no3aj.png no2.png


namelist.input
&time_control
run_days = 3,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2019, 2016, 2000,
start_month = 06, 01, 01,
start_day = 01, 01, 24,
start_hour = 00, 00, 12,
end_year = 2019, 2016, 2000,
end_month = 06, 01, 01,
end_day = 04, 04, 25,
end_hour = 00, 00, 12,
interval_seconds = 21600,
input_from_file = .true.,.true.,.true.,
history_interval = 60, 60, 60,
frames_per_outfile = 100, 100, 1000,
restart = .false.,
restart_interval = 7200,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 200
/

&domains
time_step = 30,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 130, 55, 94,
e_sn = 115, 73, 91,
e_vert = 33, 33, 33,
p_top_requested = 5000,
num_metgrid_levels = 32,
num_metgrid_soil_levels = 4,
dx = 9000, 4000, 3333.33,
dy = 9000, 4000, 3333.33,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 67, 30,
j_parent_start = 1, 84, 30,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
/

&physics
mp_physics = 2, 2, 2,
progn = 1, 1, 0,
naer = 1e9
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 2, 2, 2,
radt = 27, 9, 10,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
bldt = 1, 1, 0,
cu_physics = 5, 5, 0,
cudt = 1, 1, 1,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 3,
num_land_cat = 21,
sf_urban_physics = 0, 0, 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .true.,.true.,.true.,
cu_diag = 1, 1, 1,
cugd_avedx = 1,
mp_zero_out = 2,
mp_zero_out_thresh = 1e-8
/

&fdda
/

&dynamics
hybrid_opt = 2,
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 = 3,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 2, 2, 2,
scalar_adv_opt = 2, 2, 2,
chem_adv_opt = 2, 2, 2,
tke_adv_opt = 2, 2, 2,
tracer_opt = 0, 0, 0,
/

&bdy_control
spec_bdy_width = 5,
specified = .true., .true., .false.,
/

&grib2
/

&chem
chem_opt = 11, 11,
kemit = 19,
bioemdt = 30, 30,
photdt = 30, 30,
chemdt = 0, 0,
io_style_emissions = 0,
emiss_inpt_opt = 0, 1,
emiss_opt = 0, 3,
chem_in_opt = 0, 0,
phot_opt = 1, 1,
gas_drydep_opt = 1, 1,
aer_drydep_opt = 1, 1,
bio_emiss_opt = 0, 1,
ne_area = 83,
dust_opt = 13,
dustwd_onoff = 0, 0,
aer_op_opt = 1, 1,
dmsemis_opt = 0,
seas_opt = 0,
gas_bc_opt = 1, 1,
gas_ic_opt = 1, 1,
aer_bc_opt = 1, 1,
aer_ic_opt = 1, 1,
gaschem_onoff = 1, 1,
aerchem_onoff = 1, 1,
wetscav_onoff = 1, 1,
cldchem_onoff = 1, 1,
vertmix_onoff = 1, 1,
chem_conv_tr = 1, 1,
biomass_burn_opt = 0, 0,
plumerisefire_frq = 30, 30,
aer_ra_feedback = 1, 1,
have_bcs_chem = .true., .true.,
diagnostic_dep = 1, 1,
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
 

Attachments

  • namelist.input
    6.4 KB · Views: 0
Top