求助!!!wrfchem 臭氧模拟不出来

提妮

大家好! 我学习wrfchem有一小段时间了,想用wrfchem模拟海南地区的臭氧,但一直没有结果.PM2.5,PM10,NO都有结果.
我做了很多尝试,比如是不是我模拟的嵌套区域太小了,于是我现在改成中国地区只一层嵌套看看结果,但也模拟不出来;是不是我的排放源和chem_opt不匹配?我的排放源是用wrfchem官网的Mozart-mosaic做的,做完也是显示成功的,先开始我编译wrfchem设置的kpp=0,所选的chem_opt=201,但手册上说chem_opt=201是Mozart_mosaic并且用kpp库,于是我又尝试加kpp库结果仍是如此.非常苦恼,来家园这里求助,希望各位大佬帮帮我!!
PS:我在各种过程中都没有显示出错
这里附上我的namelist.input文件:

&time_control

run_days                            = 0,

run_hours                           = 42,

run_minutes                         = 0,

run_seconds                         = 0,

start_year                          = 2017,  2017,   2017,

start_month                         = 07,         07,         07,

start_day                           = 01,           01,           01,

start_hour                          = 00,          00,         00,

start_minute                        = 00,           00,           00,

start_second                        = 00,           00,           00,

end_year                            = 2017,  2018,   2018,

end_month                           = 07,         07,        07,

end_day                             = 02,           02,           02,

end_hour                            = 18,          00,         00,

end_minute                          = 00,           00,           00,

end_second                          = 00,           00,           00,

interval_seconds                    = 3600

input_from_file                     = .true.,.true.,.true.,

history_interval                    =  180,          60,   60,

frames_per_outfile                  =  240, 240, 240,

restart                             = .false.,

restart_interval                    = 10000,

!2=netcdf,11=pnetcdf

io_form_history                     = 2

io_form_restart                     = 2

io_form_input                       = 2

io_form_boundary                    = 2

io_form_auxinput2                   = 2,

!Antropogenic emission

io_form_auxinput5                   = 2,

auxinput5_interval                  = 3600,

!Biogenic emission

io_form_auxinput6                   = 0,

auxinput6_inname                    = 'wrfbiochemi_d<domain>',

!Biomass buring

io_form_auxinput7                   = 0,

auxinput7_inname                    = 'wrffirechemi_d<domain>_<date>',

frames_per_auxinput7                             = 1, 1, 1,

!GOCART BACKGROUND

io_form_auxinput8                   = 0,

auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>',

!Use the previous WRF-OUT

io_form_auxinput12                  = 0,

auxinput12_inname                   = 'wrf_chem_input',

!Volcanic ash emission

io_form_auxinput13                  = 0,

auxinput13_inname                   = 'wrfchemv_d<domain>',

!Aircraft emission

io_form_auxinput14                  = 0,

!GHG input

io_form_auxinput15                  = 0,

debug_level                         = 10,

/

&dfi_control

/

&domains

time_step                           = 120,

max_dom                              = 1,

e_we                                = 230,   109,    196,

e_sn                                = 220,    109,    196,

e_vert                              = 30,    30,    30,

p_top_requested                     = 5000,

num_metgrid_levels                  = 32,

num_metgrid_soil_levels             = 4,

dx                                  = 27000,  9000,  3000,

dy                                  = 27000,  9000,  3000,

grid_id                             = 1,     2,     3,

parent_id                           = 1,     1,     2,

i_parent_start                      = 1,     45,   37,

j_parent_start                      = 1,     30,   11,

parent_grid_ratio                   = 1,     3,     3,

parent_time_step_ratio              = 1,     3,     3,

feedback                            = 1,

smooth_option                       = 0

interp_type                         = 1

t_extrap_type                       = 2

force_sfc_in_vinterp                = 0

use_levels_below_ground             = .true.

use_surface                         = .true.

lagrange_order                      = 1

/

sfcp_to_sfcp                        = .true.

&physics

mp_physics                          = 10,     10,     10,

