请教大家为什么gs文件里写q ctlinfo不能返回ctl信息？

2012daqi

在用nc文件前不都是需要q ctlinfo嘛，之前都是先打开文件然后在grads里q的，然后人工把各个维度的参数算好，进行下面的计算的，今天把这个q ctlinfo写进gs文件以后，发现不能返回ctl信息，为什么呀？有没有什么方法可以在用nc文件时自动读取x，y，z的格点个数的？
这是我写的gs，就这么几行应该不会错吧，运行完在grads里再运行q ctlinfo可以返回ctl信息，说明文件打开没错呀。。。
'reinit'
filename='d:\data\OLR\olr.day.mean.nc'
'sdfopen 'filename''
result='q ctlinfo'
xinfo=sublin(result,5)
xgird=subwrd(xinfo,2)
yinfo=sublin(result,6)
ygird=subwrd(yinfo,2)
zinfo=sublin(result,7)
zgird=subwrd(zinfo,2)
tinfo=sublin(result,8)
tgird=subwrd(tinfo,2)
say 'x=' xgrid

2012daqi

哦，对了，运行完以上程序后，返回的的是 x=xgird，我在想会不会是xgrid没赋值或者xgrid是个字符变量问题？

梦纯

'reinit'
'open summer_oacres_prec.ctl'
'q ctlinfo'
record=sublin(result,4)
rd=subwrd(record,2)
rc=write(tt.txt,rd)
;
简单试了以上gs，完全可以！result='q ctlinfo'，这句改为'q ctlinfo'试试~~

2012daqi

梦纯 发表于 2016-4-22 14:45
'reinit'
'open summer_oacres_prec.ctl'
'q ctlinfo'

问题解决了，太感谢了！！