求助：绘制整层水汽通量和散度更改时间图像无变化

WOT

我想用的时间范围为1982年5月到2020年9月，更改的部分为橙色部分，更改为
i=413
while(i<=872)
更改前后图像无变化，求大神帮我看一下哪里出问题了


'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=161
while(i<=801)
'define qu5=uwnd*shum.3/9.8'
'define qv5=vwnd.2*shum.3/9.8'
'define qud5=ave(qu5,t='i',t='i+3')'
'define qvd5=ave(qv5,t='i',t='i+3')'
'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 40 120'
'set lat 20 60'
'set gxout shaded'
'set ccols 2'
'd sandu'
'set gxout vector'
'd qus;qvs'
'cbarn 1.2 0'
'draw title 1000-300hpa vapor flux and divergence'
'printim E:\GrADS\test\55.png  white x2000 y1800'
'disable print'

river

你自己算算这应该是几张图？你画图的命令又没有循环起来，所有的的图名字都是 55.png ，最终都被最后一张图覆盖了

WOT

river 发表于 2023-1-7 23:09
你自己算算这应该是几张图？你画图的命令又没有循环起来，所有的的图名字都是 55.png ，最终都被最后一张图 ...

我只想画一张图，1000hpa-300hpa整层的水汽通量和散度，所以我应该怎么修改呢

river

WOT 发表于 2023-1-8 09:27
我只想画一张图，1000hpa-300hpa整层的水汽通量和散度，所以我应该怎么修改呢

把你要的时间固定

栤唥de嶶笑

循环体内，逐月计算整层水汽通量和散度，再定义三个累加变量sumsandu=sumsandu + sandu；sumqus = sumqus + qus；sumqvs = sumqvs + qvs。

循环结束后，看累加了多少次，再对上面三个变量取平均即可

WOT

栤唥de嶶笑 发表于 2023-1-9 16:55
循环体内，逐月计算整层水汽通量和散度，再定义三个累加变量sumsandu=sumsandu + sandu；sumqus = sumqus + ...

您好，我更改了脚本，您再看一下对吗？
'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=864,12)'
'define qvd5=ave(qv5,t+0,t=864,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=864,12)'
'define qvd6=ave(qv6,t+0,t=864,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=864,12)'
'define qvd7=ave(qv7,t+0,t=864,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=864,12)'
'define qvd8=ave(qv8,t+0,t=864,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(lev=1000),qud,300)'
'define qvs=vint(pres.4(lev=1000),qvd,300)'

'define shuiqi=mag(qus,qvs)'

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


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

'set gxout shaded'
'd shuiqi'
'set gxout barb'
'set ccolor 1'
'set gxout vector'
'd qus;qvs'
'cbarn 1.2 0'

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

'printim E:\GrADS\test\44.png  white x2000 y1800'
'disable print'