爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
楼主: 悠然一梦123

WRF-CHEM沙尘模块该如何弄

[复制链接]
发表于 2018-4-20 12:18:58 | 显示全部楼层
lulu 发表于 2018-4-19 10:41
请问沙尘模式运行出来的数据,要计算哪些物理量啊,模式运行以后对结果不知所粗,希望能指点一二

你要用啥。。。可以ncdump -h wrfout* 看一下输出变量。
密码修改失败请联系微信:mofangbao
发表于 2018-4-20 15:54:45 | 显示全部楼层
astiny 发表于 2018-4-20 12:18
你要用啥。。。可以ncdump -h wrfout* 看一下输出变量。

我想提取沙尘含量和沙通量画图,想问一下前辈,我用的dust opt=1,输出值理有dust1到dust5五个量的值,我要求沙通量然后画图,应该是dust1到5相加吗?还是说dust1到5相加应该是沙尘含量呢?还有,需要乘以密度再相加才是μg/m3对吧?
密码修改失败请联系微信:mofangbao
发表于 2018-4-22 16:40:30 | 显示全部楼层
lulu 发表于 2018-4-20 15:54
我想提取沙尘含量和沙通量画图,想问一下前辈,我用的dust opt=1,输出值理有dust1到dust5五个量的值,我 ...

好像是五个相加。单位的话,我记得原来单位是ug/kg干空气。那样的话,是需要乘以空气密度,才能转化为微克每立方米的。
密码修改失败请联系微信:mofangbao
发表于 2018-4-24 09:17:35 | 显示全部楼层
astiny 发表于 2018-4-22 16:40
好像是五个相加。单位的话,我记得原来单位是ug/kg干空气。那样的话,是需要乘以空气密度,才能转化为微 ...

真的非常感谢!
密码修改失败请联系微信:mofangbao
发表于 2019-9-6 18:25:46 | 显示全部楼层
lulu 发表于 2018-4-20 15:54
我想提取沙尘含量和沙通量画图,想问一下前辈,我用的dust opt=1,输出值理有dust1到dust5五个量的值,我 ...

