无法从wrf-chem模拟结果中提取AOD数值

Stella@88

最近在学习wrf-chem,在跑完wrf-chem时遇到一些问题。我是用wrf-chemv3.5模拟pm2.5, O3, 和AOD数值。我是用RACM Chemistry and MADE/VBS aerosols (chem_opt=108)。但是模拟的结果中只有EXTCOF106, EXTCOF55, EXTCOF3数值都是Nan。
我的namelist.inp如下:
&time_control
run_days                            = 36,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2011,
start_month                         = 04,
start_day                           = 28,
start_hour                          = 00,
start_minute                        = 00,
start_second                        = 00,
end_year                            = 2011,
end_month                           = 06,
end_day                             = 02,
end_hour                            = 00,
end_minute                          = 00,
end_second                          = 00,
interval_seconds                    = 21600,
input_from_file                     = T,
history_interval                    = 60,
frames_per_outfile                  = 24,
restart                             = F,
cycling                             = F,
restart_interval                    = 1440,
io_form_history                     = 2,
io_form_restart                     = 2,
io_form_input                       = 2,
io_form_boundary                    = 2,
debug_level                         = 0,
diag_print                          = 0,
all_ic_times                        = .false.,
adjust_output_times                 = .false.,
override_restart_timers             = .false.,
write_hist_at_0h_rst                = .false.,
auxinput1_inname                    = "met_em.d<domain>.<date>"
auxinput4_inname                    = "wrflowinp_d<domain>"
auxinput4_interval                  = 360
io_form_auxinput4                   = 2
auxinput7_inname                    = "wrffirechemi_d<domain>_<date>"
io_form_auxinput7                   = 2
frames_per_auxinput7                = 1
auxinput7_interval_m                = 60
auxinput6_inname                    = 'wrfbiochemi_d01'
io_form_auxinput6                   = 2
! auxinput6_interval_h                = 24
auxinput5_interval_m                = 60
io_form_auxinput5                   = 2
/
&domains
time_step                           = 120,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
time_step_dfi                       = 60,
max_dom                             = 1,
s_we                                = 1,
e_we                                = 208,
s_sn                                = 1,  
e_sn                                = 190,
s_vert                              = 1,  
e_vert                              = 31,  
p_top_requested                     = 5000,
eta_levels = 1.000, 0.993, 0.980, 0.966, 0.950,
              0.933, 0.913, 0.892, 0.869, 0.844,
              0.816, 0.786, 0.753, 0.718, 0.680,
              0.639, 0.596, 0.550, 0.501, 0.451,
              0.398, 0.345, 0.290, 0.236, 0.188,
              0.145, 0.108, 0.075, 0.046, 0.021,
              0.000
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 18000, 10000,  3333.33,
dy                                  = 18000, 10000,  3333.33,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
interp_type                         = 2
extrap_type                         = 2
t_extrap_type                       = 2
use_levels_below_ground             = .true.
use_surface                         = .true.
lagrange_order                      = 1
zap_close_levels                    = 500
lowest_lev_from_sfc                 = .false.
force_sfc_in_vinterp                = 1
sfcp_to_sfcp                        = .false.
smooth_cg_topo                      = .false.
use_tavg_for_tsk                    = .false.
aggregate_lu                        = .false.
/
&dfi_control
dfi_opt                             = 0
dfi_nfilter                         = 7
/
&physics
mp_physics                          = 2,     3,     3,
ra_lw_physics                       = 4,     1,     1,
ra_sw_physics                       = 4,     1,     1,
radt                                = 30,    30,    30,
sf_sfclay_physics                   = 5,     1,     1,
sf_surface_physics                  = 2,     2,     2,
bl_pbl_physics                      = 5,     1,     1,
bldt                                = 0,     0,     0,
cu_physics                          = 5,     1,     0,
cu_diag                             = 1,     1,     0,
cudt                                = 0,     5,     5,
isfflx                              = 1,
ifsnow                              = 0,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,
num_land_cat                        = 20,
mp_zero_out                         = 2,
cu_rad_feedback                     = .true.
cugd_avedx                          = 1,
grav_settling                       = 1,
maxiens                             = 1,
sst_update                          = 1,
progn                               = 1,
ishallow                            = 0,
/
&fdda
grid_fdda              =      2,      1,      0,
gfdda_inname           = "wrffdda_d<domain>",
gfdda_end_h            = 99999, 99999, 99999, 99999, 99999
gfdda_interval_m       = 360,  360,  360,
io_form_gfdda          =      2,
fgdt                   =      0,      0,      0,
if_no_pbl_nudging_uv   =      0,      0,      0,
if_no_pbl_nudging_t    =      0,      0,      0,
if_no_pbl_nudging_q    =      0,      0,      0,
if_no_pbl_nudging_ph   =      0,      0,      0,
if_zfac_uv             =      0,      0,      0,
  k_zfac_uv             =     10,     10,     10,
