- 积分
- 203
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-12-10
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本人刚开始接触NCL,属于初级菜鸟级别。目前自己编了一个程序,是将.txt格式的文件转化为.nc格式的。程序如下:
begin
f = addfile("coemission.nc","c")
;creat the file
nlat = 141
nlon = 121
dim_names = (/"lat","lon"/)
dim_sizes = (/nlat,nlon/)
dim_unlimited = (/False,False/)
filedimdef(f,dim_names,dim_sizes,dim_unlimited)
;define dimensions
atts = True
atts@title = "Sample"
atts@conventions = "None"
atts@creation_date = systemfunc("date")
fileattdef(f,atts)
;creat global attributes
filevardef(f,"lat","float","lat")
filevardef(f,"lon","float","lon")
filevardef(f,"emi_co","float",(/"lat","lon"/))
;define file variables
latatts = True
latatts@long_name = "latitude"
latatts@units = "degrees_north"
filevarattdef(f,"lat",latatts)
lonatts = True
lonatts@long_name = "longitude"
lonatts@units = "degrees_east"
filevarattdef(f,"lon",lonatts)
emi_coatts = 0.0
emi_coatts@long_name ="emission of co"
emi_coatts@units = "None"
emi_coatts@_FillValue = 1
filevarattdef(f,"emi_co",emi_coatts)
;define each variable's attributes
lat = fspan(-90,90,nlat)
f->lat = (/lat/)
lon = fspan(-180,180,nlon)
f->lon = (/lon/)
;creat coordinate arrays
;generate a 2D array
z = asciiread("out.2008_D1.txt",(/141,121/),"float")
f->$"emi_co"$ = z
end
但是每次执行的时候都会提示出错
fatal:Could not create (coemission.nc)
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 11 in file txt2nc.ncl
也就是在filedimdef(f,dim_names,dim_sizes,dim_unlimited)附近出错。
我找了很多遍,都没找到错误所在,所以前来向各位求助
|
|