用fwrite读取nc资料出问题了，求指教！！

寄居蟹

求指教。我想用fwrite将nc文件的27层的数据读出来，运行也没有错误。但是运行后grads不是卡住就是直接关闭，读写其它gs文件没有任何问题。文件描述如下dset g:/180405.grd
title unset
undef 9.96921e+36
xdef  144 linear   0   2.5
ydef  96 linear   -90   1.89474
zdef   27 levels  1000  985.1  955.9  903.3  831.0  764.3  639.21  551.39  469.2
    398.7  339.17  288.5  244.7  208.49  176.93  150.3  127.37  108.4  91.3
    78.14  61.05  44.33  29.82  18.05  9.88  4.521  2.107
tdef   14 linear 00Z01MAY1804      7200mn      
VARS   81
........
q        0        -999        q
t        0        -999        t
u        0        -999        u
v        0        -999        v
........
ENDVARS


下面是我的读取资料的gs文件
'reinit'
'sdfopen G:\180405.nc'
'set gxout fwrite'
'set fwrite G:\180405.grd'
tt=1
while(tt<=14)
'set t 'tt
'set x 1 144'
'set y  1 96'
ia=1
while(ia<=27)
'set z 'ia
'd t'
ia=ia+1
endwhile

ib=1
while(ib<=27)
'set z 'ib
'd q'
ib=ib+1
endwhile


tt=tt+1
endwhile
'disable fwrite'
;

请各位高手指导下！！多层次多变量nc文件怎么读。我看了下清风的帖子,如果按照其读写法，81个变量27层，每层都定义的话就要写很多。请问这种多层多变量如何循环？谢谢

lqouc

看看你的nc维度的描述吧

dzxconan

0 360度的时候 0和360度的显示了两次~

river

ctl也有问题，一共提取了4个变量，ctl里变量却写的是11，怎么回事儿啊

isabella

还有你的'set t' tt，应该是'set t 'tt''

寄居蟹

isabella 发表于 2013-6-28 20:31

                               
登录/注册后可看大图

还有你的'set t' tt，应该是'set t 'tt''

这个貌似是木有写错的~~

寄居蟹

river 发表于 2013-6-28 18:04

                               
登录/注册后可看大图

ctl也有问题，一共提取了4个变量，ctl里变量却写的是11，怎么回事儿啊

哦哦，这个是传上来的疏忽，删了一些变量，要不传上来太多了。

寄居蟹

dzxconan 发表于 2013-6-28 17:05

                               
登录/注册后可看大图

0 360度的时候 0和360度的显示了两次~

改为"set lon 0 375.5"?

river

寄居蟹 发表于 2013-6-30 21:23

                               
登录/注册后可看大图

改为"set lon 0 375.5"?

那就经纬度设置按格点来吧，set x 1 144   set  y 1 96来试试吧

寄居蟹

river 发表于 2013-7-1 07:30

                               
登录/注册后可看大图

那就经纬度设置按格点来吧，set x 1 144   set  y 1 96来试试吧

谢谢啦，果然是格点设置的问题