求助：wrfchem运行的结果里PM的值偏高太多。。。。

绛紫

我现在运行的是wrf-chem的3层嵌套，用idv查看wrfout结果的时候，PM值偏高太多，检查过输入文件都是没有问题的，不知道是哪里出了问题，下边是我运行的namelist和结果图，想请问下有没有人可以帮忙解答的，谢谢~~
namelist如下：
&time_control
run_days                            = 0,
run_hours                           = RunHour,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = STYear,  STYear,   STYear,
start_month                         = STMonth, STMonth,   STMonth,
start_day                           = STDay,   STDay,      STDay,
start_hour                          = STHour,   STHour,    STHour,
start_minute                        = 00,  00,  00,
start_second                        = 00,  00,  00,
end_year                            = EndYear,  EndYear,   EndYear,
end_month                           = EndMonth, EndMonth,   EndMonth,
end_day                             = EndDay,  EndDay,     EndDay,
end_hour                            = EndHour, EndHour,    EndHour,
end_minute                          = 00, 00,  00,
end_second                          = 00, 00,  00,
interval_seconds                    = 21600,
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 60,   60,   60,
frames_per_outfile                  = 1,     1,    1,
restart                             = .false.,
restart_interval                    = 4000,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
io_form_auxinput5                   = 2
io_form_auxinput12                  = 2
auxinput12_inname                   = 'wrf_chem_input'
debug_level                         = 0
/
&domains
time_step                           = 180,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 3,
e_we                                = 186,    61,   64,
e_sn                                = 155,    88,   64,
e_vert                              = 30,    30,    30,
dx                                  = 27000, 9000, 3000,
dy                                  = 27000, 9000, 3000,
p_top_requested                     = 5000
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
grid_id                             = 1,     2,     3,
parent_id                           = 1,     1,     2,
i_parent_start                      = 1,     128,   19,
j_parent_start                      = 1,     76,    46,
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,     1,
naer                                = 1,
ra_lw_physics                       = 1,     1,     1,
ra_sw_physics                       = 2,     2,     2,
radt                                = 27,    27,    27,
sf_sfclay_physics                   = 2,     2,     2,
sf_surface_physics                  = 2,     2,     2,
bl_pbl_physics                      = 2,     2,     2,
bldt                                = 0,     0,     0,
cu_physics                          = 5,     5,     5,
cu_diag                             = 1,
cudt                                = 0,     0,     0,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,     0,     0,
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
cu_rad_feedback                     = .ture.,
/
&fdda
/
&dynamics
w_damping                           = 1,
diff_opt                            = 1,
km_opt                              = 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.01,    0.01,    0.01
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,
chem_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
kemit                               = 5,
chem_opt                            = 2,           2,         2,
bioemdt                             = 30,
photdt                              = 30,
chemdt                              = 1.5,
emiss_inpt_opt                      = 1,          1,          1,
io_style_emissions                  = 2,
emiss_opt                           = 3,          3,          3,
chem_in_opt                         = 0,          0,          0,
phot_opt                            = 2,          2,          2,
gas_drydep_opt                      = 1,          1,          1,
aer_drydep_opt                      = 1,          1,          1,
bio_emiss_opt                       = 1,          1,          1,
dust_opt                            = 0,
dmsemis_opt                         = 0,
seas_opt                            = 0,
gas_bc_opt                          = 1,          1,          1,
gas_ic_opt                          = 1,          1,          1,
aer_bc_opt                          = 1,          1,          1,
aer_ic_opt                          = 1,          1,          1,
gaschem_onoff                       = 1,          1,          1,
aerchem_onoff                       = 1,          1,          1,
wetscav_onoff                       = 0,          0,          0,
cldchem_onoff                       = 0,          0,          0,
vertmix_onoff                       = 1,          1,          1,
chem_conv_tr                        = 1,          1,          1,
biomass_burn_opt                    = 0,          0,          0,
plumerisefire_frq                   = 30,         30,         30,
have_bcs_chem                       = .false., .false., .false.,
aer_ra_feedback                     = 1,
diagnostic_chem                     = 1,
/
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
idv查看的结果如下图：

ytj584684235

楼主你好，我想问下里面的 naer=1 是代表1000/cm^3吗？

ytj584684235

diagnostic_chem=1可以多输出什么东西？

绛紫

ytj584684235 发表于 2015-4-2 12:40
楼主你好，我想问下里面的 naer=1 是代表1000/cm^3吗？

您好，我没看明白你说的1000/cm^3是说输出的污染物单位是吗？额。。。。

ytj584684235

这个参数不是设置气溶胶浓度的吗？指的是气溶胶浓度，你理解呢？

Alrescha

lz你的wrfchemi里面的pm值合理吗？
还有，楼主你的chem_in_opt是不是应该设为1呢？

绛紫

ytj584684235 发表于 2015-4-2 16:34
这个参数不是设置气溶胶浓度的吗？指的是气溶胶浓度，你理解呢？

是啊，它的单位是ug/m^3啊

ytj584684235

请问你在哪里看到的我怎么理解的是number aerosol的简称

ytj584684235

看到一些文章里说到的   个/cm^3   请问楼主知道在这里面如何设置吗？

绛紫

Alrescha 发表于 2015-4-2 17:46
lz你的wrfchemi里面的pm值合理吗？
还有，楼主你的chem_in_opt是不是应该设为1呢？

wrfchemi里边的PM是合理的啊，分布也是对的，所以不知道是哪里的问题，然后那个chem_in_opt第一天不是设为0吗，第二天的设为1，我贴的这个namelist是我第一天的，所以是0。