关于使用NCL将txt转换为nc格式的一些问题

fanxu

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/）

[url=]R[@88]SY]8VZP%$3TZUAJMV.png[/url] (15.51 KB, 下载次数: 0)

脚本运行现实的错误

                               
登录/注册后可看大图

眼里有星辰

请问您参考的哪篇帖子，能给我个链接吗，小白不知道怎么text转nc

对方的凤飞飞

请问楼主解决了,你这没有加经纬度信息呀