- 积分
- 392
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-2-15
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 MARZOLINI 于 2018-2-26 10:19 编辑
画功率谱图错误,不知道错在哪里,求大神帮忙看看,谢谢!
原始资料是中国地区1995-2014年共20年的夏季降水日资料,然后我在中国某一个区域上对这20年夏季日降水求区域权重平均,这步我是通过grads求得。脚本如下:'reinit''open c:\data1\ldobs\summer_precip.ctl'
'set lon 72 135.5'
'set lat 18 53.5'
'set t 1 1840'
'define regavg=aave(precip,lon=107,lon=125,lat=18,lat=30)'
'set fwrite c:\data1\ldobs\summer_regavg1.grd'
'set gxout fwrite'
'set t 1 1840'
'd regavg'
'disable fwrite'
;
下面是summer_regavg1.ctl文件:
DSET C:/data1/ldobs/summer_regavg1.grd
UNDEF -9.99e8
TITLE China Hourly Merged Precipitation Analysis
*
xdef 1 levels 1
*
ydef 1 levels 1
*
ZDEF 1 levels 1
TDEF 1840 LINEAR 00Z01jun1995 1dy
*
VARS 1
regavg 0 99 CH01 combined analysis (mm/Hour)
ENDVARS
接着我用ncl画功率谱的图,但是图画出来特别的奇怪,找不出来自己错在哪了,希望哪位大神能帮我看看,谢谢。
ncl脚本如下:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
preci1=fbindirread("C:/data1/ldobs/summer_regavg1.grd", 0,1840,"float")
preci1!0="t"
sdof_1=specx_anal(preci1,0,2,0.10)
splt_1=specx_ci(sdof_1,0.05,0.95)
printVarSummary(splt_1)
wks=gsn_open_wks("eps","spec")
r=True
r@gsnDraw = False
r@gsnFrame = False
r@tiYAxisString = "Variance/frq_interval"
f = sdof_1@frq
p = 1/f ; *highly non-linear*
p!0 = "f"
p&f = f
p@long_name = "period"
p@units = "day"
r@tiXAxisString = "Period (day)" ; xaxis
r@tiYAxisString = "Variance/freq"
plot= gsn_csm_xy(wks,p, splt_1,r)
draw(wks)
frame(wks)
end
这张图是功率谱图
而且如果选取中国另一区域进行区域权重平均得出的图与这张图一样。
|
-
|