加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

ncl写入与读取二进制文件

(2014-08-29 08:41:11)
标签:

股票

分类: NCL
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin

fnl = addfile("./uwnd.mon.mean.nc","r")
        vv = short2flt(fnl->uwnd)

        time = fnl->time
        ntime = dimsizes(time)
;        print(ntime)

        level = fnl->level
        nlevel = dimsizes(level)
;        print(nlevel)

        lat = fnl->lat
        nlat = dimsizes(lat)
;        print(nlat)

        lon = fnl->lon  
        nlon = dimsizes(lon)
;        print(nlon)

        system("rm "+"new.bin")
        do i = 0,ntime-793
             do j = 0,nlevel-17
              do k = 0,nlat-71
do m = 0,nlon-144
                                 fbindirwrite("./new.bin",vv(i,j,k,m))
end do
end do
end do
end do 
end
注明:ncl制作的目的是为了画图,在计算方面弱于fortran和matlab.如果处理比较大的数据尽量避开对ncl的使用。建议使用fortran和matlab.


读取:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
a = cbinread("./new.bin",-1,"float")
        print(a)
end

0

阅读 收藏 喜欢 打印举报/Report
前一篇:散度
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有