求助：绘制矢量箭头出错了！！！

WOT

我参考论坛里的帖子（改良了的画风场参考箭头脚本https://bbs.06climate.com/forum. ... 0196&fromuid=133783）画矢量箭头，但是grads报错了，无法执行drawvect。我已经把drawvect.gs放在制定文件夹下了，为什么还会出错呢？


grads报错内容如下：
ga-> run E:\GrADS\test\zzzhengceng.gs
No hardcopy metafile open
All files closed; all defined objects released;
All GrADS attributes have been reinitialized
Unknown command: drawvect
No hardcopy metafile open


脚本部分：
'reinit'
'sdfopen E:\GrADS\testwww\uwnd.mon.mean.nc'
'sdfopen E:\GrADS\testwww\vwnd.mon.mean.nc'
'sdfopen E:\GrADS\testwww\shum.mon.mean.nc'
'sdfopen E:\GrADS\testwww\pres.mon.mean.nc'


'set parea 1.5 10.0 1.5 8.0'
'set map 1 1 3'
'set mpdset cnworld cnriver'


'set xlopts 1 3 0.2'
'set ylopts 1 3 0.2'
'set xlint 10'
'set ylint 10'


'set t 413 869'
'set z 1 8 '
'define qu5=uwnd*shum.3/1000'
'define qv5=vwnd.2*shum.3/1000'
'set t 413'
'set z 1 8 '
'define qud5=ave(qu5,t+0,t=869,12)'
'define qvd5=ave(qv5,t+0,t=869,12)'
'define qdiv5=hdivg(qud5,qvd5)'


'set t 414 870'
'set z 1 8 '
'define qu6=uwnd*shum.3/1000'
'define qv6=vwnd.2*shum.3/1000'
'set t 414'
'set z 1 8 '
'define qud6=ave(qu6,t+0,t=870,12)'
'define qvd6=ave(qv6,t+0,t=870,12)'
'define qdiv6=hdivg(qud6,qvd6)'


'set t 415 871'
'set z 1 8 '
'define qu7=uwnd*shum.3/1000'
'define qv7=vwnd.2*shum.3/1000'
'set t 415'
'set z 1 8 '
'define qud7=ave(qu7,t+0,t=871,12)'
'define qvd7=ave(qv7,t+0,t=871,12)'
'define qdiv7=hdivg(qud7,qvd7)'


'set t 416 872'
'set z 1 8 '
'define qu8=uwnd*shum.3/1000'
'define qv8=vwnd.2*shum.3/1000'
'set t 416'
'set z 1 8 '
'define qud8=ave(qu8,t+0,t=872,12)'
'define qvd8=ave(qv8,t+0,t=872,12)'
'define qdiv8=hdivg(qud8,qvd8)'

'set z 1 8 '
'define qud=(qud5+qud6+qud7+qud8)/4'
'define qvd=(qvd5+qvd6+qvd7+qvd8)/4'
'define qdiv=(qdiv5+qdiv6+qdiv7+qdiv8)/4'


'set z 1'
'define qus=vint(pres.4(z=1),qud,300)'
'define qvs=vint(pres.4(z=1),qvd,300)'


'define shuiqi=mag(qus,qvs)'


*'define sandu=hdivg(qus,qvs)'




'set lon 40 120'
'set lat 20 60'


'set gxout shaded'
'E:\GrADS\test\text7.gs'
'd shuiqi'
'E:\GrADS\test\cbar_matlab 10 1 1'


'drawvect skip(qus,2,2) skip(qvs,2,2)  tl 0.3 10 3 12 1 0.1 1 m/s'


*'set gxout barb'
*'set ccolor 1'
*'set gxout vector'
*'set arrscl 0.3 100'
*'set arrowhead 0.06'


'd qus;qvs'
*'cbarn 1.2 0'


'draw title 1000-300hpa vapor flux and divergence'


'printim E:\GrADS\test\55.png  white x2500 y2000'
'disable print'

四叶草

看报错提示，文件都没打开。直接上全部脚本，部分脚本太迷惑了

WOT

四叶草 发表于 2023-1-17 08:25
看报错提示，文件都没打开。直接上全部脚本，部分脚本太迷惑了

全部脚本已发出，昨天修改了一下，但还是不行，请帮忙看看

四叶草

WOT 发表于 2023-1-17 09:27
全部脚本已发出，昨天修改了一下，但还是不行，请帮忙看看

昨晚修改后报错也是一样？

WOT

四叶草 发表于 2023-1-17 10:25
昨晚修改后报错也是一样？

一样的一样的一样的