各位大神求帮忙~~画多年平均时间序列图，但是不出图，不知循环哪儿出问题

摩羯小鱼虾

{:soso_e132:}{:soso_e132:}我的循环不出图，循环哪儿有问题呢？谢谢各位了{:soso_e112:}
'reinit'
'sdfopen F:\xyx\shum.mon.mean.nc'
'sdfopen F:\xyx\uwnd.mon.mean.nc'
'sdfopen F:\xyx\vwnd.mon.mean.nc'
'sdfopen F:\xyx\air.mon.mean.nc'
'set grid off'
'set grads off'
'enable print F:\xyx\sj\rain_1.gmf'
**横坐标显示时间1-12月
**画850hpa时间序列图
'set cstyle 1'
'set cmark 3'
'set ccolor 1'
'set cthick 6'
'set y 1 '
'set t 1'
'set lev 850'
**计算水汽通量
'define p=lev'
'define e=exp(17.67*(air.4-273.15)/(air.4-29.65))*(shum.1/100)*6.112'
'define q=622*e/(p-0.378*e)'
'define q1=uwnd.2*q/9.8'
'define q2=vwnd.3*q/9.8'
'define r=aave(mag(q1,q2)/1000,lon=78,lon=105,lat=28,lat=40)'
**设置循环把每个月多年平均水汽通量画出来
'define t=1'
while(t<13)
'define rr=ave(r,t=1,t=792,12)'
'd rr'
' t=t+1'
endwhile
'set display white'
'draw xlab time'  
'draw ylab Q'
'draw string 9.22 0.26  g/(s*hPa*cm)'
'print'
'disable print'

youyou

看了半天，感觉你的这个思路好乱：
1、你既然横坐标想要是时间，而你又在时间上做循环，这样子做只能得到12张图
2、你前面的那个set t 1 已经在计算前就把时间定在了第一个月，那么你后面不管你怎么做，你画出来的图应该都是1月份的
估计要动大手术你这个GS，祝你顺利！！

摩羯小鱼虾

youyou 发表于 2015-5-4 10:54
看了半天，感觉你的这个思路好乱：
1、你既然横坐标想要是时间，而你又在时间上做循环，这样子做只能得到1 ...

那我该怎么做呢？我只要一张图~~横坐标是1-12月，做多年逐月平均；{:eb511:}

yyymm

请问楼主解决了吗，我现在也遇到了类似的问题