NCL求散度之后滤波的问题

喵星教授Mamonde

我想给散度做30-60天滤波，但是滤波之后的数据全变成了fillvalue,不知道是为什么。我之前做的涡度是可以做出来的啊，一模一样的程序。
大家帮我看看喽~~
a=addfile("/data2/HOME/Frances/8phases_0113/obss/uwnd_f.nc","r")
  u1000=a->uwnd(:,{1000},{-30:45},:)
  printVarSummary(u1000)

  b=addfile("/data2/HOME/Frances/8phases_0113/obss/vwnd_f.nc","r")
  v1000=b->vwnd(:,{1000},{-30:45},:)
  printVarSummary(v1000)
  
  div = uv2dv_cfd(u1000,v1000,u1000&lat,u1000&lon, 2) ;computes divergence given the u and v components (on a fixed grid)
  printVarSummary(div)
  copy_VarMeta(u1000,div)
  f=addfile("div_nc.nc","c")
  f->div=div

  f=addfile("div_nc.nc","r")
    div=f->div(:,:,:)
  printVarSummary(div)
   ;time=f->time
   lat=f->lat
   lon=f->lon
  
  time   = div&time
  TIME   = cd_calendar(time,0)
  year   = toint(TIME(:,0))
  month  = toint(TIME(:,1))
  day    = toint(TIME(:,2))
  ddd    = day_of_year(year,month,day)
  yyyyddd= year*1000+ddd

  hclmday_div=clmDayTLL(div,yyyyddd)
  
  hclmday_sm_div=smthClmDayTLL(hclmday_div,3)
  
  hanom_div   =calcDayAnomTLL(div,yyyyddd,hclmday_div)

  hanom_sm_div=calcDayAnomTLL(div,yyyyddd,hclmday_sm_div)

  printVarSummary(hanom_sm_div)
  delete(div)
;;;;filter ***********************************
  nWgt=201
  sigma=1.0
  wgt= filwgts_lanczos (nWgt, 2, 1./60., 1./30., sigma)
  div=  wgt_runave_n_Wrap (hanom_sm_div, wgt, 0,0 )
  printVarSummary(div)
  
  ncdf4=addfile("div_3060.nc","c")
  ncdf4->div=div

；；；；***********************
做滤波之前的输出检查都对，就不知道滤波散度为什么就变成fillvalue.......

Moody's_Mood

你好，我想请教一下nwgt的值是怎么确定的?