grads中的问题

追风的阳光

以下是我的gs，我已让散度和涡度都大于0，为何最后画出的散涡比仍有正有负呢？
'reinit'
'clear'
'open D:\zssdata\fnl.ctl'
'open D:\zssdata\pre.ctl'
'enable print D:\zssdata\svbitt2.gmf'
'set dfile 1'
'set grads off'
'set grid on'
'set lev 1000 100'
'set t 9 29'
'set xlevs 00 06 12 18 00 06 12 18 00 06 18 00 06 18 00 06 18 00 06 18'
'set ylevs 10 20 30 40 50 60 70 80 90 100'
'define h=100000*hdivg(ugrdprs,vgrdprs)'
'define w=100000*hcurl(ugrdprs,vgrdprs)'
if ((h>=0) & (w>=0))
'define svbi=h/w'
'd svbi'
'define aa=tloop(aave(svbi,lon=117,lon=120,lat=30,lat=32))'
*'define bb=abs(aa)'
'set lat 29'
'set lon 117'
'set t 9 29'
'set lev 900 100'
'set zlog on'
'set xyrev off'
'set gxout contour'
'set ccolor 1'
'set cthick 4'
'set xlevs 00 06 12 18 00 06 12 18 00 06 18 00 06 18 00 06 18 00 06 18'
'set ylevs 975 925 850 700 600 500 400 300 200 100'
'set clopts -1 -1 0.18'
'set clopts -1 -1 0.18'
'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'
'set cint 5'
'set ylopts 1 5 0.2'
'set parea 1 10 1 7.5'
'd aa'
endif

mofangbao

图呢？

追风的阳光

mofangbao 发表于 2013-4-1 17:06

                               
登录/注册后可看大图

图呢？

d:\zssdata\svbitt2.gmf

追风的阳光

这就是原图，还有一部分gs上面没给出

river

追风的阳光 发表于 2013-4-1 19:14

                               
登录/注册后可看大图

图挂了······估计是这个判断语句有问题吧，还有就是楼主的时间空间都是变化的，就直接d svbi出来的是哪个时刻的图都没确定呢

追风的阳光

不好意思，这是整个gs，上面的那个是我做了个测试看看的。
'reinit'
'clear'
'open D:\zssdata\fnl.ctl'
'open D:\zssdata\pre.ctl'
'enable print D:\zssdata\svbitt2.gmf'
'set dfile 1'
'set grads off'
'set grid on'
'set lev 1000 100'
'set t 9 29'
'set xlevs 00 06 12 18 00 06 12 18 00 06 18 00 06 18 00 06 18 00 06 18'
'set ylevs 10 20 30 40 50 60 70 80 90 100'
'define h=100000*hdivg(ugrdprs,vgrdprs)'
'define w=100000*hcurl(ugrdprs,vgrdprs)'
*'set h'h''
*'set w'w''
if (h>=0)
if (w>=0)
'define svbi=h/w'
'define aa=tloop(aave(svbi,lon=117,lon=120,lat=30,lat=32))'
*'define bb=abs(aa)'
endif
endif
'set lat 29'
'set lon 117'
'set t 9 29'
'set lev 900 100'
'set zlog on'
'set xyrev off'
'set gxout contour'
'set ccolor 1'
'set cthick 4'
'set xlevs 00 06 12 18 00 06 12 18 00 06 18 00 06 18 00 06 18 00 06 18'
'set ylevs 975 925 850 700 600 500 400 300 200 100'
'set clopts -1 -1 0.18'
'set clopts -1 -1 0.18'
'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'
'set cint 5'
'set ylopts 1 5 0.2'
'set cint 5'
'set parea 1 10 1 7.5'
'd aa'
'set dfile 2'
'set gxout bar'
'set grads off'
'set bargap 60'
'set baropts outline'
'set ccolor 1'
'set grid on'
'set ylpos 9 l '
'set t 1 6'
'set xlevs 00 00 00 00 00 00'
'set ylevs 0 50 100 150 200 250 300 '
'set lon 0'
'set lat 0'
'set lev 10'
'set t 1 6'
'set clopts -1 -1 0.18'
'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'
'set cint 5'
'set parea 1 10 1 7.5'
'd a'
'print'
'disable print'
*'printim d:\zssdata\sbvitaa.gmf'
;

追风的阳光

river 发表于 2013-4-1 19:25

                               
登录/注册后可看大图

图挂了······估计是这个判断语句有问题吧，还有就是楼主的时间空间都是变化的，就直接d svbi出来的 ...

判断语句有什么问题呢？我的gs在楼下写出了

river

追风的阳光 发表于 2013-4-1 19:48

                               
登录/注册后可看大图

不好意思，这是整个gs，上面的那个是我做了个测试看看的。
'reinit'
'clear'

直接帮你把图帖上来吧

个人感觉要同时满足h>=0和w>=0的判断语句还是有问题啊

追风的阳光

river 发表于 2013-4-1 21:13

                               
登录/注册后可看大图

直接帮你把图帖上来吧

个人感觉要同时满足h>=0和w>=0的判断语句还是有问题啊

恩，h，w组合的其他情况等值线都是不显示的，我再看看if的问题

追风的阳光

mofangbao 发表于 2013-4-1 17:06

                               
登录/注册后可看大图

图呢？

完整的gs和图在楼下，能帮我看看是哪出问题了吗