关于差值循环的问题

aaaaaz

我用NCEP的2m温度场的再分析资料（192*94，不是等距的），进行差值，差成2.5*2.5的数据。但是我的gs运行之后只有第一个时间的数据，可能是时间循环有问题。我的gs文件如下，请大神给指导一下是哪里的错误。另外，如果大家有现成的2.5*2.5时间为1948-2013年的2m温度数据也请给我一下吧。十分感谢。
'reinit'
'sdfopen d:\lw\21948\air.2m.mon.mean.nc'
'set gxout fwrite'
'set sdfwrite d:\lw\hl\air2m.nc'
ti=796
j=1
while (j<=ti)
'set lon 0 360'
'set lat -90 90'
'set z 1'
'set t 'j
'define tsh=re(air,2.5)'
say j
'sdfwrite tsh'
j=j+1
endwhile

'disable fwrite'
'reinit'