if_zfac_t              =      0,      0,      0,
  k_zfac_t              =     10,     10,     10,
if_zfac_q              =      0,      0,      0,
  k_zfac_q              =     10,     10,     10,
if_zfac_ph             =      0,      0,      0,
  k_zfac_ph             =     10,     10,     10,
guv                    = 2.8e-4, 0.9e-4, 1.4e-4,
gt                     = 1.4e-4, 0.9e-4, 1.4e-4,
gq                     = 1.4e-4, 1.4e-4, 1.4e-4,
gph                    = 1.4e-4, 1.4e-4, 1.4e-4,
xwavenum               = 8
ywavenum               = 6
if_ramping             =      1,
dtramp_min             =  360.0,
obs_nudge_opt          =    0,    0,    0,
max_obs                =  90000,
fdda_start             =      0,     0,     0,     0,     0
fdda_end               =  99999, 99999, 99999, 99999, 99999
obs_nudge_wind         =      1,      1,      1,
obs_coef_wind          = 1.4E-4, 1.4E-4, 1.4E-4,
obs_nudge_temp         =      1,      1,      1,
obs_coef_temp          = 2.8E-4, 1.4E-4, 2.8E-4,
obs_nudge_mois         =      0,      0,     0,
obs_coef_mois          = 5.0E-4, 5.0E-4, 6.0E-4,
obs_rinxy              =    45.,    15.,    10.,
obs_rinsig             =   0.01,
obs_twindo             = 0.6667, 0.6667, 0.6667,
obs_npfi               =     40,
obs_ionf               =     40,     40,     40,
obs_idynin             =      0,
obs_dtramp             =    60.,
obs_ipf_errob          = .true.
obs_ipf_nudob          = .true.
obs_ipf_in4dob         = .true.
/
&dynamics
w_damping                           = 0,
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.2,    0.2
khdif                               = 0,      0,      0,
kvdif                               = 0,      0,      0,
non_hydrostatic                     = .true., .true., .true.,
moist_adv_opt                       = 2,      1,      1,     
scalar_adv_opt                      = 2,      1,      1,
chem_adv_opt                        = 2,   
/
&bdy_control
spec_bdy_width                      = 5,
spec_zone                           = 1,
relax_zone                          = 4,
specified                           = .true., .false.,.false.,
nested                              = .false., .true., .true.,
spec_exp                            = 0.33,
/
&tc
insert_bogus_storm                  = .false.
remove_storm                        = .false.
/
&grib2
/
&namelist_quilt
nio_tasks_per_group    = 0,
nio_groups             = 1,
/
&chem
kemit                  =  30,
chem_opt               =  108,    1,     1,     1,
photdt                 = 30,   30,    30,    30,
chemdt                 = 2.0, 2.5, 0.833, 0.093,   
io_style_emissions     =  2,
chem_in_opt            = 0 , 0 , 0 , 0 ,
have_bcs_chem          = .true., .true., .true., .true.,
phot_opt               =    2,    2,    1,     1,
gas_drydep_opt          =    1,    1,    1,     1,
aer_drydep_opt          =    1,    1,    1,     1,
emiss_opt              =    3,    3,    3,     3,
emiss_inpt_opt         =  1,    1,    1,     3,
bio_emiss_opt          =    3,
ne_area                = 103
bioemdt                = 30
dust_opt               =    3,
seas_opt               =    2,
depo_fact              = 0.25
biomass_burn_opt       =   2,    1,    1,    1,
plumerisefire_frq      = 30
gas_bc_opt             =    1,    1,    0,    0,
gas_ic_opt             =    1,    1,    0,    0,
aer_bc_opt             =    1,    1,    0,    0,
aer_ic_opt             =    1,    1,    0,    0,
gaschem_onoff          =    1,    1,    0,    0,
aerchem_onoff          =    1,    1,    0,    0,
wetscav_onoff          =    0,    1,    0,    0,
cldchem_onoff          =    0,    1,    0,    0,
vertmix_onoff          =    1,    1,    1,    1,
chem_conv_tr           =    1,    1,    1,    1,
aer_ra_feedback        =    1,    1,    1,    1,
aer_op_opt             =    1,    1,    1,    1,
opt_pars_out           =    1,
/
&diags
p_lev_diags            = 0,
num_press_levels       = 0,


请问是否是我的namelist设置的有没有哪里写的不对而导致的EXTCOF数值都是Nan？谢谢大家的帮忙

ll524302712

Turn on the aerosol optical properties (e.g., aer_op_opt, opt_pars_out) to obtain the extinction coefficient in the wrfout file and rerun the simulation.

ll524302712

Turn on the aerosol optical properties (e.g., aer_op_opt, opt_pars_out) to obtain the extinction coefficient in the wrfout file and rerun the simulation.