nc文件多层次数据提取问题

dsidsi

大家帮忙看一下，看看错出在哪里了，本人刚学，看了很久也没发现错在哪里。
在提取1951年-2015年1-6层次6-8月的经向风、纬向风和比湿的年平均值，用如下程序提取出来的数值与单层提取出来的数值完全对不上
'reinit'
'sdfopen E:\temp\NCEP\uwnd.mon.mean.nc'
'sdfopen E:\temp\NCEP\vwnd.mon.mean.nc'
'sdfopen E:\temp\NCEP\shum.mon.mean.nc'
'set gxout fwrite'
'set fwrite E:\temp\shum\uvq1.dat'
'set y 1 73'
'set x 1 144'


iz=1
while (iz<=6)
'set z 'iz''
i=4
while (i<=68)
aa=12*i-6
ab=12*i-4
'define um=ave(uwnd,t='aa',t='ab')'
'd um'
i=i+1
endwhile
iz=iz+1
endwhile

iz=1
while (iz<=6)
'set z 'iz''
i=4
while (i<=68)
aa=12*i-6
ab=12*i-4
'define vm=ave(vwnd.2,t='aa',t='ab')'
'd vm'
i=i+1
endwhile
iz=iz+1
endwhile


iz=1
while (iz<=6)
'set z 'iz''
i=4
while (i<=68)
aa=12*i-6
ab=12*i-4
'define shm=ave(shum.3,t='aa',t='ab')'
'd shm'
i=i+1
endwhile
iz=iz+1
endwhile



'disable fwrite'
'reinit'