留言板
- 侠骨柔情
- 楼主，你好，想问下，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,
/
- 小1龙2虾3
- 我也试过gmail的，也是一直提示说要以.edu结尾的邮箱。那我再试试，谢谢哦！
易小凯: 用gmail的就可以，你试一下
- 小1龙2虾3
- 您好，我在帖子看到您有下载过COSMIC数据，能请教一下怎样下载么？我一直注册不了账号，好像是邮箱有限制，必须要以.edu结尾的邮箱。
- youko
- 谢谢你的回复！那请问如果我计算一个指数，想看它的年际和和年代际变化，是先对原始资料场滤出年际、年代际变化，然后再计算这个指数，还是先计算好这个指数，然后对它的时间序列进行滤波呢？谢谢！
易小凯: 你好。按照我的理解，应该是先求出每年的夏季平均，然后根据每年的状况分解出年际和年代际两部分。
- youko
- 您好！我在搜索关于滤波的帖子时看到了您的回复，有个问题向向您请教。我有35年的逐日数据，如果分析某个物理量的夏季平均，想分别分析年代际变化和年际变化，请问我应该先对逐日资料滤波再进行夏季平均，还是先求出35年的夏季平均再滤波呢？谢谢，打扰您了~
- 易小凯
- 一般来说700hpa以下应该都为正值
为什么啊?
- dcx
- 我在用grads计算散度项的时候
'define vor=hcurl(UGRDprs,VGRDprs)'
'define div=hdivg(UGRDprs,VGRDprs)'
'define f=2*7.292*sin(lat*3.1416/180)*1e-5'
'define s=-(vor+f)*div*1e9'
然后对每一层求区域平均，一般来说700hpa以下应该都为正值，但是我的结果850hpa-700hpa却都为负值，不知道哪里出了问题
- 易小凯
- 这个问题我还真没有遇到过，你再去问问其他人吧
- 侠骨柔情
- 我知道问题出在哪里了，如果这个数据是我直接从网站上下载的，我就可以打开，如果是我先下到服务器上，在传到电脑上，就会出问题。请问这怎么解决？ 给你频频发了这么多，实在打扰了，谢谢！
- 侠骨柔情
- 我其实知道你的，我也是南信大的，6班的，我和5班的杨凤婷认识。
- 侠骨柔情
- 经常问你，真是麻烦了。
易小凯: 你指的是能打开ncep的数据？
- 侠骨柔情
- 我今天换了 matlab 2014a ,还是同样的错误。eCMWF 的数据应当怎么打开呀？
易小凯: 你指的是能打开ncep的数据？
- 侠骨柔情
- 是WRF模式跑出来的数据。现在打不开的数据是ECMWF 上下载的netcdf
易小凯: 你指的是能打开ncep的数据？
- 易小凯
- 你指的是能打开ncep的数据？
- 侠骨柔情
- 相同的错误，难道是我版本有问题？ 可是其他nc 数据是可以打开的
易小凯: 试试用ncread
- 易小凯
- 试试用ncread
- 侠骨柔情
- 大牛啊！想问个问题：我从ECMWF上下载的nc 数据用MATLAB （12a）打开出错，命令为：ncid = netcdf.open('05123u.nc','NC_NOWRITE') 就是按照你说的，还是出错
错误为：
Error using netcdflib
The NetCDF library encountered an error during execution of 'open' function -
'Invalid argument (22)'.
Error in netcdf.open (line 60)
[varargout{:}] = netcdflib ( 'open', filename, varargin{1} );
麻烦帮忙看看！万分感谢！
- 易小凯
- ncid = netcdf.open('HGT.nc','NC_NOWRITE')
注意后面的大小写及格式
- 侠骨柔情
- 不好意思，要再次向你求助了，我从ECMWF 下netcdf的数据用matlab打不开。 ncdisp 命令出现错误。
易小凯: 你用啥软件？
ncdisp ('20131212u.nc','nowrite')
Error using netcdflib
The NetCDF library encountered an error during execution of 'open' function -
'Invalid argument (22)'.
Error in netcdf.open (line 60)
[varargout{:}] = netcdflib ( 'open', filename, varargin{1} );
Error in internal.matlab.imagesci.nc/openToRead (line 1249)
this.ncRootid = netcdf.open(this.Filename,'NOWRITE');
Error in internal.matlab.imagesci.nc (line 122)
this.openToRead();
Error in ncdisp (line 51)
ncObj = internal.matlab.imagesci.nc(ncFile);
- 侠骨柔情
- 谢谢你耐心的回复，那netcdf需要解码吗？
易小凯: 但grib数据解码比较麻烦