求教fwrite只出来一层数据

fashion

想用fwrite重写各层数据减去前一年的值，gs和ctl文件如下，写出来的只有第一层的数据，其他层次都是0，请各位大神看看怎么解决gs文件
'reinit'
'set mpdset hires'
'set grads off'
'set grid off'
'sdfopen f:\ncdata\hgt.mon.mean.nc'
'set gxout fwrite'
'set fwrite f:\ncdata\hgt.mon.mean.new.dat'
tt=1
while(tt<=427)
'set t 'tt
zz=1
while(zz<=10)
'set x 1 144'
'set y 1 73'
'set z 'zz
'define hgt1=hgt'
'set t 'tt+12
'define hgt2=hgt'
'd hgt2-hgt1'
zz=zz+1
endwhile
tt=tt+1
endwhile
'disable fwrite'
end


ctl文件
dset f:\ncdata\hgt.mon.mean.new.dat
title hgt data
undef -9.99E+33
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 10 levels 1000 925 850 700 600 500 400 300 250 200
tdef  427 linear 00Z01JAN1980 1mo
vars 1
hgt 10 99  hgt
endvars

278803532

'set t 'tt''
'set z 'zz''
......
it=tt+12
'define hgt1=hgt(t='tt')'
'define hgt2=hgt(t='it')'
试试这个，话说，你set那两句是不是少了一个分号啊。。

学而不厌

看来两遍，没发现问题！！楼主解决后，指出问题在哪里吧

是不是这的问题
reinit'
'set mpdset hires'
'set grads off'
'set grid off'
'sdfopen f:\ncdata\hgt.mon.mean.nc'
'set gxout fwrite'
'set fwrite f:\ncdata\hgt.mon.mean.new.dat'
tt=1
while(tt<=427)
zz=1
while(zz<=10)
'set t 'tt
'set x 1 144'
'set y 1 73'
'set z 'zz
'define hgt1=hgt'
'set t 'tt+12
'define hgt2=hgt'
'd hgt2-hgt1'
zz=zz+1
endwhile
tt=tt+1
endwhile
'disable fwrite'
end

fashion

278803532 发表于 2016-4-5 14:12
'set t 'tt''
'set z 'zz''
......

分号没问题的啊，你这样写跟我原来本质是一样的，出来的还是只有第一层的数据

fashion

278803532 发表于 2016-4-5 14:12
'set t 'tt''
'set z 'zz''
......

分号没问题的啊，你这样写跟我原来本质是一样的，出来的还是只有第一层的数据

星雨

如果实在找不出问题，可以先一句一句的来，这样就知道哪里出问题了。

fashion

谢谢大家，3楼已经把问题找出来了