关于gs读写dat的问题【已解决】

希尔博拉斯

    小白想问一个问题，用gs把nc或者dat写到另一个dat里面的时候，似乎只能写三维的数据，写了经纬度跟时间就没法写多层次，请问这种情况是不是只能用fortran去解决？

lqouc

显然fwrite可以输出五维。

希尔博拉斯

lqouc 发表于 2014-4-4 20:21
显然fwrite可以输出五维。

真的啊。。为毛我每次写时间经纬度层次的数据会出错，而只要一个经度或者一个层次就没问题啊。。。
'reinit'
'set gxout fwrite'
'sdfopen f:\data\omega.mon.mean.nc'
'set fwrite f:\wind\wnd\wnd4.dat'



'set x 1 144'
'set z 1 12'
'set y 1 73'

i=0
while(i<=64)

'define a=ave(omega,t='12+i*12',t='14+12*i)'
'd a'
i=i+1
endwhile
'disable fwrite'

这样是可以的吗？

lqouc

希尔博拉斯 发表于 2014-4-4 20:24
真的啊。。为毛我每次写时间经纬度层次的数据会出错，而只要一个经度或者一个层次就没问题啊。 ...

不可以，垂直空间要单独循环，至于为什么，你要搞清楚grads的存储方式。一个x*y可以是一个最大的记录，所以输出的时候，x和y可以循环也可以不用循环，但是z，t，var必须循环。

边缘人糊涂乐

举个例子：
'reinit'
'open e:\air\air.1983.nc'
'set gxout fwrite'
'set fwrite e:\air\air.dat'
'set x 1 144'
'set y 1 73'
id=1
while(id<=365)
'set t 'id''
llev=1
while(llev<=12)
'set z 'llev''
'd air'
llev=llev+1
endwhile
id=id+1
endwhile

'disable fwrite'
'reinit'
;

希尔博拉斯

lqouc 发表于 2014-4-4 20:28
不可以，垂直空间要单独循环，至于为什么，你要搞清楚grads的存储方式。一个x*y可以是一个最大的记录，所 ...

哦。这样啊，好的，谢谢了

希尔博拉斯

边缘人糊涂乐 发表于 2014-4-4 20:40
举个例子：
'reinit'
'open e:\air\air.1983.nc'

懂了。谢谢啦

katrina

最近在处理多层次的海温一直出错，非常感谢！