NCL把变量写入nc文件

nwafuhcz

我想把ERA5的降水资料，1979年全年，一天一个文件，写入一个nc文件当中。首先批量把变量读取出来，然后写入nc。写入nc文件的代码完全是copy《NCL数据处理与绘图基础教程》的，但是试了好多次就是fileattdef(fout, fAtt)这一句有错误，错误提示如下：fatal:Could not create (G:\ERA5\daily\1979.nc)warning:FileSetFileOption: invalid file or format
fatal:["Execute.c":8637]:Execute: Error occurred at or near line 31 in file practice.ncl
有高人的话救救孩子吧，把百度找遍了也不知道是哪里的错误
下面是我的代码：

1. begin
   
2. 
   
3.   list_of_files = systemfunc("ls G:/ERA5/daily/precipitation/1979/era5.precip.197901*.nc")
   
4.   file_list = addfiles(list_of_files, "r")
   
5.   outfile = "G:\ERA5\daily\1979.nc"
   
6.   time = file_list[:]->time
   
7.   lat = file_list[0]->lat
   
8.   lon = file_list[0]->lon
   
9.   tprate = file_list[:]->tprate
   
10.   ;printVarSummary(time)
    
11.   ;printVarSummary(lat)
    
12.   ;printVarSummary(lon)
    
13.   ;printVarSummary(pre)
    
14.   ntim = dimsizes(time)
    
15.   nlat = dimsizes(lat)
    
16.   nlon = dimsizes(lon)
    
17. 
    
18.   pre = tprate
    
19.   pre = tprate*1000*60*60*24
    
20.   pre@units = "mm"
    
21. 
    
22.   fout = addfile(outfile, "c")
    
23.   setfileoption(fout, "DefineMode", True)
    
24.   fAtt = True
    
25.   fAtt@title = "NCL Efficient Approach to netCDF Creation"
    
26.   fAtt@souce_file = file_list
    
27.   fAtt@Conventions = "CF"
    
28.   fAtt@creation_date = systemfunc("date")
    
29.   fAtt@history = "NCL scripts: practice.ncl"
    
30.   fAtt@comment = "Convert variable tprate from pricipitation rate to amount"
    
31.   fileattdef(fout, fAtt)
    
32. 
    
33.   dimNames = (/"time", "lat", "lon"/)
    
34.   dimSizes = (/"-1", "nlat", "nlon"/)
    
35.   dimUnlim = (/True, False, False/)
    
36.   filedimdef(fout,dimNames,dimSizes,dimUnlim)
    
37.   filevardef(fout,"time",typeof(time),getvardims(time))
    
38.   filevardef(fout,"lat",typeof(lat),getvardims(lat))
    
39.   filevardef(fout,"lon",typeof(lon),getvardims(lon))
    
40.   filevardef(fout,"pre",typeof(tprate),getvardims(tprate))
    
41.   filevarattdef(fout,"time",time)
    
42.   filevarattdef(fout,"lat",lat)
    
43.   filevarattdef(fout,"lon",lon)
    
44.   filevarattdef(fout,"pre",pre)
    
45.   setfileoption(fout,"DefineMode",False)
    
46. 
    
47.   fout->pre = (/pre/)
    
48.   fout->time = (/time/)
    
49.   fout->lat = (/lat/)
    
50.   fout->lon = (/lon/)
    
51.   delete(time)
    
52.   delete(lat)
    
53.   delete(lon)
    
54.   delete(tprate)
    
55.   delete(pre)
    
56. 
    
57. 
    
58. end
    
59. 
    

复制代码

nwafuhcz

就是提示第31行出错

nwafuhcz

救救孩子吧

小羊趴在云朵上

ncl写成nc文件：修改出来了嘛？

杨松

你的outfile路径可能错了

来自安卓APP客户端

iovedys

outfile路径的斜杠是不是反了

小魔孤

一般输出文件的创建和写入：
finame_out = "./*****.nc"
system("rm -rf "+finame_out)
fiout = addfile(finame_out,"c")
fiout->u = data_anal
你可以参考下哈，看是不是缺了什么.

nwafuhcz

小羊趴在云朵上 发表于 2021-4-11 23:10
ncl写成nc文件：修改出来了嘛？

没有,,,,,,,,,,

nwafuhcz

iovedys 发表于 2021-4-15 22:06
outfile路径的斜杠是不是反了

改了之后还是一样的错误提示

不知道起啥名字

我产生这个错误然后修正成功。错误就在output路径。修改前：diro = "data03/***/xinli/models/anom/JMA/t2m/"+year+"_xinli/",修改后：diro = "/data03/***/xinli/models/anom/JMA/t2m/"+year+"_xinli/"。错误就在data03前面那个“/”。供参考。