- 积分
- 11125
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-1-9
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
想对流函数做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)
|
|