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.
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,
/
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.
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,
/