请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3431|回复: 0

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

[复制链接]

新浪微博达人勋

发表于 2019-8-26 17:51:30 | 显示全部楼层 |阅读模式
GrADS
系统平台:
问题截图: -
问题概况: wrfchem模拟,O3和NO2没有结果
我看过提问的智慧: 看过
自己思考时长(天): 20

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

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

x
大家好! 我学习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的截图:
SO2.png
NO.png
O3的结果是这样的:
O3.png

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

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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