grads绘制整层水汽通量出错了

WOT

我使用的是这个脚本https://bbs.06climate.com/forum. ... 2300&fromuid=133783绘制整层水汽通量，自己简单修改了一下，但是出错了，只能画出右上角一小部分的水汽通量和散度，求求大佬们帮我看一下哪里出错了
以下是我改写后的脚本
'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 z 1 8'
*'set t 161'
*'set t  801'
'define qud5=ave(uwnd*shum.3/9.8,t=1,t=859)'
'define qvd5=ave(vwnd.2*shum.3/9.8,t=1,t=859)'
*statr integral from pres
'set z 1'
'define qus=vint(pres.4(lev=0),qud5,300)'
'define qvs=vint(pres.4(lev=0),qvd5,300)'
'define div=hdivg(qus,qvs)'
'set parea 1.5 15 1.5 20'
'set mpdset cnworld cnriver'
'set map 1 1 3'
'set xlopts 1 3 0.2'
'set ylopts 1 3 0.2'
'set xlint 10'
'set ylint 10'
'set lon 60 150'
'set lat 10 60'
'set gxout shaded'
'd div'
'set gxout vector'
'd qus;qvs'
*adjust the colorbar location
'cbarn 1.2 0'
'draw title 1000-300hpa vapor flux and divergence'
'printim E:\GrADS\test\7.png white x2000 y1800'
'disable print'

WOT

{:eb303:}{:eb303:}

烟雨镖局

是不是涉及到了来回切换环境的问题？

WOT

烟雨镖局 发表于 2023-1-6 00:09
是不是涉及到了来回切换环境的问题？

可以详细说一下哪里需要改进吗？

老A

你看看数据的经纬度范围对不对

WOT

老A 发表于 2023-1-11 09:39
你看看数据的经纬度范围对不对

对的。范围没错

老A

从脚本看没什么不对，只有可能是
'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'
这几个数据的经纬度范围或者时间有问题，你全部打开看看数据的经纬度范围对不对

WOT

老A 发表于 2023-1-13 08:55
从脚本看没什么不对，只有可能是
'sdfopen E:\GrADS\testwww%uwnd.mon.mean.nc'
'sdfopen E:\GrADS\testw ...

好的好的，谢谢