NCL 滤波

鲤鱼么

求助！NCL滤波，10-20天带通，OLR-day-mean，算出来全是NaN.....请问有没有大佬教教我
脚本如下：
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("e:/olrjason.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    =   201                ; 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(lat|:,lon|:,time|:), wgt, 0 )     
  copy_VarCoords(olr(lat|:,lon|:,time|:),qBPF)


;-------------Output Data-------------------
system("/bin/rm -f "+"qbwo.nc")
  fout = addfile("e:/qbwo.nc","c")
  filedimdef(fout,"time",-1,True)
  fout->qfil = qBPF(time|:,lat|:,lon|:)

  end

我还是叫小确幸

如果用的lanczos滤波，首尾会丢失数据的。可以先用长时间序列滤波，之后再挑选你关注的时段。

一大碗年糕

最近刚学习的滤波和ncl的滤波函数记录了一下，不知有没有用哈哈
http://mp.weixin.qq.com/s?__biz= ... 69b4f41066568ef2#rd

鲤鱼么

我还是叫小确幸 发表于 2021-5-25 09:47
如果用的lanczos滤波，首尾会丢失数据的。可以先用长时间序列滤波，之后再挑选你关注的时段。

谢谢你~我参考了NCL官网的定义也是先算了长时间之后选了自己的时段，但就是得不出数据....我再瞅瞅 谢谢哦

鲤鱼么

谢谢各位回复哦...但是我还是搞不懂

鲤鱼么

一大碗年糕 发表于 2021-5-25 09:17
最近刚学习的滤波和ncl的滤波函数记录了一下，不知有没有用哈哈
http://mp.weixin.qq.com/s?__biz=MzkyODI ...

谢谢你哦...但是我这个新手有点看不出脚本有啥问题...是不是需要循环呀

南气斑织逍遥

olr(lat|:,lon|:,time|:) 改为  olr(time|:,lat|:,lon|:),

鲤鱼么

南气斑织逍遥 发表于 2021-5-28 11:58
olr(lat|:,lon|:,time|:) 改为  olr(time|:,lat|:,lon|:),

谢谢~滤波成功啦 和python差不多~

什夜0.0

你好，想问最后怎么滤波成功的呢，我试了一下你的程序，但提示我维度不匹配，求解答，感恩