您好!打扰了。请问您是怎么提取沙尘含量和沙通量画图的?如果是ncl想请您告知修改的官网哪个脚本好吗?非常感谢~(初学所以很想及时看看跑出来的可视化结果~
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2019-10-5 16:36:13 | 显示全部楼层
本帖最后由 lynn1993 于 2019-10-5 16:37 编辑
xlilium 发表于 2019-9-6 18:25
您好!打扰了。请问您是怎么提取沙尘含量和沙通量画图的?如果是ncl想请您告知修改的官网哪个脚本好吗? ...

请问,您是如何设置namelist.input的参数呢?只做沙尘呢。我做出来的结果都是0呢,如果可以,麻烦您加我QQ一下呀。1165602270
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2019-10-10 11:24:01 | 显示全部楼层
lynn1993 发表于 2019-10-5 16:36
请问,您是如何设置namelist.input的参数呢?只做沙尘呢。我做出来的结果都是0呢,如果可以,麻烦您加我Q ...

不好意思才看到,以下是我的namelist.input参数,供参考。

&time_control
run_days                            = 2,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2010,
start_month                         = 07,
start_day                           = 14,
start_hour                          = 00,
start_minute                        = 00,
start_second                        = 00,
end_year                            = 2010,
end_month                           = 07,
end_day                             = 16,
end_hour                            = 00,
end_minute                          = 00,
end_second                          = 00,
interval_seconds                    = 21600,
input_from_file                     = .true.,
history_interval                    = 60,
frames_per_outfile                  = 72,
restart                             = .false.,
restart_interval                    = 0,
io_form_history                     = 2,
io_form_restart                     = 2,
io_form_input                       = 2,
io_form_boundary                    = 2,
auxinput6_inname                    = 'wrfbiochemi_d01',
auxinput7_inname                    = 'wrffirechemi_d<domain>',
auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>',
auxinput12_inname                   = 'wrf_chem_input',
auxinput5_interval_m                = 86400,
auxinput7_interval_m                = 86400,
auxinput8_interval_m                = 86400,
io_form_auxinput2                   = 2,
io_form_auxinput5                   = 0,
io_form_auxinput6                   = 0,
io_form_auxinput7                   = 0,
io_form_auxinput8                   = 0,
io_form_auxinput12                  = 0,
debug_level                         = 0,
auxinput1_inname                    = "met_em.d<domain>.<date>",
/
auxinput13_inname                   = 'wrfchemv_d<domain>',
auxinput13_interval_m               = 86400,
io_form_auxinput13                  = 0,

&dfi_control
/

&domains
time_step                           = 600,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
s_we                                = 1,
e_we                                = 41,
s_sn                                = 1,
e_sn                                = 41,
e_vert                              = 40,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 100000,
dy                                  = 100000,
grid_id                             = 1,
parent_id                           = 0,
i_parent_start                      = 1,
j_parent_start                      = 1,
parent_grid_ratio                   = 1,
parent_time_step_ratio              = 1,
p_top_requested                     = 5000,
feedback                            = 1,
smooth_option                       = 0,
p_top_requested                     = 5000,
zap_close_levels                    = 50,
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
num_land_cat                        = 21,
mp_physics                          = 4,
progn                               = 0,
ra_lw_physics                       = 1,
ra_sw_physics                       = 2,
radt                                = 30,
sf_sfclay_physics                   = 1,
sf_surface_physics                  = 2,
bl_pbl_physics                      = 1,
bldt                                = 0,
cu_physics                          = 5,
cu_diag                             = 1,
cudt                                = 0,
ishallow                            = 0,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,
mp_zero_out                         = 2,
mp_zero_out_thresh                  = 1.e-12,
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,
km_opt                              = 4,
diff_6th_opt                        = 0,
diff_6th_factor                     = 0.12,
base_temp                           = 290.
damp_opt                            = 0,
zdamp                               = 5000.,
dampcoef                            = 0.01,
khdif                               = 0,
kvdif                               = 0,
non_hydrostatic                     = .true.,
moist_adv_opt                       = 2,
scalar_adv_opt                      = 2,
chem_adv_opt                        = 2,
tke_adv_opt                         = 2,
time_step_sound                     = 4,
h_mom_adv_order                     = 5,
v_mom_adv_order                     = 3,
h_sca_adv_order                     = 5,
v_sca_adv_order                     = 3,
/

&bdy_control
spec_bdy_width                      = 5,
spec_zone                           = 1,
relax_zone                          = 4,
specified                           = .true.,
nested                              = .false.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/


&chem
kemit                               = 1,
chem_opt                            = 401,
bioemdt                             = 0,
photdt                              = 0,
chemdt                              = 10,
io_style_emissions                  = 0,
emiss_opt                           = 3,
emiss_opt_vol                       = 0,
emiss_ash_hgt                       = 20000.,
chem_in_opt                         = 0,
phot_opt                            = 0,
gas_drydep_opt                      = 0,
aer_drydep_opt                      = 1,
bio_emiss_opt                       = 0,
ne_area                             = 0,
dust_opt                            = 1,
dmsemis_opt                         = 0,
seas_opt                            = 0,
depo_fact                           = 0.25,
gas_bc_opt                          = 0,
gas_ic_opt                          = 0,
aer_bc_opt                          = 1,
aer_ic_opt                          = 1,
gaschem_onoff                       = 0,
aerchem_onoff                       = 1,
wetscav_onoff                       = 0,
cldchem_onoff                       = 0,
vertmix_onoff                       = 1,
chem_conv_tr                        = 0,
conv_tr_wetscav                     = 0,
conv_tr_aqchem                      = 0,
biomass_burn_opt                    = 0,
plumerisefire_frq                   = 30,
have_bcs_chem                       = .false.,
aer_ra_feedback                     = 0,
aer_op_opt                          = 0,
opt_pars_out                        = 0,
diagnostic_chem                     = 0,
/
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2019-10-30 22:17:01 | 显示全部楼层
谢谢谢谢,已经解决这个问题
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2019-11-27 21:26:02 | 显示全部楼层
lynn1993 发表于 2019-10-30 22:17
谢谢谢谢,已经解决这个问题

请问你是怎么解决的呢?我也遇到了同样的问题?对比上面的namelist没有错啊
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2020-3-29 10:46:40 | 显示全部楼层
xlilium 发表于 2019-10-10 11:24
不好意思才看到,以下是我的namelist.input参数,供参考。

&time_control

您好,能请教一下chem里的起沙机制该怎么设置吗?是用dust schem这个参量吗?最近想用shao方案,但不知道该在哪个变量里选。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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