NCL对经纬度进行循环

什夜0.0

请教各位大佬，我要对一个范围内的olr和风场做滤波，已经对一个点滤波成功了，但我应该对经纬度进行循环来得到整个空间的滤波，怎么用do具体写出循环呢？请求解答，感激不尽！

什夜0.0

顶顶 有没有好心人

雨子哥

同求，怎样实现最快循环遍历经纬度格点呢？

什夜0.0

雨子哥 发表于 2022-4-14 12:02
同求，怎样实现最快循环遍历经纬度格点呢？

这是我对1979-2020年北纬0-50°，东经100-180°滤波的循环程序，在循环之前先筛选一下你需要的范围速度会快很多。
begin

f = addfile("/home/mars/Liubs/lunwen/data/olr.day.mean.nc","r")

    time = (f->time)
    time@units  = "hours since 1800-01-01 00:00:0.0"

; Convert to UTC time.

   time_ut  = cd_calendar(time,0)

; Store return information into more meaningful variables.

    year   = time_ut(:,0)  ; Convert to integer for

; print(year)

   indtime = ind((year.ge.1979).and.(year.le.2020))
   ;print(indtime)
   
   olr = (f->olr(indtime,{0:50},{100:180}))
   ;print(olr)
   xBPF2 = new((/15341,21,33/),float)
;-------------do loop -------------------
do i = 0,32
    do j = 0,20

  olr2  = olr(:, j, i)
  
  ihp     = 2                             ; band pass
  sigma   = 1.0                           ; Lanczos sigma

  nWgt    = 31                          ; loose 100 each end                           
  fca     = 1./60.                       ; start freq
  fcb     = 1./30.                        ; last  freq
  wgt     = filwgts_lanczos (nWgt, ihp, fca, fcb, sigma )
  opt     = 0
  ;printVarSummary(wgt)
  xBPF    = wgt_runave_n_Wrap ( olr2, wgt, opt, 0 )      ; 10-20 day
  xBPF@units = "W/m^2"
  xBPF2(:, j, i) = xBPF

   end do
end do ; the stride is not optional if end < start

;-------------output -------------------
copy_VarCoords(olr(:, :, :),xBPF2)
system("/bin/rm -f "+"30-60olr.nc")
  fout = addfile("30-60olr.nc","c")
  filedimdef(fout,"time",-1,True)
  fout->olr = xBPF2(:,:,:)

  end

雨子哥

谢谢楼主！NCL的处理速度还是很快的，我也已经通过直接写循环实现了。之前提问是因为我以为需要并行运行。