NCL流函数滤波

辣辣辣辣辣椒粉

想对流函数做30-60d的带通滤波来着，流函数中没有缺测值，第0维是时间维数，但是结果迷之bug，滤波结果全是1e+20，看了好久不知道为啥会这样，求指点
time = cd_calendar(u&time,0)
year = toint(time(:,0))
mon = toint(time(:,1))
day  = toint(time(:,2))
ddd=day_of_year (year,mon,day)
yyyyddd=year*1000+ddd

copy_VarCoords(u,sf  )
copy_VarCoords(u,vort)
sf_anom  =new(dimsizes(u),typeof(u),u@_FillValue)
vort_anom=new(dimsizes(u),typeof(u),u@_FillValue)
do i=0,2
  sf_clmday  :=clmDayTLL     (sf(:,i,:,:),  yyyyddd)
  vort_clmday:=clmDayTLL     (vort(:,i,:,:),yyyyddd)
  sf_anom(:,i,:,:)    =calcDayAnomTLL(sf(:,i,:,:),  yyyyddd,sf_clmday)
  vort_anom(:,i,:,:)  =calcDayAnomTLL(vort(:,i,:,:),yyyyddd,sf_clmday)
end do
print("ok2 anom")
;3 butterworth of 30-60d
ca=60
cb=30
dt=1
fca=dt/ca
fcb=dt/cb
opt   = True        ; options to set
opt@return_envelope = True ; time series of filtered and envelope values
sf_bf  =bw_bandpass_filter(sf_anom(0:354,0,1:71,:)  ,fca,fcb,True,0)
vort_bf=bw_bandpass_filter(vort_anom(:,:,1:71,:),fca,fcb,True,0)

yeah...

请问楼主后来滤波结果怎么样啊  我也要做流函数滤波 可是不知道怎么搞完全

辣辣辣辣辣椒粉

yeah... 发表于 2017-4-8 18:52
请问楼主后来滤波结果怎么样啊  我也要做流函数滤波 可是不知道怎么搞完全

求流函数的时候纬度要按照从南到北放就正常了，我是之前求距平的时候除了点问题

linhaiqiu123

学习了，大侠