ncl读取grd文件并作图的问题

麦田_smile

首先，我用的grd文件是fortran计算的相关系数，180*89*25，资料的ctl如下：

dset E:\lww\rxg_res\sst_uv_precip\res_sst.grd
undef -9.99e+2
title
xdef  180 linear 0 2
ydef  89 linear -88 2
zdef 1 levels 0 1
tdef 25 linear 01jan1979 1MO
vars 2
st1 1 -999 timeseries rain data
st2 1 -999 timeseries rain data
endvars
;

我现在在学习ncl，所以那这个来作图练习下，我写的ncl脚本如下：

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
begin
path = "E:/lww/rxg_res/sst_uv_precip/res_sst.grd"
dims = (/25,89,180/)
file = fbindirread(path,0,dims,"float")
wks  = gsn_open_wks("x11","01")
ton    = file(12,:,:)
file!1 = "lat"
file!2 = "lon"              
plot = gsn_contour(wks,ton,False)
delete(plot)
end

但运行脚本之后 一直有问题，请教大家一下，截图如下：

dataiyang

file = fbindirread(path,0,dims,"float")
这个中“0”代表什么意思？
另外，如果本身是grb格式的话ncl应该可以直接使用addfile读取数据的！
给个建议哈！

麦田_smile

dataiyang 发表于 2013-11-29 11:46

                               
登录/注册后可看大图

file = fbindirread(path,0,dims,"float")
这个中“0”代表什么意思？
另外，如果本身是grb格式的话ncl应 ...

不是grb的 是grd的

dataiyang

麦田_smile 发表于 2013-11-29 16:55

                               
登录/注册后可看大图

不是grb的 是grd的

是grib格式的吗？

麦田_smile

dataiyang 发表于 2013-11-30 09:51

                               
登录/注册后可看大图

是grib格式的吗？

是grd格式的

sun_shine_Xia

If written by fortran, the file must have been been written to a file created by a fortran "open" statement of the form
    open (..., access="direct", form="unformatted", recl=..., ...)
官网上有这么一句话，不知道是不是跟这个有关~

abc258886

好啊，太感谢了！

黄小仙儿

请问您的问题解决了吗?我也遇到了这个问题

liulinaacc

请问下你这个grd数据有两个变量，你这个是画的那个变量？如何体现的呢？

我就是-猴子

是不是grads格式的不能用ncl 直接读取啊