爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6029|回复: 1

请问WRF Chem 为什么没有进行化学反应呢

[复制链接]

新浪微博达人勋

发表于 2022-7-17 18:23:23 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
请问各位老师,我用mozbc驱动的wrf chem为什么输出的化学反应得到的臭氧浓度始终为0呢,也就是模拟过程中没有进行化学反应,namelist如下
&time_control
run_days                            = 3,
run_hours                           = 12,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2008, 2008, 2008,
start_month                         = 02,   02,   02,
start_day                           = 24,   24,   24,
start_hour                          = 00    00,   00,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2008, 2008, 2008,
end_month                           = 02,   02,   02,
end_day                             = 27,   27,   27,
end_hour                            = 12,   12,   12,
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,  1000,
restart                             = .false.,
restart_interval                    = 0,
io_form_history                     = 2,
io_form_restart                     = 2,
io_form_input                       = 2,
io_form_boundary                    = 2,
auxinput5_interval_m                = 1440, 1440, 1440,
auxinput7_interval_m                = 1440, 1440, 1440,
auxinput8_interval_m                = 1440, 1440, 1440,
auxinput13_interval_m               = 1440, 1440, 1440,
io_form_auxinput2                   = 2,
io_form_auxinput5                   = 0,
io_form_auxinput6                   = 0,
io_form_auxinput7                   = 0,
io_form_auxinput8                   = 0,
io_form_auxinput12                  = 0,
io_form_auxinput13                  = 0,
frames_per_auxinput5                = 24, 24,24,
frames_per_auxinput6                = 24, 24,24,
frames_per_auxinput7                = 24, 24,24,
frames_per_auxinput8                = 24, 24,24,
frames_per_auxinput12               = 24, 24,24,
frames_per_auxinput13               = 24, 24,24,
debug_level                         = 0,
auxinput1_inname                    = "met_em.d<domain>.<date>",
/
&dfi_control
/

&domains
time_step                           = 100,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 2,
s_we                                = 1,     1,     1,
e_we                                = 109,   176,   91,
s_sn                                = 1,     1,     1,
e_sn                                = 104 ,  161,
e_vert                              = 45,    45,    48,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 30000, 6000, 4000,
dy                                  = 30000, 6000, 4000,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     31,   44,
j_parent_start                      = 1,     36,   33,
parent_grid_ratio                   = 1,     5,     3,
parent_time_step_ratio              = 1,     5,     3,
p_top_requested                     = 5000,
feedback                            = 1,
smooth_option                       = 1,
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                          = 2,     2,     2,
progn                               = 1,     1,     1,
ra_lw_physics                       = 1,     1,     1,
ra_sw_physics                       = 2,     2,     2,
radt                                = 30,    30,       4,
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,     0,
cu_diag                             = 1,     1,     0,
cudt                                = 5,     5,     5,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
num_land_cat                        = 21,
sf_urban_physics                    = 0,     0,     0,
mp_zero_out                         = 2,
mp_zero_out_thresh                  = 1.e-12
cu_rad_feedback                     = .true.,
/


&fdda
/

&dynamics
rk_ord                              = 3,
w_damping                           = 0,
diff_opt                            = 1,1,1,
km_opt                              = 4,4,4,
diff_6th_opt                        = 0,0,0,
diff_6th_factor                     = 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,      2,      2,
scalar_adv_opt                      = 2,      2,      2,
chem_adv_opt                        = 2,      2,      2,
tke_adv_opt                         = 2,      2,      2,
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                            = 1,   1,    301,
bioemdt                             = 30,  30,
photdt                              = 30,  30,
chemdt                              = 10,    2,      1,
io_style_emissions                  = 1,
emiss_opt                           = 7,     7,      0,
emiss_inpt_opt                      = 111,  111,      0,
emiss_opt_vol                       = 0,     0,      0,
emiss_ash_hgt                       = 20000.,
chem_in_opt                         = 1,      1,     0,
phot_opt                            = 2,      2,     2,
gas_drydep_opt                      = 0,      0,     0,
aer_drydep_opt                      = 1,      1,     1,
bio_emiss_opt                       = 2,      2,     0,
dust_opt                            = 1,
dust_schme                          = 1,
dmsemis_opt                         = 1,
seas_opt                            = 1,
depo_fact                           = 0.25,
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                       = 0,     0,      0,
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,      0,
conv_tr_wetscav                     = 1,     1,      1,
conv_tr_aqchem                      = 1,     1,      1,
biomass_burn_opt                    = 3,     3,      0,
plumerisefire_frq                   = 30,    30,    30,
have_bcs_chem                       = .true., .false., .false.,
aer_ra_feedback                     = 0,
aer_op_opt                          = 1,
opt_pars_out                        = 0,
diagnostic_chem                     = 1,     1,
chemdiag                            = 1,     1,
/

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2022-7-18 16:27:12 | 显示全部楼层
有没有大佬帮下忙
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表