用grads来提取fnl数据，运行gs文件为什么提示Invalid operand，新手求指教

laonazaici

我跟着兰溪的那个教学一步步来做，ctl等文件都能生成，但是运行那个GS文件的时候就总是报错，求看看问题出在哪里。我是想提取fnl文件里一个定点位置的一个变量的数据。
gs文件如下：
'reinit'
'open e:/fnl_20101013_00_00_c.ctl'
'set gxout fwrite'
'set fwrite e:/test.dat'

'set lat 15'
'set lon 110'

z=1
while(z<=26)
'set z'z
'd TMP1829m'
z=z+1
endwhile

'disable fwrite'
'reinit'
;

laonazaici

对应的ctl文件开头部分如下：
dset e:/fnl_20101013_00_00_c
index e:/fnl_20101013_00_00_c.idx
undef 9.999E+20
title e:/fnl_20101013_00_00_c
*  produced by grib2ctl v0.9.12.5p33k
dtype grib 3
options yrev
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 1 linear 00Z13oct2010 1mo
zdef 26 levels
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10
vars 127
no4LFTXsfc  0 132,1,0  ** surface Best (4-layer) lifted index [K]
no5WAVAprs  0 230,100,500  ** 5-wave geopot. height anomaly [gpm]
no5WAVHprs  0 222,100,500  ** 5-wave geopotential height [gpm]
ABSVprs 26 41,100,0 ** Absolute vorticity
ACPCPsfc  0 63,1,0  ** surface Convective precipitation [kg/m^2]
ALBDOsfc  0 84,1,0  ** surface Albedo [%]
。。。。

laonazaici

运行的时候每做一次循环都会提示错误Invalid operand，求教是什么原因

laonazaici

如果是把GS文件改成如下形式：
'reinit'
'open e:/fnl_20101013_00_00_c.ctl'
'set gxout fwrite'
'set fwrite e:/test.dat'

'set lat 15'
'set lon 110'

zz=1
while(zz<=26)
'set z 'zz
'd TMP1829m'
zz=zz+1
endwhile

'disable fwrite'
'reinit'
；
报错就是
All file closed;all defined objects released;
All DrADs attributes have been reinitialized

river

laonazaici 发表于 2016-11-26 08:38
如果是把GS文件改成如下形式：
'reinit'
'open e:/fnl_20101013_00_00_c.ctl'

那是正常的提示，不是错误

laonazaici

river 发表于 2016-11-26 17:39
那是正常的提示，不是错误

按照这个编写的ctl文件也打不开，求大神看看是不是ctl文件哪里编写错误了。
dset e:/test.dat
undef 9.999E+20
title e:/fnl_20101013_00_00_c
*  produced by grib2ctl v0.9.12.5p33k
ydef 1 linear 15.00000
xdef 1 linear 110.000000
tdef 1 linear 00Z13oct2010 1mo
zdef 26 levels
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10
vars 1
TMP1829m  0 11,103,1829 ** 1829 m Temp. [K]
ENDVARS

river

laonazaici 发表于 2016-12-1 14:17
按照这个编写的ctl文件也打不开，求大神看看是不是ctl文件哪里编写错误了。
dset e:/test.dat
undef 9. ...

你的资料只有一个点？一个时次？但是有很多层？

laonazaici

river 发表于 2016-12-1 21:08
你的资料只有一个点？一个时次？但是有很多层？

资料是一个时刻的fnl资料，但是我想提取的是其中某一个变量在一个点的不同层的值

laonazaici

river 发表于 2016-12-1 21:08
你的资料只有一个点？一个时次？但是有很多层？

上面那个ctl是我跟着兰溪的教学里面，对应我之前的那个gs文件编写的

temperature

有用有用，谢谢楼主