关于grads话水汽通量散度距平的问题。

a292295

我在用水汽通量散画距平的时候出现了logic error 8 in gxshad的错误，我看别人出现这个错误是因为海洋或陆地数据其中一个为缺测值，比如海温数据之类的，才会画不出来。但是rhum、air、uwnd、vwnd理论上grads不是能自动跳过缺省值吗？应该能直接用grads画出来啊。希望大神求救！下面是我grads的gs文件的一部分。
这个是分特定年份求6月的水汽通量散度平均值
'define mm=0'
'define aqdiv6= rhum.1*0 '
tt=6
yr=1948
while(tt<=831)
'set t 'tt
*鼎盛年
*if(yr=1951|yr=1953|yr=1958|yr=1963|yr=1965|yr=1969|yr=1972|yr=1977|yr=1979|yr=1982|yr=1987|yr=1991|yr=1997|yr=2002|yr=2004|yr=2006|yr=2009|yr=2015)  
*衰弱年
if(yr=1952|yr=1954|yr=1959|yr=1964|yr=1966|yr=1970|yr=1973|yr=1978|yr=1980|yr=1983|yr=1988|yr=1992|yr=1998|yr=2003|yr=2005|yr=2007|yr=2010|yr=2016)   
*正常年
*if(yr=1955|yr=1956|yr=1957|yr=1960|yr=1961|yr=1962|yr=1967|yr=1968|yr=1971|yr=1974|yr=1975|yr=1976|yr=1981|yr=1984|yr=1985|yr=1986|yr=1989|yr=1990|yr=1993|yr=1994|yr=1995|yr=1996|yr=1999|yr=2000|yr=2001|yr=2008|yr=2011|yr=2012|yr=2013|yr=2014)   
'define p6= lev'
'define e6=(6.112*exp((17.67*(air.4-273.16))/(air.4-29.65)))'
'define qs6=(0.622*e6/(p6-0.378*e6))'
'define q6=qs6*rhum.1/100'
'define qdiv6=hdivg(q6*uwnd.2/9.8,q6*vwnd.3/9.8)'
'aqdiv6 = aqdiv6 + qdiv6'
'mm=mm+1'
endif
tt=tt+12
yr=yr+1
endwhile
'aqdiv6 = aqdiv6/mm'

这个是分类别求66年的水汽通量散度平均值，还没有用气候态。。。。
'define mm=0'
'define aqdiv66= rhum.1*0 '
tt = 42
while(tt<=831)
'set t 'tt
'define p66= lev'
'define e66=(6.112*exp((17.67*(air.4-273.16))/(air.4-29.65)))'
'define qs66=(0.622*e66/(p66-0.378*e66))'
'define q66=qs66*rhum.1/100'
'define qdiv66=hdivg(q66*uwnd.2/9.8,q66*vwnd.3/9.8)'
'define aqdiv66 = aqdiv66 + qdiv66'
tt =tt +12
endwhile
'aqdiv66 = aqdiv66 / 66'



希望大神们指教，我确实有点懵。。

又是那隻貓

计算的时候最好还是用matlab fortran python这样的工具吧

z.biye

请问您这个问题是怎么解决的