WRF-CHEM沙尘预报

汤圆YYE

求助大家！跑wrf-chem，加了人为源和沙尘模式，chem_opt=10,dust_opt=3,结果很不好，西北沙尘成片高值，想请教下大家可能是什么导致的这个结果？
&time_control
run_days                            = 0,
run_hours                           = 24,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2024, 2023, 2000,
start_month                         = 02, 07,   01,
start_day                           = 23, 01,   24,
start_hour                          = 15,   00,   12,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2024, 2023, 2000,
end_month                           = 02, 08,   01,
end_day                             = 24, 01,   25,
end_hour                            = 15,   00,   12,
end_minute                          = 00,   00,   00,
end_second                          = 00,   00,   00,
interval_seconds                    = 10800
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 60,  60,   60,
frames_per_outfile                  = 1000, 1000, 1000,
write_hist_at_0h_rst                = .true.,
restart                             = .false.,
restart_interval                    = 14400,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
io_form_auxinput5                   = 2
auxinput5_interval                  = 60
debug_level                         = 0
history_outname = 'wrfout_d<domain>'


&domains
time_step                           = 60,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 630,    86,   94,
e_sn                                = 400,    111,    91,
e_vert                              = 51,    29,    30,
p_top_requested                     = 1000,
num_metgrid_levels                  = 34,
num_metgrid_soil_levels             = 4,
dx                                  = 9000, 4000,  3333.33,
dy                                  = 9000, 4000,  3333.33,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     28,    30,
j_parent_start                      = 1,     30,    30,
parent_grid_ratio                   = 1,     5,     3,
parent_time_step_ratio              = 1,     5,     3,
feedback                            = 0,
smooth_option                       = 0,
rh2qv_wrt_liquid                    = .true.,
eta_levels = 1.00000,0.99381,0.98643,0.97786,0.96815,
    0.95731,0.94538,0.93122,0.91490,0.89653,0.87621,0.85405,0.82911,
    0.80160,0.77175,0.73981,0.70509,0.66798,0.62889,0.58823,0.54957,
    0.51281,0.47788,0.44471,0.41323,0.38336,0.35503,0.32819,0.30276,
    0.27869,0.25592,0.23439,0.21405,0.19484,0.17672,0.15963,0.14352,
    0.12836,0.11410,0.10070,0.08811,0.07630,0.06523,0.05487,0.04517,
    0.03611,0.02765,0.01977,0.01243,0.00560,0.00000,
/


&physics
mp_physics                          = 10,
num_land_cat                        = 20,
ra_lw_physics                       = 4,
ra_sw_physics                       = 4,
!radt                                = 10,
radt                                = 0,
!sf_sfclay_physics                   = 2,
sf_sfclay_physics                   = 5,
sf_surface_physics                  = 2,
bl_pbl_physics                      = 6,
bldt                                = 0,
cu_physics                          = 5,
cudt                                = 5,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,
topo_wind                           = 0,
progn                               = 1,
/

&fdda
/

&dynamics
w_damping                           = 1,
diff_opt                            = 1,   1,
km_opt                              = 4,   4,
use_theta_m                         = 1,
use_q_diabatic                      = 1,
diff_6th_opt                        = 1,   0,
diff_6th_factor                     = 0.2, 0.2,
epssm                               = 0.5, 0.3,
base_temp                           = 290.,
damp_opt                            = 1,
zdamp                               = 7000, 5000,
dampcoef                            = 0.01,
khdif                               = 0,
kvdif                               = 0,
non_hydrostatic                     = .true.,
moist_adv_opt                       = 1,
scalar_adv_opt                      = 1,
gwd_opt                             = 0,
/

&bdy_control
spec_bdy_width                      = 7,
spec_zone                           = 1,
relax_zone                          = 6,
specified                           = .true., .false.,.false.,
nested                              = .false., .true., .true.,
/

&grib2
/