【已解决】grads绘图总有warning出现，但是能画出图来，为什么，图是正确的吗？

shirly

我用grads的tregr函数绘图，time.ctl是我造的时间数据，共60个时次的空间场，每一个时次的所有点的值都等于该时次（如第二时次的所有点的time就是2）。sp.p.ave.ctl是时次为60的空间场
time.ctl如下：
dset d:/lw/dat/season/regr_by_grads/time.dat
title Tim
undef  9.96921e+36
xdef 127 linear 72.75 0.5
ydef 79 linear 16.75 0.5
zdef 1 linear 0 1
tdef 60 linear 00Z16Jan1950 1mo
vars 1
tim 0 99 *tim
endvars
sp.p.ave.ctl如下：
dset d:/lw/dat/season/sp/sp.p.ave.dat
title CRU TS 3.10 china
undef  -9.99e+8
xdef 127 linear 72.75 0.5
ydef 79 linear 16.75 0.5
zdef 1 linear 0 1
tdef 60 linear 00Z16Dec1950 12mo
vars 1
p 0 99 china-p
endvars

画图程序如下：
'reinit'
'open d:/lw/dat/season/sp/sp.p.ave.ctl'
'open d:/lw/dat/season/regr_by_grads/time.ctl'
'set x 1 '
'set y 1 '
'set t 1 60'
'define a=aave(tim.2,lon=72.75,lon=135.75,lat=16.75,lat=55.75)'
'set parea 1 5 4.6 8.2'
'set x 1 127'
'set y 1 79'
'set z 1'
'set t 1 '
'define regr1=tregr(a,p.1,t=1,t=60)'
'set gxout shaded'
'set clevs -0.20 -0.16 -0.12 -0.08 -0.04  -0.02 0  0.02 0.04 0.08 0.12 0.16 0.20'
'set ccols  48 47 46 45 44 43 42 22 23 24 25 26 27 28 '
'cnbasemap regr1*10'
'draw ylab Regre of P in spring'
'printim d:/lw/dat/season/regr_by_grads/tu/regr_p.sp.season.gif white'

上面的是计算春天的时间倾向率的，还有其它三个季节的就不粘了，运行结果可以出图，但是总有很多warning出现（具体如下），图对吗？
Data Request Warning:Request is completely outside file limits
……
regr1*10
NO hardcopy metafile open
All files closed ; all defined objects released;
file:///C:/Users/dell/AppData/Roaming/Tencent/Users/943000942/QQ/WinTemp/RichOle/VVIEY_%25G~WXWEVOH8LNJWT6.jpg

应该说清楚了吧，希望能得到指点~^^

river

应该是两个ctl里面的时间设置不一样引起的吧
tdef 60 linear 00Z16Jan1950 1mo
tdef 60 linear 00Z16Dec1950 12mo
第一个是你自己造的，那后面那个的时间就不能乱改，所以我觉得把第一个时间改成tdef 60 linear 00Z16Dec1950 12mo，和第二个统一就行了

王卡卡在

好！！好！！！！

river

错误提示还是warning啊,楼主贴个完整的啊,整个省略号别人也看不出来啊.这本来是个很好的排除错误的东西，不完整的话看不出来问题

做个霸气的木头

warning有不同，要看你的完整代码才能分清楚，还有楼主发错版块了~

mofangbao

图没出来

shirly

做个霸气的木头 发表于 2013-3-31 09:06

                               
登录/注册后可看大图

warning有不同，要看你的完整代码才能分清楚，还有楼主发错版块了~

完整的代码是四个季节的，如果都贴上来的话比较啰嗦。所以只贴了一个季节。
我现在把整个代码贴上来，再把运行的结果warning贴出来。（time.ctl和sp.p.ave.ctl代码如上贴出）

整个代码：
'reinit'
'open d:/lw/dat/season/sp/sp.p.ave.ctl'
'open d:/lw/dat/season/sum/sum.p.ave.ctl'
'open d:/lw/dat/season/aut/aut.p.ave.ctl'
'open d:/lw/dat/season/win/win.p.ave.ctl'
'open d:/lw/dat/season/regr_by_grads/time.ctl'
'define_colors'
'set grads off'
'set grid off'
'set mpdset hires'
'set mpdset hires'
'set x 1 '
'set y 1 '
'set t 1 60'
'define a=aave(tim.5,lon=72.75,lon=135.75,lat=16.75,lat=55.75)'

