- 积分
- 576
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2022-1-21
- 最后登录
- 1970-1-1
|
GrADS
系统平台: |
|
问题截图: |
|
问题概况: |
使用Lanczos带通滤波器对olr进行10-20天滤波时效果较好,30-60天滤波时未过滤出,想请教大佬这是什么原因呢 |
我看过提问的智慧: |
看过 |
自己思考时长(天): |
7 |
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
请教各位大佬,相同的脚本使用Lanczos带通滤波器对olr进行10-20天滤波时效果较好,但进行30-60天滤波时未过滤出任何时间维度数据,这是什么原因呢。脚本如下:
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
; -------------Input Data-------------
f = addfile("/home/u/data/olr.day.mean.nc","r")
olr = f->olr
; -----------Lanczos Filter-----------
ihp = 2
;1 for highfrequency, 2 for bandpass and 0 for lower pass
sigma = 1.0 ; Lanczos sigma
nWgt = 9 ; loose 200days each end
fca = 1./20. ; start freq
fcb = 1./10. ; last freq
wgt = filwgts_lanczos (nWgt, ihp, fca,fcb, sigma )
qBPF = wgt_runave_Wrap( olr(time|:,lat|:,lon|:), wgt, 0 )
copy_VarCoords(olr(time|:,lat|:,lon|:),qBPF)
;-------------Output Data-------------------
system("/bin/rm -f "+"10-20olr.nc")
fout = addfile("10-20olr.nc","c")
filedimdef(fout,"time",-1,True)
fout->qfil = qBPF(time|:,lat|:,lon|:)
end
|
|