最近接觸單柱模式，需要處理輸出NC文件，即NetCDF文件(Network Common Data Format)，簡(jiǎn)單了解了NCL和NCO工具。一個(gè)NetCDF數(shù)據(jù)集包含維(dimensions)、變量(variables)和屬性(attributes)三種描述類型，每種類型都會(huì)被分配一個(gè)名字和ID，這些類型共同描述了一個(gè)數(shù)據(jù)集，NetCDF庫(kù)可以同時(shí)訪問(wèn)多個(gè)數(shù)據(jù)集，用ID來(lái)識(shí)別不同數(shù)據(jù)集。變量存儲(chǔ)實(shí)際數(shù)據(jù)，維給出了變量維度信息，屬性則給出了變量或數(shù)據(jù)集本身的輔助信息屬性，又可以分為適用于整個(gè)文件的全局屬性和適用于特定變量的局部屬性，全局屬性則描述了數(shù)據(jù)集的基本屬性以及數(shù)據(jù)集的來(lái)源。
流行的處理軟件有很多，對(duì)比如下：

Paste_Image.png

Paste_Image.png

NCL

NCL全稱為NCAR Command Language，是由美國(guó)大氣研究中心（NCAR）開發(fā)的 專門用于科學(xué)資料分析與可視化的一門解釋型語(yǔ)言。
NCL支持的數(shù)據(jù)類型包括NetCDF 3/4，GRIB 1/2, HDF 4/5，HDF-EOS 2/5， shapefile，ASCII，binary。內(nèi)建了大量的分析函數(shù)。可以很方便的產(chǎn)生 高質(zhì)量的圖形并且有大量的圖形源（resource）可被用于自定義圖形。大量 示例腳本和相應(yīng)的圖形可在官網(wǎng)（ncl.ucar.edu）學(xué)習(xí)。

基本命令

參考資料

官網(wǎng)

Paste_Image.png

NCO工具

netCDF Operators (NCO) 工具用來(lái)操作和分析netCDF自描述數(shù)據(jù)存儲(chǔ)格式。
nco 使用頻率較高的兩中處理方式就是average與concatenation
添加(Appending)是針對(duì)變量的。
file1 = T,U,V
file2 = P,C
file3 = T,U,V,P,C
連接(Concatenation)是針對(duì)（如：時(shí)間）維度的。
file1 = T(0-12,:,:)
file2 = T(13-24,:,:)
file3 = T(0-24,:,:)

基礎(chǔ)命令

ncks filename | more #查看文件
ncks -v variable filename
ncrcat #將已有record集合為一個(gè)文件
ncecat#將ensemble集合為一個(gè)文件
ncra #將record取平均
ncwa#將record的某些部分取平均， -a 是去掉某些變量
ncdiff #取不同
ncea #變量取平均
ncrename #變量重命名
cdo #變量插值
ncpdq -U #變量類型轉(zhuǎn)換

參考資料

官網(wǎng)網(wǎng)站http://nco.sourceforge.net/nco.html
http://bosomfriend.blog.edu.cn/home.php?mod=space&uid=1555165&do=blog&id=75336
http://hxfcalf.blog.163.com/blog/static/21575548201051191420134

實(shí)踐

1. 在我的實(shí)驗(yàn)中只用到了幾個(gè)命令
   ncrename ##rename dimensions, variables, and attributes in netCDF file.
   ncdump -h XXX.nc #查看文件結(jié)構(gòu)
   ncks XXX.nc | more##查看數(shù)據(jù)內(nèi)容
   ncks -H -v u -d lev,15 -d time,56,1136 ARM95_4scam.nc U850.nc ##查看搜索lev數(shù)組第15元素，time數(shù)組在56，1136標(biāo)號(hào)內(nèi)，變量為u的數(shù)據(jù)并輸出到U850.nc文件中
   ncwa -a time,lev U850.nc -O U850.nc ##平均掉time和氣壓水平

2. cnl文件
   <pre>
   load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
   load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/ut_string.ncl"
   begin
   ;var_names = (/"LWCF", "SWCF", "PRECT", "Q850", "T850"/)
   var_names = (/"U200","U850","V200","V850","PRECT","Q850","FLUT","FSNTOA","T850","T200","CLDHGH","CLDLOW","CLDMED","CLDTOT","PS"/)
   obs_path = "/home/huangxin/CSM-tunner/SCAM/cesm1_2_2/SCMtest_11111/run/metrics/climo/obs"
   model_base = "/home/huangxin/CSM-tunner/SCAM/cesm1_2_2/SCMtest_11111/run/metrics/climo/model/tune/"
   climo_path = "/home/huangxin/CSM-tunner/SCAM/cesm1_2_2/SCMtest_11111/run/metrics/climo/"

   num_vars = dimsizes(var_names)
   var_obs = new((/num_vars/), float)
   var_obs@_FillValue = 1.0e33
   var_obs = var_obs@_FillValue
   var_model = new((/num_vars/), float)
   var_model@_FillValue = 1.0e33
   var_model = var_model@_FillValue
   var_rms = new((/num_vars/), float)
   var_rms@_FillValue = 1.0e33
   var_rms = var_rms@_FillValue
   do k = 0, num_vars - 1
   obs_file = obs_path + "/" + var_names(k) +".nc"
   model_file = model_base + "tmp.nc"
   fp_obs = addfile(obs_file, "r")
   ;print(fp_obs->$var_names(k)$)
   var_obs(k) = fp_obs->$var_names(k)$(0,0)
   fp_model = addfile(model_file, "r")
   var_model(k) = fp_model ->$var_names(k)$(0,0,0)
   end do
   var_rms = (var_obs - var_model)^2
   </pre>

3. 未解決：利用bash的while循環(huán)按行讀文件(如下)，每讀一行處理一次ncl，但會(huì)遇到語(yǔ)法問(wèn)題:syntax: line -1
   <pre>
   while read variable
   do
   ...
   ncl XXX.ncl
   ...
   done
   </pre>
   目前用csh來(lái)處理，或bash的其他效率低的讀文件方式。