【求助】WRF-chem模拟出的primary PM2.5和BC都几乎为零

葡萄不是哭包

如题，我用wrfchem模拟颗粒物，但最终输出的初次PM2.5(P25)和黑碳的浓度为十的-13次方其他颗粒物还正常（硫酸盐，硝酸盐，铵盐等）

我的namelist设置为：
&time_control
!run_days                            = 18,
!run_hours                           = 06,
!run_minutes                         = 0,
!run_seconds                         = 0,
start_year                          = 2016, 2016, 2008,
start_month                         = 09,   09,   02,
start_day                           = 25,   25,   25,
start_hour                          = 18,   18,   00,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2016, 2016, 2008,
end_month                           = 11,   11,   05,
end_day                             = 01,   01,   31,
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                    = 180,  60,   60,
!frames_per_outfile                  = 80, 240, 8,
frames_per_outfile                   =10000,10000
restart                             = .false.,
restart_interval                    = 5000,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
io_form_auxinput5                   = 2
io_form_auxinput6                   = 0
io_form_auxinput8                   = 0
debug_level                         = 10
auxinput8_inname                    ="wrfchemi_gocart_bg_d<domain>"
!auxinput7_interval_m                = 21600, 21600, 21600,
auxinput8_interval_m                = 21600, 21600, 21600,
/

&domains
time_step                           = 90,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 65,  85,    316,
e_sn                                = 57, 58,    316,
e_vert                              =  30,    30,     30,
p_top_requested                     = 5000,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
dx                                  = 27000,9000,  3333.33,
dy                                  = 27000,9000,  3333.33,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     19,    32,
j_parent_start                      = 1,     24,    75,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     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
mp_physics                          = 2,     2,     3,
progn                               = 1,     1,
ra_lw_physics                       = 1,     1,     1,
ra_sw_physics                       = 2,     2,     1,
radt                                = 30,    30,    30,
sf_sfclay_physics                   = 1,     1,     1,
sf_surface_physics                  = 2,     2,     2,
sf_urban_physics                    = 2, 2,
bl_pbl_physics                      = 1,     1,     1,
bldt                                = 0,     0,     0,
cu_physics                          = 5,     5,     0,
cudt                                = 5,     5,     5,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
num_land_cat                        = 21,
cu_rad_feedback                     = .true.,.true.,
sf_urban_physics                    = 0,     0,     0,
/

&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.,
moist_adv_opt                       = 1,      1,      1,     
scalar_adv_opt                      = 1,      1,      1,     
/

&bdy_control
spec_bdy_width                      = 5,
spec_zone                           = 1,
relax_zone                          = 4,
specified                           = .true., .false.,.false.,
nested                              = .false., .true., .true.,
/

&grib2
/
&chem
chem_opt                            = 35,35,
kemit                               = 2,
io_style_emissions                  = 2,
bioemdt                             = 30,30,
photdt                              = 30,30,
chemdt                              = 1.5,1.5,
emiss_inpt_opt                      = 102,102,
emiss_opt                           = 8,8,
emiss_opt_vol                       = 0,0,
emiss_ash_hgt                       = 20000.,
chem_in_opt                         = 1,1,
phot_opt                            = 1,1,
gas_drydep_opt                      = 1,1,
aer_drydep_opt                      = 1,1,
bio_emiss_opt                       = 3,3,
dust_opt                            = 1,
depo_fact                           = 0.25,
dmsemis_opt                         = 1,
seas_opt                            = 1,
depo_fact                           = 0.25,0.25,
gas_bc_opt                          = 1,1,
gas_ic_opt                          = 1,1,
aer_ic_opt                          = 1,1,
aer_bc_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,
conv_tr_wetscav                     = 1,1,
conv_tr_aqchem                      = 1,1,
biomass_burn_opt                    = 1,1,
plumerisefire_frq                   = 30,30,
have_bcs_chem                       = .false.,.false.,
have_bcs_tracer                     = .false.,.false.,
aer_ra_feedback                     = 1,1,
aer_op_opt                          = 2,2,
opt_pars_out                        = 1,
diagnostic_chem                     = 0,
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/



请问这是怎么回事？？？

葡萄不是哭包

在顶一下帖子 请问大家有啥建议吗

wjmzbmr

可能是化学机制的选择有问题？

s646634007

auxinput_5 -> Anthropogenic emissions
auxinput_6 -> Biogenic emissions (BEIS, MEGAN)
auxinput_7 -> Surface biomass burning fields
auxinput_8 -> GOCART background fields
auxinput_12 -> Chemistry initial fields
auxinput_13 -> Volcanic ash emissions
auxinput_14 -> Aircraft emissions
auxinput_15 -> Green House Gas emissions
这是官方文档中的资料。冒昧的问一下，你的auxinput好像只有auxinput_8。不需要增加Anthropogenic emissions等其他排放数据吗？