progn                               = 0,

ra_lw_physics                       = 4,     4,     4,

ra_sw_physics                       = 4,     4,     4,

radt                                = 30,    10,    3.33334,

sf_sfclay_physics                   = 1,     1,     1,

sf_surface_physics                  = 1,     1,     1,

bl_pbl_physics                      = 5,     5,     5,

bldt                                = 0,     0,     0,

cu_physics                          = 5,     5,     0,

cu_diag                             = 0,

cudt                                = 0,     0,     0,

ishallow                            = 0,

isfflx                              = 1,

ifsnow                              = 1,

icloud                              = 1,

surface_input_source                = 1,

num_soil_layers                     = 4,

sf_urban_physics                    = 0,     1,     1,

num_land_cat                        = 21,

mp_zero_out                         = 2,

mp_zero_out_thresh                  = 1.e-12

!GRELL Only

maxiens                             = 1,

maxens                              = 3,

maxens2                             = 3,

maxens3                             = 16,

ensdim                              = 144,

cu_rad_feedback                     = .true.,

/

&fdda

/

&dynamics

rk_ord                              = 3,

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.01,   0.01,   0.01

khdif                               = 0,      0,      0,

kvdif                               = 0,      0,      0,

non_hydrostatic                     = .true., .true., .true.,

moist_adv_opt                       = 2,      0,      0,

scalar_adv_opt                      = 2,      0,      0,

chem_adv_opt                        = 2,      0,      0,

tke_adv_opt                         = 2,      0,      0,

time_step_sound                     = 4,      4,      4,

h_mom_adv_order                     = 5,      5,      5,

v_mom_adv_order                     = 3,      3,      3,

h_sca_adv_order                     = 5,      5,      5,

v_sca_adv_order                     = 3,      3,      3,

/

&bdy_control

spec_bdy_width                      = 5,

spec_zone                           = 1,

relax_zone                          = 4,

specified                           = .true., .false.,.false.,

nested                              = .false.,.true.,.true.,

/

&grib2

/

&namelist_quilt

nio_tasks_per_group = 0,

nio_groups = 1,

/

&chem

kemit                               = 1,

chem_opt                            = 201,     201,    201,

chem_in_opt                         = 0,             0,        0,

chemdt                              = 3,        1,        0.333334,

phot_opt                            = 3,            3,        3,

photdt                              = 30,        30,        30,

aer_ra_feedback                     = 0,        0,        0,

io_style_emissions                  = 1,

emiss_opt                           = 10,             10,        10,

emiss_inpt_opt                             = 102,        102,        102,

bio_emiss_opt                       = 0,        0,        0,

bioemdt                             = 30,        6,        6,

ne_area                             = 1000,

biomass_burn_opt                    = 0,             0,        0,

plumerisefire_frq                   = 30,             30,        30,

scale_fire_emiss                     = .true.,        .true.,        .true.,

emiss_opt_vol                       = 0,              0,        0,

emiss_ash_hgt                       = 20000.,  

aircraft_emiss_opt                     = 0,        0,        0,

dust_opt                            = 3,

seas_opt                            = 3,

dmsemis_opt                         = 0,

gas_drydep_opt                      = 1,             1,        1,

aer_drydep_opt                      = 1,             1,        1,

depo_fact                           = 0.25

!bc,ic 1=default bc/ic

have_bcs_upper                      = .false., .false., .false.,

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

cldchem_onoff                       = 0,       0,        0,

vertmix_onoff                       = 1,       0,        0,

chem_conv_tr                        = 1,              0,        0,

conv_tr_wetscav                     = 0,              0,        0,

conv_tr_aqchem                      = 0,       0,        0,

fixed_ubc_inname                     = "ubvals_b40.20th.track1_1996-2005.nc"

aer_op_opt                          = 5,

opt_pars_out                        = 0,

diagnostic_chem                     = 0,

/

贴上ncview的截图:

O3的结果是这样的:

是不是我的namelist里头哪个数出了问题,恳请大家给我指点,谢谢!