使用ANTHRO（anthro_emis）工具处理EDGAR排放清单，WRF-CHEM模拟值极低，几乎无浓度

15104562864

各位同学，有个问题请教大家。我使用WRF-CHEM（v4.0）默认清单处理工具ANTHRO（anthro_emis）处理EDGAR排放清单（使用EDGAR2020套壳EDGAR_HTAP 2010），输出小时值清单，比如wrfchemi_d01_2020-02-08_20:00:00。24核并行运算不报错，模拟结果中PM2.5、NH3等浓度比实际值低3~4个数量级，几乎没有浓度。排查以后我觉得问题在两个方面:

第一，anthro_emis处理过程中，因受wrf-chem、edgar等版本变化影响，清单物种的单位转换是否存在数量级差异。

第二，namelist.input设置错误，没有读取清单，清单仅仅使用极低的初始浓度，越积分越接近于0。阅读了一些气象家园帖子，其中有人说auxinput5_interval_m、chem_opt、io_form_auxinput5、emiss_inpt_opt等的设置都会影响模拟，比如我把auxinput5_interval_m从60改成100000，模拟值就变为了比较正常的浓度；再比如之前报错不识别排放清单第二天的第一个小时，我按照一些帖子把force_use_old_data=.true.加上去，可以模拟不报错，但浓度极低。

有没有做过类似工作的同学，求赐教。

图片是我的anthro_emis的分配文件

附namelist.input，这个namelist.input已经被我改的面目全非了，请大家看看有哪些问题。

若同学可提供成功的namelist.input，有偿。

&time_control

run_days                            = 10,

run_hours                           = 0,

run_minutes                         = 0,

run_seconds                         = 0,

start_year                          = 2020, 2020, 2020,

start_month                         = 01,    01,   01,

start_day                           = 01,    01,   01,

start_hour                          = 00,   00,   00,

start_minute                        = 00,   00,   00,

start_second                        = 00,   00,   00,

end_year                            = 2020, 2020, 2020,

end_month                           = 01,    01,   01,

end_day                             = 10,    10,   10,

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,  180,   60,

frames_per_outfile                  = 50000, 50000, 50000,

restart                             = .false.,

write_hist_at_0h_rst                = .true.

restart_interval                    = 1440,

io_form_history                     = 2,

io_form_restart                     = 2,

io_form_input                       = 2,

io_form_boundary                    = 2,

io_form_auxinput2                   = 2,

auxinput1_inname                    = "met_em.d<domain>.<date>",

auxinput5_inname                    = "wrfchemi_d<domain>_<date>",

auxinput6_inname                    = 'wrfbiochemi_d<domain>',

!auxinput7_inname                    = 'wrffirechemi_d<domain>',

!auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>',

auxinput12_inname                   = 'wrf_chem_input',

auxinput5_interval_m                = 100000, 100000, 100000,

! auxinput7_interval                  = 21600, 21600, 21600,

! auxinput8_interval                  = 21600, 21600, 21600,

io_form_auxinput5                   = 2,

io_form_auxinput6                   = 2,

io_form_auxinput7                   = 0,

io_form_auxinput8                   = 0,

io_form_auxinput12                  = 0,

io_form_auxinput13                  = 0,

io_form_auxinput14                  = 0,

io_form_auxinput15                  = 0,

frames_per_auxinput5                = 1,

debug_level                         = 0,

force_use_old_data=.true.

/

&domains

time_step                           = 180,

time_step_fract_num                 = 0,

time_step_fract_den                 = 1,

max_dom                             = 2,

e_we                     = 118,      211,       73,

e_sn                     = 85,      169,       73,

e_vert                              = 35,    35,    35,

p_top_requested                     = 5000,

num_metgrid_levels                  = 34,

num_metgrid_soil_levels             = 4,

dx                       = 81000,     27000,   9000,

dy                       = 81000,     27000,   9000,

grid_id                 = 1,     2,     3,

parent_id                = 1,        1,        2,

i_parent_start           = 1,       23,       26,

j_parent_start           = 1,       18,       26,

parent_grid_ratio       = 1,     3,     3,

parent_time_step_ratio  = 1,     3,     3,

feedback                            = 0,

smooth_option                       = 0

/

&physics

mp_physics                          = 10,     2,     2,

progn                               = 0,     0,     0,

naer                                = 1e9

ra_lw_physics                       = 1,     1,     1,

ra_sw_physics                       = 2,     2,     2,

radt                                = 30,    30,    30,

sf_sfclay_physics                   = 1,     1,     1,

sf_surface_physics                  = 2,     2,     2,

bl_pbl_physics                      = 1,     1,     1,

bldt                                = 0,     0,     0,

cu_physics                          = 5,     5,     5,

cudt                                = 5,     5,     5,

isfflx                              = 1,

ifsnow                              = 0,

icloud                              = 1,

surface_input_source                = 3,

num_soil_layers                     = 4,

num_land_cat                        = 21,

sf_urban_physics                    = 0,     0,     0,

maxiens                             = 1,

maxens                              = 3,

maxens2                             = 3,

maxens3                             = 16,

ensdim                              = 144,

cu_rad_feedback                     = .true.,

cu_diag                             = 1

/

&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,     

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                               = 1,

chem_opt                            = 2,         2,       2,

photdt                              = 30,       30,      30,

chemdt                              = 2.,  0.66666,   0.22222,

!frames_per_emissfile               = 1,         1,       1,

io_style_emissions                  = 2,

emiss_inpt_opt                      = 1,        1,        1,

emiss_opt                           = 3,        3,        3,

chem_in_opt                         = 0,        0,        0,

phot_opt                            = 1,        1,        1,

gas_drydep_opt                      = 0,        0,        0,

aer_drydep_opt                      = 0,        0,        0,

bio_emiss_opt                       = 0,        0,        0,

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                        = 0,        0,        0,

biomass_burn_opt                    = 0,        0,        0,

plumerisefire_frq                   = 30,       30,      30,

aer_ra_feedback                     = 0,        0,        0,

have_bcs_chem                       = .false., .false., .false.,

/

&namelist_quilt

nio_tasks_per_group = 0,

nio_groups = 1,

/