将nc资料转变为*.grd的二进制文件出错

这个也有-

我用的数据是NOAA 2'*2‘的地形数据，没有高度这一维，时间维是1，下面是我编写的导出gs
'reinit'
'xdfopen E:\data2\data\high\high_nc.ctl'
'set gxout fwrite'
'set fwrite E:\data2\data\high\high.grd'
it=1
while(it<=1)
'set t 'it''
'set x 1 10801'
'set y 1 5401'
'set lev 0'
'd z'
  it=it+1
endwhile
'disable fwrite'
'reinit'
下面是我自己编写相对应的ctl
dset     E:\data2\data\high\high.grd
title    monthly mean high
undef    nan
dtype    netcdf
xdef    10801 linear -180 0.033333
ydef    5401 linear -90 0.033333
zdef    1 levels 0 1
tdef    1 linear 00Z01FEB2005 1mo
vars     1
z 1  99  high
endvars
但是打开总有错误：
E:/123.png

这个也有-

不知道怎么上传图片，下面是错误显示的
Error:nc_open failed to open file E:\...
detCDF:Unknoun file format
gadsdf:Couldn't ingest SDF metadata
SET Error:No files open yet
SET Error:No files open yet

这个也有-

还有一点说明不知道有没有关系，这个nc资料用sdfopen打不开，是我自己编写了一个ctl，用xdfopen打开的，然后可以很很好作图。

夏朗的芒果

如果是一些非常规的nc数据，可以考虑用fortran读写变量，这样比较方便。

这个也有-

夏朗的芒果 发表于 2014-12-4 20:53
如果是一些非常规的nc数据，可以考虑用fortran读写变量，这样比较方便。

嗯，fortran用的不熟练，所以想走个捷径。

夏朗的芒果

这个也有- 发表于 2014-12-4 21:49
嗯，fortran用的不熟练，所以想走个捷径。

除了fortran，netcdf还提供了其他的库，比如c,java。如果你会matlab的话，应该更方便，matlab也是可以提取的。

lqouc

ctl写的不对吧，明明都导出grd了，为什么还在dtype写着nc呢？

这个也有-

lqouc 发表于 2014-12-5 09:49
ctl写的不对吧，明明都导出grd了，为什么还在dtype写着nc呢？

是改成grid么？

lqouc

这个也有- 发表于 2014-12-5 20:09
是改成grid么？

直接去掉那一行，还有grads应该不能识别nan作为缺省吧，还有zdef明明就一层level为啥写个0还有1？是多少就写多少。

adviser85guwen

我遇到和你一样的问题。你能用sdfopen把NC数据转成二进制了么？