画高温天数图，出现缺失值，请问怎么处理

zyj0625

想统计各个点的高温天数图，已经产生各点的高温阈值为tline，想统计各个点高温天数在hd中
这个ncl编程的代码
begin
f=addfile("C:\Users\deii\Desktop\tmax.1980.nc", "r")
tmax=f->tmax
fline=addfile("C:\Users\deii\Desktop\tmax.line.nc", "r")
tline=fline->tline
hd=new((/360,720/),integer)
hd!0="lat"
hd!1="lon"
lon_pts=fspan(0.0,359.5,720)
lat_pts=fspan(-89.75,89.75,360)
lon_pts@units="degrees_east"
lat_pts@unist="degrees_north"
hd&lon=lon_pts
hd&lat=lat_pts
do m=0,359
    do n=0,719
        if(ismissing(tmax(0,m,n)))then
            hd(m,n)=0
            continue
        else
          do i=0,365
            if(tmax(i,m,n).ge.tline(m,n))then
                hd(m,n)=hd(m,n)+1
            end if
          end do
        end if
    end do
end do运行后报错为 The result of the conditional expression yields a missing value. NCL can not determine branch, see ismissing function
请问这种情况可以怎么解决？希望可以得到大家的解答！！

dh5002

请问楼主现在还有这个程序吗，有的话可以发我一份吗