ncl如何处理grd数据啊

水霖玲 · 发表于 2014-6-14 11:15:50

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

求教ncl处理grd数据方法。

sun_shine_Xia · 发表于 2014-6-14 17:36:43

水霖玲发表于 2014-6-14 16:48
ncl能直接读取gad啊，怎么读取的啊，求教

  u_file    =path+"U/EC-Uwind(1979-2013)FMA.grd"
  a       = fbindirread(u_file,0,(/35,14,73,144/),"float")
像这种啊~然后如果你需要给变量附上属性的话再另外定义就行了！比如说：
   time=ispan(1979,2013,1)
      lat=fspan(-90,90,73)
      lon=fspan(0,360,144)
      time@units="years since 1979-11-01 00:00:0.0"
      lat@units="degrees_north"
      lon@units="degrees_east"
   u=a(:,10,:,:)
      u!0="time"
      u!1="lat"
      u!2="lon"
      u&time=time
      u&lat=lat
      u&lon=lon
上面这个就没有定义Z方向了，你可以自己改一改

sun_shine_Xia · 发表于 2014-6-14 11:44:58

直接读进去不就好了么，想要咋处理？

水霖玲 · 发表于 2014-6-14 16:48:18

sun_shine_Xia 发表于 2014-6-14 11:44
直接读进去不就好了么，想要咋处理？

ncl能直接读取gad啊，怎么读取的啊，求教

水霖玲 · 发表于 2014-6-14 16:48:37

sun_shine_Xia 发表于 2014-6-14 11:44
直接读进去不就好了么，想要咋处理？

ncl能直接读取grd啊，怎么读取的啊，求教

水霖玲 · 发表于 2014-6-14 16:48:50

sun_shine_Xia 发表于 2014-6-14 11:44
直接读进去不就好了么，想要咋处理？

ncl能直接读取grd啊，怎么读取的啊，求教

JensenCJ · 发表于 2014-6-14 17:26:39

水霖玲发表于 2014-6-14 16:48
ncl能直接读取grd啊，怎么读取的啊，求教

你这数据是gfs的还是ec还是fnl的呀，你是想问数据的读取与赋值吗

水霖玲 · 发表于 2014-6-20 15:23:21

sun_shine_Xia 发表于 2014-6-14 17:36
u_file =path+"U/EC-Uwind(1979-2013)FMA.grd"
a = fbind ...

谢谢。我大致搞定了

King@Control · 发表于 2014-6-22 12:01:08

不知楼主的脚本可否share下。。看看如何处理grd

水霖玲 · 发表于 2014-6-25 13:52:29

King@Control 发表于 2014-6-22 12:01
不知楼主的脚本可否share下。。看看如何处理grd

begin
path = "$NCARG_ROOT/lib/ncarg/data/cdf/rDjuping.grd"
f=fbindirread(path,0,-1,"float")
;names=getfilevarnames(f)
;print(f)
f!0="time"
f&time=ispan(1960,2011,1)

y1=runave_n(f,9,0,0)
我这是一维的

		自动登录	找回密码
密码			立即注册