请教NCL Butterworth带通滤波问题

什夜0.0

在ncl官网下载的程序，修改参数如下：
diri  = "/home/u/uwnd."
    fili  = "850.1979-2020.nc"
    fi    = addfile(diri+fili, "r")
    ua    = fi->uwnd            ; ua(time)    ; read from one grid point
    ca    = 60.0        ; band start (longer period)
    cb    = 30.0        ; band end

    fca   = 1.0/ca      ; 'left'  frequency
    fcb   = 1.0/cb      ; 'right' frequency

    dims  = 0           ; 'time' dimension  

    opt   = True        ; options to set
    opt@return_envelope = True ; time series of filtered and envelope values

    ua_bf = bw_bandpass_filter (ua,fca,fcb,opt,dims)       ; (ua,fca,fcb,opt,dims)
    ua_bf!0  =  "a"
    copy_VarMata(ua, ua_bf(0,:,:,:,:))

    ua_bf@long_name = "Band Pass: "+cb+"-"+ca+" day"

    printVarSummary(ua_bf)

将ua维度复制给ua_bf时（倒数第三行）出现错误，想问问应该怎么修改呢

Zhangyuanping

分别输出原始ua的纬度大小看看，再输出ua_bf的看看，赋信息的时候对应大小赋就行了

什夜0.0

Zhangyuanping 发表于 2022-11-15 15:29
分别输出原始ua的纬度大小看看，再输出ua_bf的看看，赋信息的时候对应大小赋就行了

哈哈哈谢谢您 已经解决啦

zhangdapao

我也遇到了，最后发现是原本的Varmeta打错了哈哈哈 打成了Varmata