用GrADS画剖面的时候，出图有问题

Eegle

以下gs是用来画沿着117°E螺旋度的剖面图，但是在出图的时候 出现了不管我怎么改变t值，画出来的图 都是一样的。请教各位大神，是不是我的循环设置有错，我尝试过把'define vor=hcurl(uwnd.1,vwnd.2)'  'define lx=-(vor*omega.3*10.)/12.64'这两句加在set t后面 但出现了错误 提示说是vor和lx定义出错，麻烦大家帮我看看，不胜感激。
'reinit'
'enable print e:\zhuang\hel\hel_117E1608.gmf'
'sdfopen E:\zhuang\uwnd.2014.nc'
'sdfopen E:\zhuang\vwnd.2014.nc'
'sdfopen E:\zhuang\omega.2014.nc'
'sdfopen E:\zhuang\rhum.2014.nc'
'sdfopen E:\zhuang\air.2014.nc'

'set lev 1000 100'
'define vor=hcurl(uwnd.1,vwnd.2)'
'define lx=-(vor*omega.3*10.)/12.64'

'set grads off'
'set grid off'
'set csmooth on'
'set map 15 1 9'
'set parea 1 10.2 1.5 8'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'

'set lon 117'
'set lat 20 30'
'set lev 1000 100'
'set zlog on'
'set t 663'
'set gxout shaded'
'set cint 0.5'
'set cmin 0'
'd lx*1e+6'
'set ylpos 0 r'

'set gxout contour'
'set cint 0.5'
'set cmin 0'
'set cthick 11'
'cbarn'
'd lx*1e+6'

'q time'
tt=subwrd(result,3)
hour=substr(tt,1,2)
day=substr(tt,4,2)
mon='6'
year='2014'
BT=hour+8      
if(BT>24)
day=day+1
BT=BT-24
endif
'draw title 117E hel 'year'-'mon'-'day' 'BT':00'
   
'print'
'disable print'
'printim e:\zhuang\hel\hel_117E_1608.gif white'
'c'
;

lqouc

define变量前先set t start end
或者直接把你后面的set t 663放在define的前面

Eegle

lqouc 发表于 2014-10-6 17:58
define变量前先set t start end
或者直接把你后面的set t 663放在define的前面

用了第一种方法 问题得到解决 谢谢 另外我想问一下 我之前的错误原理是什么呢 我想下次的时候 避免再出现类似的错误 最后再次万分感谢。
后来 我尝试在define之前加上循环语句：
h=656
while(h<671)
set t 'h''
define.....
define......
.
.
.
.
'print'
'c'
h=h+1
endwhile
'disable print'
;
但提示错误

                               
登录/注册后可看大图

lqouc

Eegle 发表于 2014-10-6 18:52
用了第一种方法 问题得到解决 谢谢 另外我想问一下 我之前的错误原理是什么呢 我想下次的时候 避免再出现 ...

因为在你define变量的时候grads默认的维度是x(start,end),y(start,end),z(z=1),t(t=1)
你没有设置时间维，所以define出来的数据就只有t=1，后面不管你在怎么设置，全都只有t=1时刻的结果。
你用我建议的两种方法应该都可以实现，但是第一种方法比较消耗计算机的内存，因为不管你用哪个时刻，计算机会吧所有时刻都算一遍。
至于你自己改的时候出现的维度的错误是因为curl这个函数必须是对空间有两维的数组运算。

Eegle

lqouc 发表于 2014-10-6 19:12
因为在你define变量的时候grads默认的维度是x(start,end),y(start,end),z(z=1),t(t=1)
你没有设置时间维 ...

谢谢  后来尝试了一下用循环去画图 可能错误的原理跟之前的一样。

eroson

遇到同样的问题，进来学习一下~~

Eegle

eroson 发表于 2014-10-7 11:01
遇到同样的问题，进来学习一下~~

看完之后 问题得到解决了吗？

xibao

我也遇到了同样的问题呢，受教了

xibao

感谢楼主的无私奉献，好人有好报