WRF-Chem模拟PM10和PM2.5浓度离谱

燃烧的雪花

各位大佬，我在使用WRF-Chem模拟PM2.5的时候，发现在沙漠地区和沿海地区会出现浓度异常高的情况，甚至能达到10^8。经过我的分析，发现在开着沙尘方案的时候会出现这样的情况，而关闭沙尘模拟就不会有这种情况，以下是我的namelist，请问有知道为什么会出现这种现象以及如何解决的嘛？
&time_control
start_year                          = 2025, 2025, 2020, 2020,
start_month                         = 03,   03,   10,   06
start_day                           = 24,   24,   13,   10
start_hour                          = 00,   00,   00,   00,
start_minute                        = 00,   00,   00,   00,
start_second                        = 00,   00,   00,   00,
end_year                            = 2025, 2025, 2020, 2018,
end_month                           = 03,   03,   10,   06
end_day                             = 28,   28,   16,   16
end_hour                            = 00,   00,   00,   00,
end_minute                          = 00,   00,   00,   00,
end_second                          = 00,   00,   00,   00,
interval_seconds                    = 21600,
input_from_file                     = .true.,.true.,.true.,.true.,
history_interval                    = 60,   60,   60,  60,
frames_per_outfile                  = 1,    1,    1,   1,
history_outname                     = './output/test3/wrfout_d<domain>_<date>'
restart                             = .false.,
restart_interval                    = 1440
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
debug_level                         = 300
iofields_filename = "my_file_d01.txt","my_file_d02.txt","my_file_d03.txt"
auxinput4_interval                  = 360, 360, 360 360,
io_form_auxinput4                   = 2
auxinput4_inname                    = 'wrflowinp_d<domain>',
write_hist_at_0h_rst                = .false.
auxinput5_inname                    = 'wrfchemi_d<domain>_<date>',
frames_per_auxinput5                = 24, 24, 24, 24,
auxinput5_interval_m                = 60, 60, 60, 60,
io_form_auxinput5                   = 2,
auxinput6_inname                    = 'wrfbiochemi_d<domain>',
io_form_auxinput6                   = 2,
!auxinput7_inname                    = 'wrffirechemi_d<domain>_<date>',
!io_form_auxinput7                   = 2,
!frames_per_auxinput7                = 1,  1,  1,  1,
!auxinput7_interval_m                = 60, 60, 60, 60,
/

&domains
time_step                           = 60,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 2
e_we                                = 108, 121,   301, 94,
e_sn                                = 92,  121,   251, 97,
e_vert                              = 28,   28,    28,  28,
p_top_requested                     = 10000,
eta_levels                          = 1,0.993,0.983,0.97,0.954,0.934,0.909,0.880,0.835,0.790,0.745,0.700, 0.621,0.548,0.482,0.420,0.364,0.313,0.266,0.224,0.185,0.145, 0.118,0.089,0.063,0.040,0.020,0
num_metgrid_levels                  = 34,
num_metgrid_soil_levels             = 4,
dx                                  = 25000, 5000, 1000, 1000,
dy                                  = 25000, 5000, 1000, 1000,
grid_id                             = 1,     2,     3,  4,
parent_id                           = 0,     1,     2,  3,
i_parent_start                      = 1,    55,    36, 119,
j_parent_start                      = 1,    48,    22, 123,
parent_grid_ratio                   = 1,     5,     5,  3,
parent_time_step_ratio              = 1,     3,     3,  3,
feedback                            = 1,
smooth_option                       = 0
interp_type                         = 1
lowest_lev_from_sfc                 = .false.
lagrange_order                      = 1
force_sfc_in_vinterp                = 1
zap_close_levels                    = 500
sfcp_to_sfcp                        = .false.
adjust_heights                      = .false.

/

&physics
sst_update                          = 1
mp_physics                          = 2,     2,     2,   2,
ra_lw_physics                       = 5,     5,     5,   5,
ra_sw_physics                       = 2,     2,     2,   2,
radt                                = 27,    27,    27,  27,
sf_sfclay_physics                   = 2,     2,     2,   2,
sf_surface_physics                  = 2,     2,     2,   2,
bl_pbl_physics                      = 2,     2,     2,   2,
bldt                                = 0,     0,     0,   0,
cu_physics                          = 2,     2,     0,   3,
!cu_diag                             = 1,     1,     1,   1,
cudt                                = 5,     5,     5,   5,
isfflx                              = 1,
ifsnow                              = 0,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
num_land_cat                        = 21,
surface_input_source                = 1,
sf_urban_physics                    = 1,    1,     1,    1,
tmn_update                          = 1,
sst_skin                            = 1,
progn                               = 1,    1,     1,    1,
mp_zero_out                         = 0,
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
/
&fdda
grid_fdda                           = 1,    1,     1,   1,
gfdda_inname                        = "wrffdda_d<domain>",
gfdda_end_h                         = 9999,  9999,  9999,  9999,
gfdda_interval_m                    = 360,   360,   360,   360,
fgdt                                = 0,     0,     0,     0,
if_no_pbl_nudging_uv                = 0,     0,     0,     0,
if_no_pbl_nudging_t                 = 0,     0,     0,     0,
if_no_pbl_nudging_q                 = 0,     0,     0,     0,
if_zfac_uv                          = 0,     0,     0,     0,
  k_zfac_uv                          = 10,   10,    10,     10,
