- 积分
- 2537
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-10-21
- 最后登录
- 1970-1-1
|
发表于 2018-9-16 19:54:11
|
显示全部楼层
我这样设置输出之后就是nc数据了 直接用也可以
path1="/home/qusl/data/ncep1/evp/EVA_NCEP1_1948-2017_daily_144.grd"
nrec=0 ;从第一个数据开始读取
dims=(/25550,73,144/) ;数据的维数
evp = fbindirread(path1,nrec,dims,"float");读取二进制文件
evp@_FillValue= 32766
time=ispan(0,25549,1) ;1948-2017每年365天
lat=fspan(-89.5,89.5,73) ;纬度序列
lat@units = "degrees_north" ; Attach units attributes
lon=fspan(0.5,360.5,144)
lon@units = "degrees_east"
evp!0="time"
evp!1="lat"
evp!2="lon"
evp&time = time
evp&lat = lat ; Assign coordinate arrays
evp&lon = lon |
|