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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 13087|回复: 5

wrfchem输出如何得到污染物地表浓度或者柱浓度

[复制链接]

新浪微博达人勋

发表于 2020-12-17 09:34:47 | 显示全部楼层 |阅读模式

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

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

x
大家好,我想请问一下,我在wrf-chem的输出结果里找到SO2、NOx、O3的浓度都是四维的,就是有垂直层,一般分析的好像都是地面浓度或者柱浓度,这是怎么通过模式输出结果得到的呀?还是要自己计算,怎么个计算法呢?
有没有了解chem方面的同学帮忙解答一下吖,谢谢!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2020-12-17 17:01:23 | 显示全部楼层
第一层就代表地表

点评

请问楼主怎么计算柱浓度呢?  发表于 2021-5-17 14:44
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-5-17 10:34:38 | 显示全部楼层
楼主,想问下,wrfchem输出的PM2.5和PM10均为零,其他变量都有数值,这是怎么回事?namelist.input 设置如下:
&time_control
run_days                            = 3,
run_hours                           = 13,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2013,
start_month                         = 06,
start_day                           = 25,
start_hour                          = 12,
start_minute                        = 00,
start_second                        = 00,
end_year                            = 2013,
end_month                           = 06,
end_day                             = 29,
end_hour                            = 00,
end_minute                          = 00,
end_second                          = 00,
interval_seconds                    = 21600,
input_from_file                     = .true.,
history_interval                    = 60,
frames_per_outfile                  = 1,
restart                             = .false.,
restart_interval                    = 43200,  
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
auxinput5_interval_m                = 60,
io_form_auxinput5                   = 2
auxinput6_interval_m                = 43200,
auxinput6_inname                    ='wrfbiochemi_d01',
io_form_auxinput6                   = 2,
debug_level                         = 0,
auxinput11_interval_s               = 120, 120, 120, 120,
auxinput11_end_h                    = 192, 30, 30, 30,
/

&domains
time_step                           = 60,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 301,
e_sn                                = 301,
e_vert                              = 35,
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 6000,
dy                                  = 6000,
grid_id                             = 1,
parent_id                           = 0,
i_parent_start                      = 1,
j_parent_start                      = 1,
parent_grid_ratio                   = 1,
parent_time_step_ratio              = 1,
feedback                            = 1,
smooth_option                       = 0
/

&physics
mp_physics                          = 10,
progn                               = 1,
ra_lw_physics                       = 1,
ra_sw_physics                       = 1,
radt                                = 10,
sf_sfclay_physics                   = 1,
sf_surface_physics                  = 1,
bl_pbl_physics                      = 1,
bldt                                = 1,
cu_physics                          = 1,
cudt                                = 0,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
num_land_cat                        = 21,
sf_urban_physics                    = 0,
num_urban_layers                    = 1200,
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
cu_rad_feedback                     = .false.
/

&fdda
obs_nudge_opt                       = 0,0,0,0,0
max_obs                             = 150000,
fdda_start                          =   0.,0.,     0.,     0.,     0.
fdda_end                            = 14400.
obs_nudge_wind                      = 1,1,1,0,0
obs_coef_wind                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_temp                      = 1,0,0,0,0
obs_coef_temp                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_mois                      = 1,1,1,1,1
obs_coef_mois                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_rinxy                           = 1080.,360.,90.,30.,180
obs_rinsig                          = 0.1,
obs_twindo                          = 0.6666667,0.6666667,0.6666667,0.6666667,0.6666667,
obs_npfi                            = 10,
obs_ionf                            = 2,2, 2, 2, 2,
obs_idynin                          = 0,
obs_dtramp                          = 40.,
obs_prt_freq                        = 10, 10, 10, 10, 10,
obs_prt_max                         = 10
obs_ipf_errob                       = .true.
obs_ipf_nudob                       = .true.
obs_ipf_in4dob                      = .true.
obs_ipf_init                        = .true.
/

&dynamics
w_damping                           = 1,
diff_opt                            = 1,
km_opt                              = 4,
diff_6th_opt                        = 2,
base_temp                           = 290.,
damp_opt                            = 0,
zdamp                               = 5000.,
dampcoef                            = 0.2,
khdif                               = 0,
kvdif                               = 0,
non_hydrostatic                     = .true.,
moist_adv_opt                       = 1,
scalar_adv_opt                      = 1,
chem_adv_opt                        = 1,
tracer_adv_opt                      = 1,
tke_adv_opt                         = 1,
/

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

&grib2
/

&chem
kemit = 1,
chem_opt = 5,
bioemdt = 30,
photdt = 30.,
chemdt = 10.0,
io_style_emissions = 2,
emiss_inpt_opt = 1,
emiss_opt = 13,
chem_in_opt = 1,
phot_opt = 0,
gas_drydep_opt = 1,
aer_drydep_opt = 1,
bio_emiss_opt = 0,
gas_bc_opt = 1,
gas_ic_opt = 1,
aer_bc_opt = 1,
aer_ic_opt = 1,
gaschem_onoff = 1,
aerchem_onoff = 1,
wetscav_onoff = 0,
cldchem_onoff = 0,
vertmix_onoff = 1,
chem_conv_tr = 0,
conv_tr_wetscav= 1,
conv_tr_aqchem = 1,
biomass_burn_opt = 0,
aer_ra_feedback = 0,
ne_area = 300,
opt_pars_out = 1,
have_bcs_chem = .true.,
/


&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2023-7-12 16:01:33 | 显示全部楼层
请问柱浓度应该怎么算呐
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2024-2-23 11:44:48 | 显示全部楼层
同问,请问楼主知道柱浓度怎么算了吗
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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