ncl 缺测值处理问题求助

qizidou

请问大家下，我想从t1中筛选部分数据放到tt1中（t1中有缺测数据），以便画图使用，用的是以下循环语句            
tt1= new((/73,144/),float)
tt1@_FillValue=9.96921e+36
do i=0,72,1
  do j=0,143,1
    if (t1(i,j).ge.2.1.or.t1(i,j).le.-2.1) then
         tt1(i,j)=t1(i,j)
    else if (t(i,j).eq.ismissing(t1)) then
         tt1(i,j)=t1(i,j)
    else if (t1(i,j).gt.-2.1.and.t1(i,j).lt.2.1) then
         tt1(i,j)=tt1@_FillValue
    end if  
   end if
  end if
end do
end do
但总报错The result of the conditional expression yields a missing value. NCL can not determine branch, see ismissing function，我去官网看了也不懂是哪里有问题，还望各位指点！

qizidou

参考了气象家园上别人的程序，改为：
tt1    = new((/73,144/),float)
tt1@_FillValue=9.96921e+36
do i=0,72,1
  do j=0,143,1
    if (ismissing(t1(i,j))) then  
       tt1(i,j)=tt1@_FillValue
     else  if((.not.ismissing(t1(i,j))).and.(t1(i,j).ge.2.1.or.t1(i,j).le.-2.1)) then
       tt1(i,j)=t1(i,j)
     end if
    end if
  end do
end do
就不会再报错了，具体原理还不是太清楚，供大家参考~

爱知天气

qizidou 发表于 2018-8-24 14:34
参考了气象家园上别人的程序，改为：
tt1    = new((/73,144/),float)
tt1@_FillValue=9.96921e+36

请问楼主这个t1是什么   为什么之前没有定义t1  判断语句就出现了t1

qizidou

爱知天气 发表于 2018-12-5 20:56
请问楼主这个t1是什么   为什么之前没有定义t1  判断语句就出现了t1

t1是原始数据，一开始有说明“从t1中筛选部分数据放到tt1中”

SYutin

诶嘛激动，谢谢楼主分享，也帮我解决问题了

peetree

感谢！及时解决掉问题

四月雨之音

楼主有没有试过,t1@_FillValue然后再用你的原来的脚本，还出错么

mmwj1990

感谢分享，很有用

duanxian

{:eb502:}出现了和楼主一样的错误，按照楼主的方式修改居然对了，非常有用