*******spring
'set parea 1 5 4.6 8.2'
'set x 1 127'
'set y 1 79'
'set z 1'
'set t 1 '
'define regr1=tregr(a,p.1,t=1,t=60)'
'set gxout shaded'
'set clevs -0.20 -0.16 -0.12 -0.08 -0.04  -0.02 0  0.02 0.04 0.08 0.12 0.16 0.20'
'set ccols  48 47 46 45 44 43 42 22 23 24 25 26 27 28 '
'cnbasemap regr1*10'
*'d regr1*10'
'draw ylab Regre of P in spring'
'printim d:/lw/dat/season/regr_by_grads/tu/regr_p.sp.season.gif white'

*********summer
'set parea 6 10 4.6 8.2'
'set x 1 127'
'set y 1 79'
'set z 1'
'set t 1 '
'define regr2=tregr(a,p.2,t=1,t=60)'
'set gxout shaded'
'set clevs -0.20 -0.16 -0.12 -0.08 -0.04  -0.02 0  0.02 0.04 0.08 0.12 0.16 0.20'
'set ccols  48 47 46 45 44 43 42 22 23 24 25 26 27 28 '
'cnbasemap regr2*10'
'draw ylab  Regre of P in summer'

*********autum
'set parea 1 5 0.6 4.2'
'set lon 72.75 135.75'
'set lat 16.75 55.75'
'set lev 1'
'set t 1 '
'define regr3=tregr(a,p.3,t=1,t=60)'
'set gxout shaded'
'set clevs -0.20 -0.16 -0.12 -0.08 -0.04  -0.02 0  0.02 0.04 0.08 0.12 0.16 0.20'
'set ccols  48 47 46 45 44 43 42 22 23 24 25 26 27 28 '
'cnbasemap regr3*10'
'draw ylab Regr of P in autum'

*******winter
'set parea 6 10 0.6 4.2'
'set lon 72.75 135.75'
'set lat 16.75 55.75'
'set lev 1'
'set t 1 '
'define regr4=tregr(a,p.4,t=1,t=60)'
'set gxout shaded'
'set clevs -0.20 -0.16 -0.12 -0.08 -0.04  -0.02 0  0.02 0.04 0.08 0.12 0.16 0.20'
'set ccols  48 47 46 45 44 43 42 22 23 24 25 26 27 28 '
'cnbasemap regr4*10'
'cbarn'
'draw ylab Regr of P in winter'
'printim d:/lw/dat/season/regr_by_grads/tu/regr_p.x.season.gif white'
'reinit'

运行的j结果（warning）:
ga->run d:\lw\dat\season\regr_by_grads\regr.p.season_tu.gs
All Grads attributes have been reinitialized
Data Request Warning:Request is completely outside file limits
……（共有55行上面的warning,完了就是下面的）
regr1*10
regr2*10
regr3*10
regr4*10
No hardcopy metafile open
All files closed ; all defined objects released;
All Grads attributes have been reinitialized
ga->

shirly

mofangbao 发表于 2013-3-31 13:00

                               
登录/注册后可看大图

图没出来

我之前也发过一篇帖子，图片出不来，我是用QQ截图再发的，是不是贴图的方法不对？

mofangbao

shirly 发表于 2013-3-31 15:23

                               
登录/注册后可看大图

我之前也发过一篇帖子，图片出不来，我是用QQ截图再发的，是不是贴图的方法不对？

QQ截图是在你本地 你要用图片上传的方式传上来

shirly

mofangbao 发表于 2013-3-31 18:45

                               
登录/注册后可看大图

QQ截图是在你本地 你要用图片上传的方式传上来

嗯，谢谢清风大哥

做个霸气的木头

shirly 发表于 2013-3-31 15:21

                               
登录/注册后可看大图

完整的代码是四个季节的，如果都贴上来的话比较啰嗦。所以只贴了一个季节。
我现在把整个代码贴上来，再 ...

惭愧了，木有想到gs这么长。话说warning翻译过来是请求完全超出文件范围，应该是某个设置没有在数据的范围内，既然出现了55条，你从头到尾排查下，看哪里的设置没有在数据的范围内。