grads新手寻求帮助

瑾珂

'reinit'
'open K:/ncep/fnl_2017.ctl'
'set gxout contour'
'set ccolor 1'
'set cthick 5'
'set display white color'
'set mpdset cnworld'
'set xlopts -10'
'set ylopts -10'
'set grid off'
'set lon 80 120'
'set lat 10 30'
'set xlint 10'
'set ylint 5'

i=1
while (i<=44)
'set t 'i''
'define season=UGRDprs(lev=850)-UGRDprs(lev=200)'
i=i+1
endwhile




'display ave(season,t=25,t=37,-b)'
'run axis.gs -type b -position i -label on -interval 10 -sinterval 2 -lsize 0.17 -langle 45 -voffset -0.2'
'run axis.gs -type l -position i -label on -interval 5 -sinterval 2 -lsize 0.17'
'draw title seasonwind'
'printim K:/ncep/fnl_seasonwind2.png white'
以上脚本，'display ave(season,t=25,t=37,-b)'无论t=xx怎么变化，为什么出来的图总是一模一样啊？到底是哪里的问题?

男紫汗

数值是一摸一样吗？可能数值有微小差异在图上不容易看出来

guzhong

最致命的在这里
i=1
while (i<=44)
'set t 'i''
'define season=UGRDprs(lev=850)-UGRDprs(lev=200)'
i=i+1
endwhile
这一段语句运行出来的结果是season是t=44时的u风差。
我比较愚笨，提供一个不高明的解决方案：
i=25
j=0
'define s=0'
'define season=0'
while (i<=37)
'set t 'i''
'define ss=season+s'
'define s=season'
'define season=UGRDprs(lev=850)-UGRDprs(lev=200)'

i=i+1
j=j+1
endwhile
'd ss/j'
再根据你需要修改 i 的取值范围

DAISYQH

应该要先清除之前的图像，才能出现下一张，在printim那行命令之后加'c'试试

瑾珂

男紫汗 发表于 2017-12-2 13:15
数值是一摸一样吗？可能数值有微小差异在图上不容易看出来

数值是不一样的

瑾珂

guzhong 发表于 2017-12-2 15:02
最致命的在这里
i=1
while (i

好的，谢谢，我去看看

rgyxtq2009

好的，谢谢咯扑住