NCL读取EC新资料问题

anlizheng

EC新资料是0.125度分辨率的，我用NCL读取转码后的GRIB文件总是出错，提示Segmentation fault (core dumped)，一开始在cygwin平台下运行的，以为内存不够，现在换到FEDORA系统下还是这么提示，请问是什么原因啊
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"   
;load "./myfunctionec.ncl"
begin
dir = "./011712-convert/"
filelist= "W_NAFP_C_ECMF_20150117174525_P_C1D01171200011712001.bin"
in = addfile(dir+filelist+".grib","r")
t2 = in ->2T_GDS0_SFC({30:40},{100:110})
;lat2 = in ->g0_lat_0({30:60})
printVarSummary (t2)
print (in)
end

又是那隻貓

版本问题吧

anlizheng

又是那隻貓 发表于 2015-1-18 14:18
版本问题吧

感觉应该不是，EC老文件这么读取没问题的

尽头的尽头

同一台电脑？那内存不还是一样的。试试大型机上运行

anlizheng

尽头的尽头 发表于 2015-1-18 18:02
同一台电脑？那内存不还是一样的。试试大型机上运行

不是同一台，Fedora那台8G的，我要读的数组是1441*2880的，感觉应该够了。

悠幽

Segmentation fault不一定是内存不够，数组溢出也有可能报错，还有其他挺多情况也会报错，再好好查查程序吧，遇到过几次这个错误，很麻烦的。

anlizheng

悠幽 发表于 2015-1-18 20:08
Segmentation fault不一定是内存不够，数组溢出也有可能报错，还有其他挺多情况也会报错，再好好查查程序吧 ...

我把代码贴出来了，只是简单的读取一下变量，麻烦大神们看看

aongyd

你这段代码我自己机器上运行是可以通过的，是不是你ncl的问题？

anlizheng

aongyd 发表于 2015-1-27 16:48
你这段代码我自己机器上运行是可以通过的，是不是你ncl的问题？

弄好了，是新EC资料二级压缩过了，要先转才能用

蓝胖子

anlizheng 发表于 2015-1-28 16:23
弄好了，是新EC资料二级压缩过了，要先转才能用

楼主，这个二级压缩怎么处理的，我是菜鸟