grads绘图出错

WOT

'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'
i=567
while(i<=678)
'define qu5=uwnd*shum.3/9.8'
'define qv5=vwnd.2*shum.3/9.8'
'define qud5=ave(qu5,t='i',t='i+4')'
'define qvd5=ave(qv5,t='i',t='i+4')'
'define qus=vint(pres.4(lev=0),qud5,300)'
'define qvs=vint(pres.4(lev=0),qvd5,300)'
'define sandu=hdivg(qus,qvs)'
i=i+12
endwhile
'set z 1'
'set parea 1.5 10.0 1.5 8.0'
'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 0 60'
'set gxout shaded'
'set ccols 2'
'd sandu'
'cbarn 1.2 0'
'draw title 1000-300hpa vapor flux and divergence'
'printim E:\GrADS\test\7.png  white x2000 y1800'
'disable print'
以上是我在论坛里找到的脚本，之后进行改写，不知道哪里出错了，求解答！！！

river

'set parea 1.5 10.0 1.5 8.0'   去掉这句试试

WOT

river 发表于 2022-12-26 23:02
'set parea 1.5 10.0 1.5 8.0'   去掉这句试试

还是不行，图像没变化

ayzqs

WOT 发表于 2022-12-27 10:12
还是不行，图像没变化

'set parea 1.5 10.0 1.5 8.0'去掉还是不行吗？那看不出哪里呀

WOT

ayzqs 发表于 2022-12-28 12:37
'set parea 1.5 10.0 1.5 8.0'去掉还是不行吗？那看不出哪里呀

是否还有其它问题呢？

jlmjlm2007

pres.4(lev=0)改为pres.4(lev=1000)试试

WOT

jlmjlm2007 发表于 2022-12-29 22:05
pres.4(lev=0)改为pres.4(lev=1000)试试

还是不行

jlmjlm2007

WOT 发表于 2023-1-4 10:37
还是不行

'define qu=ave(vint(lev(z=1),q*1000*u,300),t=25,t=37)'
'define qv=ave(vint(lev(z=1),q*1000*v,300),t=25,t=37)'
'define quv=hdivg(qu,qv)'
这是我当时写得整层水汽通量散度，总感觉跟你写的括号位置不一样。

WOT

jlmjlm2007 发表于 2023-1-7 17:06
'define qu=ave(vint(lev(z=1),q*1000*u,300),t=25,t=37)'
'define qv=ave(vint(lev(z=1),q*1000*v,300) ...

你好，能求一个完整版的脚本吗？