把nc转换成dat时出现的维度问题

羊羊羊

       我在做异常场时发现一个问题，一直没想通，我先用nc做完平均态，写成dat后配上相应的ctl，结果打开ctl在同一个时间内，无论怎么设置z维，图都不变，只有改变时间维才会变，是什么原因呢？
       下面给出我的gs(内含ctl)
'reinit'
'set gxout fwrite'
'set fwrite f:/tbbdata/uwnd-5.dat'
it=1979
i=1
while(it<2014)
'sdfopen f:/data/uwnd.'it'.nc'
'set dfile 'i''
'set z 1 17'
'set x 1 144'
'set y 1 73'
if(math_mod(i,4)=0)
'define a=ave(uwnd,t=122,t=152)'
while(iz<=17)
'd a'
iz=iz+1
endwhile
endif
if(math_mod(i,4)!=0)
'define a=ave(uwnd,t=121,t=151)'
while(iz<=17)
'd a'
iz=iz+1
endwhile
endif
it=it+1
i=i+1
endwhile
'disable fwrite'

fname = 'f:/tbbdata/uwnd-5.ctl'
rc = write(fname,'dset f:\tbbdata\uwnd-5.dat')
rc = write(fname,'undef -9.99E+8')
rc = write(fname,'title u')
rc = write(fname,'xdef 144 linear 0 2.5')
rc = write(fname,'ydef 73 linear -90 2.5')
rc = write(fname,'zdef 17 levels 1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30 20 10')
rc = write(fname,'tdef 35  linear MAY1979 1yr')
rc = write(fname,'vars 1')
rc = write(fname,'u 0 99 Time Mask')
rc = write(fname,'endvars')
rc = close(fname)

羊羊羊

操作如图，看到两个层次完全重合，线条变白

river

你去参考一下吧http://bbs.06climate.com/forum.php?mod=viewthread&tid=14459