ncl过多循环会使得fatal产生咩

四月雨之音

begin
;----------------------------------------  GPCP -----------------------------
nnlat = 180
nnlon = 360
f1 = addfile("/home/yangsong3/data-observation/PRECPT/gpcp_1dd_v1.2_p1d.19961001-20150531.daily.nc","r")
tc_date = ut_calendar(f1&time,-2)
print(tc_date)
pStart  = 20000101
pEnd    = 20141231
tstart  = ind(tc_date.eq.pStart)
tend    = ind(tc_date.eq.pEnd)
temp = f1->precip(tstart:tend,::-1,:)
gpcp  = temp(0:5474,:,:)
mmm=0
do tt = 0,5478
if(tt.ne.59.and.tt.ne.1520.and.tt.ne.2981.and.tt.ne.4442)
gpcp(mmm,:,:)  = temp(tt,:,:)
mmm=mmm+1
end if
end do
set_default_fillvalue("double",0)
do tt=0,5474
  do jj=0,179
     do ii=0,359
     if(gpcp(tt,jj,ii).eq.0)
     mmm=mmm-1
     end if
    end do
  end do
end do
print(mmm)
我想大概了解一下缺测的数据多不多，不过运行脚本以后说错误如图，想问问大伙是啥造成的？

藕茭白空心菜

这个循环不算多呀

四月雨之音

藕茭白空心菜 发表于 2016-10-27 20:38
这个循环不算多呀

恩恩，我也觉得三个不算多，会跑起来很慢但理论上应该不至于出错

Ning_1f

您好，请问您在哪里下载到了v1.2的nc版本呢？能不能给个链接，非常感谢！

cqlqlq

有缺测值，if ( .not.ismissing(r1(k,j,i)) ) then 去掉缺测值