还是关于合成求降水平均的问题，求大神分析一下。

pb880991

在fortran里求平均太烦，因为缺测值不好办，我是在grads里弄的：
'reinit'
'open D:/lunwen/stationrain/grid.ctl'
'open D:/lunwen/stationrain/stn_rain.ctl'
'enable print D:/lunwen/stationrain/La.gmf'
la='1998 1999 2000 2001 2006 2008'
iii=1
'rainsum=0'
while(iii<=6)
year=subwrd(la,iii)
mon1=(year-1951)*12+6
mon2=(year-1951)*12+7
mon3=(year-1951)*12+8
'rainsum=rainsum+rain(t='mon1')+rain(t='mon2')+rain(t='mon3')'
iii=iii+1
endwhile
'rain(t=1)=rainsum/18'
'set t 1'
'define a=oacres(g,rain.2,1)'
'define b=maskout(a,g-0.5)'
'set mpdset cnscs'
'set grads off'
'set grid off'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set gxout shaded'
'set t 1'
'set lon 70 136'
'set lat 15 55'
'set parea 0 11 1 7.5'
'd b'
'cbarn 1 0'
'draw title La'
'print'
'disable print'
'reinit'
结果出现这种问题：
说我'rain'不是变量或函数名
怎么办

pb880991

而且grads好像卡死了

mofangbao

rain是哪儿来的？站点ctl里面的？站点数据经过插值了吗？如果rain是ctl里面的，难道还能给现有变量重新赋值？楼主对grads脚本的理解有些不对

river

你打开了两个ctl，先不管你的那个累加对不对，但好歹得告诉grads你的rain变量是哪个ctl里面的，显然应该不是第一个ctl里的吧。再有哪有先给二进制资料重新赋值再插值的呢，根本思路就不对

kongfeng0824

是站点数据吗？有多少个站？