利用函数事先定义NetCDF属性实现nc文件的输出

轻拢烟云

想要以NetCDF格式保留滤波后的数据，但在数据处理之后如果像之前合并多年数据那样直接写入，无法保留时间维度信息，然而我还想把这个数据保存好以便后续使用，因此采取了用函数事先定义的方法，在这里记录一下。
读入文件：

1. f1 = addfile(pt+"pres.sfc.1979-2019.nc", "r")
   
2. 
   
3. pres = f1->pres/100
   
4. printVarSummary(pres)
   
5. lat = f1->lat
   
6. lon = f1->lon
   
7. time = f1->time
   
8. 
   
9. nlat = dimsizes(lat)
   
10. nlon = dimsizes(lon)

复制代码

进行一些数据处理（进行7-90天带通滤波）

1. ihp = 2                                         ; band pass
   
2. sigma = 1.0                                     ; Lanczos sigma
   
3. nWgt = 201                                      ; loose 100 each end
   
4. fca = 1./90.                                   ; start freq
   
5. fcb = 1./7.                                    ; last freq
   
6. wgt = filwgts_lanczos (nWgt, ihp, fca, fcb, sigma )
   
7.                
   
8. plcz = wgt_runave_n_Wrap(pres, wgt, 0, 0)                ; 7-90 day

复制代码

输出变量

1.    system("rm -rf ls F:/data/NCAR-Reanalysis1/pres.sfc.1979-2019_lb.nc")    ;删除可能存在的文件
   
2.    fout = addfile("F:/data/NCAR-Reanalysis1/pres.sfc.1979-2019_lb.nc","c")  ;创建文件
   
3. 
   
4.    setfileoption(fout, "DefineMode", True)
   
5. 
   
6.    ;定义文件的全局属性
   
7.    fAtt   = True
   
8.    fAtt@title  = "NetCDF output using method 2"
   
9.    fAtt@creation_date = systemfunc("date")
   
10.    fileattdef(fout, fAtt)
    
11. 
    
12.    ;定义坐标变量
    
13.    dimNames = (/ "time","lat","lon"/)       ;命名维
    
14.    dimSizes = (/   -1,   nlat, nlon/)       ;维大小，-1表示未知
    
15.    dimUnlim = (/   True, False,False/)      ;定义是否不限定维
    
16.    filedimdef(fout, dimNames, dimSizes, dimUnlim)
    
17. 
    
18.    ;定义变量名，变量属性和维度
    
19.    filevardef(fout, "time", typeof(time), getVarDimNames(time))
    
20.    filevardef(fout, "lat", typeof(lat), getVarDimNames(lat))
    
21.    filevardef(fout, "lon", typeof(lon), getVarDimNames(lon))
    
22.    filevardef(fout, "plcz", "float", (/ "time","lat","lon"/))
    
23.    ;给每个变量复制属性
    
24.    filevarattdef(fout, "time", time)
    
25.    filevarattdef(fout, "lat", lat)
    
26.    filevarattdef(fout, "lon", lon)
    
27.    filevarattdef(fout, "plcz", plcz)
    
28.    ;将数据写入文件，仅需写入数据，变量属性已经在前面预先写入文件
    
29.    fout->time = (/ time /)          ;用(/.../) 去除变量的元数据
    
30.    fout->lat = (/ lat /)
    
31.    fout->lon = (/ lon /)
    
32.    fout->plcz = (/ plcz /)
    

复制代码