急：画图出错，求帮忙

fairyes

我是学grads的新手，老师让画个图要求得出中国区域几个时间点pdsi平均值的前后几年的变化曲线，我写的gs文件如下所示：
'sdfopen d:/grads/mada/data.cdf'
'set lon 73.66'
'set lat 4'
'set lev 0 0'
'set t 152 157'
'set e 1 1'
'define mpdsi=aave(pdsi,lon=73.66,lon=135.08,lat=4,lat=53.52)'
'd (mpdsi(t)+mpdsi(t+148)+mpdsi(t+363)+mpdsi(t+383)+mpdsi(t+511)+mpdsi(t+539))/6'
;

显示出错信息为syntax error:invalid dimension expression
           expecting +/-/= operator,found)
           error eccurred at column 2
          display error: invalid expression
          expression =(mpdsi(t)+mpdsi(t+148)+mpdsi(t+363)+mpdsi(t+383)+mpdsi(t+511)+mpdsi(t+539))/6
看不太懂是什么意思，求高手帮忙
ps:怎么截图上来能看清啊，grads是黑底红色出错信息，截图上来都看不清

善人/jw

你那个T应该是做循环来算吧？

fairyes

善人/jw 发表于 2013-1-24 00:15

                               
登录/注册后可看大图

你那个T应该是做循环来算吧？

哦，多谢提醒，我刚学fortran和grads，试了下写循环，
'sdfopen d:/grads/mada/data.cdf'
'set lon 73.66'
'set lat 4'
'set lev 0 0'
'set t 152 157'
'set e 1 1'
'define mpdsi=tloop(aave(pdsi,lon=73.66,lon=135.08,lat=4,lat=53.52))'
while(t<=157)
'set t 152'
'define data=(mpdsi(t)+mpdsi(t+148)+mpdsi(t+363)+mpdsi(t+383)+mpdsi(t+511)+mpdsi(t+539))/6'
t=t+1
endwhile
'set gxout line'
'd data'
;
但还是有问题，出错信息说data not a variable or function name, 那我该怎么设置这个变量呢？求高手指点，谢谢啦~~~

善人/jw

循环写错了，看看手册循环那部分

river

fairyes 发表于 2013-1-24 16:56

                               
登录/注册后可看大图

哦，多谢提醒，我刚学fortran和grads，试了下写循环，
'sdfopen d:/grads/mada/data.cdf'
'set lon 73. ...

楼主刚开始学习就要写这些比较复杂的东西啊，先看看基础的练习一下吧

fairyes

river 发表于 2013-1-26 19:34

                               
登录/注册后可看大图

楼主刚开始学习就要写这些比较复杂的东西啊，先看看基础的练习一下吧

嗯，多谢鼓励，主要是老板当时想看数据分析结果，比较急，我会慢慢学习的