grads简单的绘制风场，但一个层次画出来是一样的..

cycycycycy

我需要画全球1960年12月-2016年11月的矢量风场图（200hpa，150hpa，100hpa），出来的图一个层次上的风场是一样的。。
自己看不出问题来。。。
有问题的gs：
'reinit'
'sdfopen D:\data\uwnd.mon.mean.nc'
'sdfopen D:\data\vwnd.mon.mean.nc'
'set mpdset lowres'
'set mproj latlon'
'set map auto'
'set map 50 1 1'
'draw map'
'set xlopts 1 4 0.22'
'set ylopts 1 4 0.16'
'set x 1 144'
'set y 1 73'
it=0
while(it<=11)
iz=10
while(iz<=12)
'set z 'iz''
'define uuave=ave(uwnd.1,t=156+'it',t=816+'it',12)'
'define vvave=ave(vwnd.2,t=156+'it',t=816+'it',12)'
'set x 1 144'
'set y 1 73'
'set gxout shaded'
'set csmooth on'
'set black -24 24'
'set ccolor 8'
'set cint 5'
'd uuave'
'cbarn'
'set gxout contour'
'set csmooth on'
'set clab masked'
'set ccolor 8'
'set cint 10'
'd uuave'
'set gxout vector'
'set cthick 10'
'set arrowhead 0.1'
'd skip(uuave,5,0);skip(vvave,0,5);mag(uuave,vvave)'
'set font 5'
'draw title MON.'it'-H'iz' ave-uvwind(1961-2016)'
'printim D:\data\wind\uvvector\NO'it'H'iz'uvwindave1961-2016.png white'
'c'
iz=iz+1
endwhile
'c'
it=it+1
endwhile
;
很奇怪，原来我就画纬向风场出的图就是对的
gs：
'reinit'
'sdfopen F:\data\uwnd.mon.mean.nc'
*'set grads off'
*'set grid off'
'set mpdset lowres'
'set mproj latlon'
'set map 50 1 1'
'draw map'
'set xlopts 1 4 0.22'
'set ylopts 1 4 0.16'
'set x 1 144'
'set y 1 73'
iz=10
while(iz<=12)
'set z 'iz''
it=1
while(it<=12)
'define uuave=ave(uwnd,t='it'+155,t='it'+155+12*55,12)'
'set gxout shaded'
'set csmooth on'
'set clevs 30'
'set ccols 0 6'
'd uuave'

'set gxout contour'
'set cterp on'
'set csmooth on'
'd uuave'
'draw title NO.'it'H'iz' ave-uwind(1961-2016)'
'printim F:\data\NO.'it'.H'iz'uwindave1961-2016.png white'
'c'
it=it+1
endwhile
'c'
iz=iz+1
endwhile
；

cycycycycy

好吧，用下面的gs改了就对了。。。。
虽然还是看不出我第一个gs有什么问题。。。想删帖不会删

lqouc

156+'it'最好写成‘156+it’

cycycycycy

lqouc 发表于 2017-4-20 20:40
156+'it'最好写成‘156+it’

我如果都引号起来，会报错

18393810639

看看你们的讨论哦