爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 2798|回复: 0

使用ncl创建nc格式文件求助

[复制链接]

新浪微博达人勋

发表于 2017-9-21 13:00:04 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
ncl入门小白,需要将探空资料(txt格式)整理为nc格式文件,按照网上一篇帖子进行了修改,但是运行还是出现了很多问题,主要集中于输出文件(fout)浮点数据写入方面,请各位大神帮忙解决一下,红色部分为运行出错部分。
filename ="/mnt/g/201001021915.txt"
ncols = numAsciiCol(filename)
nrows = numAsciiRow(filename)
data = asciiread(filename,(/nrows,ncols/),"float")
pres = data(:,3)
tdry = data(:,2)
rh   = data(:,4)
alt  = data(:,0)
;txt file read

system("rm -rf /mnt/g/output.nc")
fout = addfile("/mnt/g/output.nc","c")

setfileoption(fout,"DefineMode",True)

fileAtt                = True
fileAtt@title          = "NetCDF output using method 2"
fileAtt@creation_date  = systemfunc("date")
fileattdef(fout,fileAtt)

dimNames = (/"pres","tdry","rh","alt"/)
dimSizes = (/nrows,nrows,nrows,nrows/)
dimUnlim = (/False,False,False,False/)
filedimdef(fout,dimNames,dimSizes,dimUnlim)

filevardef(fout,"pres",typeof(pres),getvardims(pres))
filevardef(fout,"tdry",typeof(tdry),getvardims(tdry))
filevardef(fout,"rh",typeof(rh),getvardims(rh))
filevardef(fout,"alt",typeof(alt),getvardims(alt))


filevarattdef(fout,"pres",pres)
filevarattdef(fout,"tdry",tdry)
filevarattdef(fout,"rh",rh)
filevarattdef(fout,"alt",alt)

fout->pres = (/pres/)
fout->tdry = (/tdry/)
fout->rh   = (/rh/)
fout->alt  = (/alt/)


脚本运行现实的错误

脚本运行现实的错误
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表