找了一下午没找到，决定发帖问：grads中时间循环的问题

cch513

各位大侠，请求帮助。本人想通过query命令，来得到时间信息，然后用时间信息来命名文件名。
我的gs附上，请大侠们指点迷津。
'reinit'
'sdfopen f:\fnl\nc\uwnd.1973.nc'
'sdfopen f:\fnl\nc\vwnd.1973.nc'
'sdfopen f:\fnl\nc\hgt.1973.nc'
'sdfopen f:\fnl\nc\air.1973.nc'
'set grads off'
'set grid off'
'set mpdset hires cnworld'
'set lon 90 150'
'set lat 10 50'
*
*
'c'
d1=1973
b1=775
while(b1<=843)
'set lev 850'
'set t 'b1''
result= 'query  time'
tt1=sublin(result,1)
t1=subwrd(tt1,3)
say 'time=' tt1
'define w=sqrt(uwnd.1*uwnd.1+vwnd.2*vwnd.2)'
'set gxout shaded'
'set cterp on'
'set csmooth on'
'set cint 4'
'set cmin 12'
'd w'
'run cbarn.gs'
*
'draw title ' t1 ' 850w&h'
'enable print f:\fnl\nc\gmf\'t1'850str-hgt.gmf'
'print'
'disable print'
；

tt1的结果就不对，直接显示是tt1=querytime
把引号去掉也不行，请教呀

lqouc

result= 'query  time'这句
直接写'q time'试试

mofangbao

楼上正解，q了之后结果就已经在result里面，不需要自己赋值，其实还是楼主没有仔细看官方手册的脚本那一章节

cch513

致敬楼上两位大侠！一语中的！非常感谢！直接按照大侠的指点修改了一下，图就出来啦！

cch513

'reinit'
'sdfopen f:\fnl\nc\uwnd.1973.nc'
'sdfopen f:\fnl\nc\vwnd.1973.nc'
'sdfopen f:\fnl\nc\hgt.1973.nc'
'sdfopen f:\fnl\nc\air.1973.nc'
'set grads off'
'set grid off'
'set mpdset hires cnworld'
'set lon 90 150'
'set lat 10 50'
*
*
'c'
d1=1973
b1=775
while(b1<=843)
'set lev 850'
'set t 'b1''
'query  time'
tt1=sublin(result,1)
t1=subwrd(tt1,3)
say 'time=' tt1
'define w=sqrt(uwnd.1*uwnd.1+vwnd.2*vwnd.2)'
'set gxout shaded'
'set cterp on'
'set csmooth on'
'set cint 4'
'set cmin 12'
'd w'
'run cbarn.gs'
*
'draw title ' t1 ' 850w&h'
'enable print f:\fnl\nc\gmf\'t1'850str-hgt.gmf'
'print'
'disable print'
；

cch513

这是正确的，现将他发在这里，说不定会有人用到呢。

cch513

在这里同时提个建议：就是希望发帖的人有始有终，有问题，问题解决了，也希望能把解决问题的方法或方案或细节，最后放在帖子里，方便后来人查看学习。

river

不用多说了，连清风大神都来了

cch513

lqouc 发表于 2016-2-17 16:57
result= 'query  time'这句
直接写'q time'试试

大侠一语中的，问题解决啦，非常感谢！

cch513

mofangbao 发表于 2016-2-17 20:27
楼上正解，q了之后结果就已经在result里面，不需要自己赋值，其实还是楼主没有仔细看官方手册的脚本那一章 ...

感谢大侠赐教，送朵花，略表心意！