爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 2918|回复: 0

Lanczos带通滤波器对olr滤波

[复制链接]

新浪微博达人勋

发表于 2022-2-19 15:02:10 | 显示全部楼层 |阅读模式
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


olr.30-60d.png
olr.mean.png
olr.10-20d.png
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表