求助：fwrite读出的数据只有一层

hz92

我要求出nc数据64年以来17层夏季位势高度，但是运行grads以后grd文件，自己写了ctl文件，再画图设置z为1-17的任何层次，画出图和z=1的的图都是一样的，大家帮忙看看怎么回事。
这是grads文件:
'reinit'
'sdfopen E:\ncep\hgt.mon.mean.ls.nc'
'set gxout fwrite'
'set fwrite e:\shgt.grd'
t=6  
while(t<=765)
z=1
while(z<=17)
'set z 'z''
'set lon 0  359'
'set lat -90 90'
'set t 't''
'define shgt=ave(hgt,t='t',t='t+2')'
'd shgt'
z=z+1
endwhile
t=t+12
endwhile
'disable fwrite'
这是我写的ctl文件
dset ^shgt.grd
title shgt-summer
undef -9999
xdef  145 linear  0  2.5
ydef   73 linear -90 2.5
zdef   17 levels
1000 925 850 700 600 500 400 300 250 200  150 100 70 50 30 20 10
tdef  64 linear   00Jun1948    1yr
vars   1
shgt 0 99 summer-hgt
Endvars

xiaoting

试试将ctl中
vars   1
shgt 0 99 summer-hgt
改为
vars   1
shgt 17 99 summer-hgt

hanqianli

我也遇到同样的问题啦。。。苦恼中。。。

虾米go

有解决么？