nc转grd的四维数组用两个while还是有问题

末端

'reinit'

year=1947
while(year<1948)
'set gxout fwrite'
'set fwrite h:\data\soda\t\T_mon_'%year%'.grd'
'sdfopen h:\data\soda\SODASST'%year%'.nc'
'set x 1 720'
'set y 1 330'
i=1
while(i<13)
'set t 'i''
j=1
while(j<41)
'set z 'j''
'd temp'
j=j+1
endwhile
i=i+1
endwhile


'disable fwrite'
'close 1'
say year
year=year+1
endwhile

'reinit'
ctl文件：
dset h:\data\soda\t\T_mon_1947.grd
undef -9.99e+08
title SODA monthly S
xdef 720 linear -179.75 0.5
ydef 330 linear -75.25 0.5
zdef  40 levels 5374 5124 4874 4624 4374 4124 3874 3624 3374         3124 2874 2624 2375 2125 1875 1625 1378 1139 918 729 579 465 381 317 268 229 197 171 148 129 112 96 82 70 57 46 35 25 15 5
tdef  12 linear 00z15JAN1947 1mo
vars 1
temp 40 99 t,z,y,x unit:degC
endvars

无法画出第一层以外的数据

末端

补充问题，如果单选择任意的某一层数据写的话，是可以成功画出图来

藍鶄

楼主 你这个是不是画出来的图都和第一张一样？？

末端

藍鶄 发表于 2014-6-19 15:48
楼主 你这个是不是画出来的图都和第一张一样？？

不是和第一站一样，而是除了第一张之外，其他层都没有数据

藍鶄

末端 发表于 2014-6-19 17:18
不是和第一站一样，而是除了第一张之外，其他层都没有数据

楼主 你在
dset h:\data\soda\t\T_mon_1947.grd
undef -9.99e+08
着两行之间加一句options template
应该就能解决了 解决了的话告诉我下

末端

藍鶄 发表于 2014-6-19 17:32
楼主 你在
dset h:\data\soda\t\T_mon_1947.grd
undef -9.99e+08

终于解决了，可能是电脑的问题，我把电脑重启了以后这个问题就没有了。用options template应该是不行的，我这里只写了一个文件，不是批量的文件。
哎，鬼电脑这奇葩