- 积分
- 75
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-7-30
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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/)
|
-
脚本运行现实的错误
|