astiny的留言板_气象家园_气象人自己的家园

侠骨柔情 2021-5-20 09:46

楼主，你好，想问下，wrfchem输出的PM2.5和PM10均为零，其他变量比如no2、O3都有数值，这是怎么回事？我用的WRF3.9,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,
/

海蓝蓝 2020-12-27 12:55

astiny: 网不好。很显然。

其实网速还可以啊。迅雷下载嗖嗖的

海蓝蓝 2020-10-31 00:59

hi,请教python下载era5的问题。利用官网的下载实例代码retrieve这一种下载方法，我能实现下载，但是下载速度一开始还可以，但几分钟后就龟速了。请问您知道原因呢吗?

帅气的走在路上 2019-4-2 19:24

前辈您好，我用anthro_emis处理过的meic清单，跑wtfchem，pm2.5特别低，个位数。同时，还遇到一些其他的问题，我QQ327039389，期待得到你的指导，感谢。
南京大学建院，研三。

guodp 2018-12-27 18:03

我在运行的时候提示这个错误是怎么回事,谢谢
WRAPIT Version: 120209
COMPILING /home/IncubatorShokuhou/wrfchem_emission_maker/ver02.f90
LINKING
END WRAPIT

Copyright (C) 1995-2015 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.3.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.
At line 78 of file /home/IncubatorShokuhou/wrfchem_emission_maker/ver02.f90 (unit = 15, file = '/home/IncubatorShokuhou/wrfchem_emission_maker/emidata/201601/2016_01__agriculture__BC.asc')
Fortran runtime error: End of file
是不是ncl没有安装好呀

chickai1109 2018-9-23 18:38

您好，我想请问下用extcof55计算AOD时怎么算的？是用extcof55乘以对应的eta层的高度，再把所有层的值加起来么？

weiniex 2018-8-3 12:08

您好，又来麻烦您了。。。我根据WRF官网NCL例子，修改后自己绘制PM2.5等污染源的图，效果很差。不清楚是绘图程序问题，还是生成的数据有问题。请问您有关于绘制污染物浓度分布的NCL吗？感谢感谢

weiniex 2018-7-24 21:04

您好，prep_chem_source生成的backg 44M，convert_emiss.exe后只有7.5M。emissopt由7M变为1.8M，emissfire由9.1M变为3M。这是正常的吗？wrf结果数据偏差大是不是因为这个哈。

weiniex 2018-7-23 20:03

astiny: 是不是机制不对？chemopt应该选301.

我在生成wrfchemi_d01的时候设置为301了，那三个文件也成功生成了。后面运行wrf.exe的时候chem_opt=11。最后用NCL画图的时候数据差很多，例如画一个PM10的图，陆地部分没有值，海面数值看着正常。或者只有地图边界上有数据，这明显不对。应该是方案没设置对的缘故，我再多试试，NCL也刚开始学习当中。

weiniex 2018-7-23 11:08

astiny: 可能版本问题。需要3。6.1以下版本，以上版本该程序不能用。

嗯，是版本问题，我重装了就好了。生成的污染值存在问题不知道是什么问题。。。

weiniex 2018-7-21 20:18

尝试了好多次.bin文件大小看着没问题，convert也出来文件了但是大小不对，wrfchemi_d01只有2.5kb。不知道您的有多大呢？是不是我的chem没有装好，KPP=0的原因吗？

weiniex 2018-7-21 11:20

astiny: 链接: https://pan.baidu.com/s/1NDJcRF5aeqPF_uiCKaVy8A 密码: d6q1

非常感谢您的宝贵资料

，我正卡在.convert_emiss.exe上，很救急。再次表示由衷的感谢！

weiniex 2018-7-18 16:17

你好，请问我用prep_chem_sources提取PM2.5和相关空气污染指标数据。配置.inp文件该如何配置呢？哪里能找到配置参数的说明哈？

weiniex 2018-7-17 18:35

astiny: 照理来说应该不需要。你报关于HDF5库的错，应该不是因为没有HDF5的库。是其他的原因造成的。我以前也遇到过类似问题，但是我忘了自己是怎么解决的了。

已经安好了

。hdf5.mod 那个错误我是重新安装了HDF5：
“FC=gfortran CC=gcc ./configure --prefix=/usr/local/HDF5 --with-zlib=/usr/local/zlib --enable-fortran”比原来多了FC、CC部分。后来再次编译出现edgar_emissions.f90错误，在一个链接找到了答案：http://forum.wrfforum.com/viewtopic.php?f=39&t=9199

weiniex 2018-7-17 16:03

你好哈,请问我在安装prep_chem_source时候报关于HDF5库的错，那些库我是在装WRF的时候装了，是需要重新装一遍吗？

穆斯塔法 2017-10-20 13:03

问一下，你有用过prep_chem_sources吗，我在运行rep_chem_sources_RADM_WRF_FIM.exe的时候一直处bug:Program received signal SIGSEGV: Segmentation fault - invalid memory reference. 你知道怎么解决吗？前面编译的时候没有问题。

denglin19901030 2017-6-14 08:28

astiny: 应该不用写。。

谢谢啦~

denglin19901030 2017-6-13 12:47

astiny: 找了个学长帮我做了一个源。。。

那请问auxinput10_inname是怎么写的呢？还是可以不用写啊？

denglin19901030 2017-6-11 17:29

astiny: 不要意思。。。完全忘光了。。。而且现在大型机登不上去。。。看不到以前的input。。。

谢谢回复~那想请问下后面你还是用的ANTHRO_emis来编辑人为源的吗？我看你的帖子后面因为单位的问题说换了一个源，想问下是什么呢？谢谢啦

denglin19901030 2017-6-11 12:24

你好，想请问下anthro_emis的input的数据格式应该是怎样的呢？我看readMe里说只要有一个一维的lon&lat,有一个data和datasec以及排放源的变量就可以了，但是我转的时候一直会提示要读取next_flnm，意思是每个小时都要写一个input文件吗？

astiny的个人空间 http://bbs.06climate.com/?28025 [收藏] [复制] [分享] [RSS]

留言板

astiny