if_zfac_t                           = 0,     0,     0,     0,
  k_zfac_t                           = 10,   10,    10,     10,
if_zfac_q                           = 0,     0,     0,     0,
  k_zfac_q                           = 10,   10,    10,     10
guv                                 = 0.0003,     0.0003,     0.0003,   0.0003,
gt                                  = 0.0003,     0.0003,     0.0003,   0.0003,
gq                                  = 0.005,      0.005,      0.005,    0.005,
if_ramping                          = 1,
dtramp_min                          = 60.0,
io_form_gfdda                       = 2,
/
&dynamics
w_damping                           = 0,
diff_opt                            = 1,
km_opt                              = 4,
diff_6th_opt                        = 0,      0,      0,     0,
diff_6th_factor                     = 0.12,   0.12,   0.12,  0.12,
base_temp                           = 290.
damp_opt                            = 0,
zdamp                               = 5000.,  5000.,  5000., 5000.,
dampcoef                            = 0.2,    0.2,    0.2,   0.2,
khdif                               = 0,      0,      0,     0,
kvdif                               = 0,      0,      0,     0,
non_hydrostatic                     = .true., .true., .true., .true.,
h_mom_adv_order                     = 5,  5,  5,  5,
v_mom_adv_order                     = 3,  3,  3,  3,
h_sca_adv_order                     = 5,  5,  5,  5,
v_sca_adv_order                     = 3,  3,  3,  3,
moist_adv_opt                       = 1,  1,  1,  1,
scalar_adv_opt                      = 1,  1,  1,  1,
chem_adv_opt                        = 2,  2,  2,  2,
/
&bdy_control
spec_bdy_width                      = 5,
spec_zone                           = 1,
relax_zone                          = 4,
spec_exp                            = 0.33,
specified                           = .true.,.false.,.false.,.false.,
nested                              = .false.,.true.,.true.,.true.,
/

&scm
/
&noah_mp
/
&fire
/
&dfi_control
/
&tc
/
&diags
/
&logging
/

&grib2
/

&chem
kemit                               = 7,
chem_opt                            = 9,  9,  9,  9,
chem_in_opt                         = 1,  1,  1,  1,
bio_emiss_opt                       = 3,  3,  3,  3,
photdt                              = 10, 10, 10, 10,
chemdt                              = 1.5,  1.5,  1,  1,
io_style_emissions                  = 2,
emiss_inpt_opt                      = 101, 101, 101, 101,
emiss_opt                           = 4,   4,   4,   4,
phot_opt                            = 3,   3,   3,  3,
gas_drydep_opt                      = 1,   1,   1,  1,
aer_drydep_opt                      = 1,   1,   1,  1,
ne_area                             = 200,
biomass_burn_opt                    = 2,   2,   2,  2,
plumerisefire_frq                   = 60,  60,  60, 60,
dust_opt                            = 4,
dust_schme                          = 3,
dustwd_onoff                        = 1,
dmsemis_opt                         = 0,
seas_opt                            = 0,
opt_pars_out                        = 1,
gas_bc_opt                          = 1,   1,   1,  1,
gas_ic_opt                          = 1,   1,   1,  1,
aer_bc_opt                          = 1,   1,   1,  1,
aer_ic_opt                          = 1,   1,   1,  1,
gaschem_onoff                       = 1,   1,   1,  1,
aerchem_onoff                       = 1,   1,   1,  1,
wetscav_onoff                       = 1,   1,   1,  1,
cldchem_onoff                       = 1,   1,   1,  1,
vertmix_onoff                       = 1,   1,   1,  1,
chem_conv_tr                        = 0,   0,   0,  0,
aer_ra_feedback                     = 1,   1,   1,  1,
have_bcs_chem                       = .true., .true., .true.,.true.,
scale_fire_emiss                    = .false.,.false.,.false.,.false.,
! mozart_ph_diag                      = 1
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/