Cdo command for NetCDF

Cdo command for NetCDF weather file

CDO DOES NOT HAVE a graphic environment! It only processes data. The commands run directly on the Linux terminal.

Note: For NetCDF data support, the installation of the NetCDF libraries are required;

Syntax: cdo <operator> input.nc output.nc or
cdo <operator> input.nc.

首先:export PATH=/public/software/anaconda2/bin:$PATH

1. File Information:

1 infon: information about the dataset listed by the variable name.
2 sinfon: information on the dataset in a reduced form.
3 griddes: shows information about the ?le domain.
4 zaxisdes: shows information about the vertical coordinate.
5 pardes: description of the variable or variables of the ?le.

cdo infon standard.nc | head
cdo infov ifile 
cdo showtimestamp ~
cdo sellonlatbox  ~

2. File Manipulation:

1 copy: duplicate the ?le and/or convert to another format.
2 cat: increases the number of ?le times.
3 merge: joins ?les with di?erent variables.
4 split: divide the ?le (vertical level or time).

copy: duplicates the ?le and/or converts from one format to another.
Note:The ?les must have the same dimensions (latitude, longitude, vertical levels).

Example1: Converts from NetCDF to GRIB:
cdo -f grb copy pr.GPCP.1996.2005.nc pr.GPCP.1996.2005.grb

Example2: Converts from GRIB to NetCDF:
cdo -f nc copy pr.GPCP.1996.2005.grb pr.GPCP.1996.2005.nc

Example of concatenation with the copy operator:
cdo copy pr.GPCP.1996.1999.nc pr.GPCP.2000.2002.nc pr.1996.2002.nc

cat: increases the number of ?le times.
cdo cat pr.GPCP.1996.1999.nc pr.GPCP.2000.2002.nc pr.1996.2002.nc

merge: joins ?les with di?erent variables.
cdo merge u.NCEP.R2.1996.2005.nc v.NCEP.R2.1996.2005.nc uv.nc
Note:Files do not necessarily have to have the same dimensions (latitude, longitude, vertical levels).

split: divide or separate the ?le.
Operator: splitlevel, splithour, splitday, splitmon, splitseas and splityear
Example: Separate ?le in months:
cdo splitmon pr.GPCP.1996.2005.nc month.

Where: month. represents a pre?x or any name de?ned by user and pr.GPCP.1996.2005.nc is the input ?le.There will be 12 ?les created: month.01.nc, month.02.nc, ... ,month.12.nc.

3. Selection Operators:

1 sellevel: selects the vertical level.
2 selmon: selects the month of the ?le.
3 selseason: selects the season of the year.
4 selyear: selects the year ?le.
5 sellonlatbox: cut the data in the domain of interest

cdo selyear,2000/2002 tar.NCEP.R2.1996.2005.nc tar.2000.2002.nc
cdo selyear,1996/1998,2000 tar.NCEP.R2.1996.2005.nc tar.1996a1998.2000.nc
cdo seldate,2014-12-12T12:00:00, 2015-01-31T18:00:00 infile.nc outfile.nc 
sellonlatbox: cut the data according to the area of interest.
cdo sellonlatbox,-90,-30,-60,10 pr.GPCP.1996.1999.nc pr.AS.nc

The convention will always be: longitude west, longitude east, latitude south and north latitude

Longitude values can be reported in either the 0 to 360 format and in the -180 to 180 format.

How to convert the global data to -180 to 180? Only works with global data.
cdo sellonlatbox,-180,180,-90,90 pr.GPCP.1996.1999.nc pr.global.nc

setcalendar: sets the calendar type of the ?le.
Possibilities: standard, proleptic_gregorian, 360_day, 365_day e366_day.

ncdump -h stan.nc

cdo setcalendar,365_day stan.nc st.nc

setmissval: sets a new missing value (UNDEF).
cdo setmissval,-999 tar.GFDL-ESM2M.1996.2005.nc tar.GFDL-ESM2M.nc

4. Operations with constants:

1 addc: sum to a constant.
2 subc: subtracts from a constant.
3 mulc: multiplies by a constant.
4 divc: divided by a constant.

乘法
Converting the precipitation of [kg/(m2 s)] to (mm/day).
cdo mulc,86400 pr.CANESM2.1996.2005.nc pr.mensal.nc
減法
Converting Kelvin temperature (K) to degrees Celsius (oC)
cdo subc,273.15 tar.NCEP.R2.1996.2005.nc tc.nc

5. Operations using two sets of data (uni or two-dimensional).

1 add: adds two ?elds.
2 sub: subtract two ?elds.
3 mul: multiplies two ?elds.
4 div: divide dois campos.

The ?les must have the same dimensions (latitude, longitude, vertical levels).

cdo sub pr.CANESM2.2005.nc pr.GPCP.2005.nc dif.nc
Annual statistical calculation. Converts the data to annual time scale
independent of the temporal resolution.
For example, if the data is hourly, daily or monthly, it will be converted
for an annual archive.
1 yearsum: annual sum.
2 yearmean: annual average.
3 yearavg: annual average.

Seasonal statistical calculation::
1 seassum: seasonal sum.
2 seasmean: seasonal average.
3 seasavg: seasonal average.

Statistical calculation - converts to monthly scale..
1 ymonsum: monthly sum.
2 ymonmean: monthly average.
3 ymonavg: monthly average.

Operators to calculate correlation.
?dcor ? correlates all grid points for each time of the two
variables. The result is a time series of correlation between them.
timcor ? performs the correlation at each grid point for all times
of the two variables. The result is a ?le spatial correlation with
only one time.

Example1: Correlation between precipitation and temperature. The
data are monthly referring to the year of 2004 (12 months or times).
cdo ?dcor prec.mensal.nc temp.mensal.nc output.nc

Examples of other applications:

1 Calculate the anomaly.
2 Create masking.
3 Filling missing data.
4 Change NetCDF ?le values.
5 Text ?le conversion to NetCDF.
6 Conversion of a radiosonde to the NetCDF format.
最后編輯于
?著作權歸作者所有,轉載或內容合作請聯系作者
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發布,文章內容僅代表作者本人觀點,簡書系信息發布平臺,僅提供信息存儲服務。
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市,隨后出現的幾起案子,更是在濱河造成了極大的恐慌,老刑警劉巖,帶你破解...
    沈念sama閱讀 228,505評論 6 533
  • 序言:濱河連續發生了三起死亡事件,死亡現場離奇詭異,居然都是意外死亡,警方通過查閱死者的電腦和手機,發現死者居然都...
    沈念sama閱讀 98,556評論 3 418
  • 文/潘曉璐 我一進店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人,你說我怎么就攤上這事。” “怎么了?”我有些...
    開封第一講書人閱讀 176,463評論 0 376
  • 文/不壞的土叔 我叫張陵,是天一觀的道長。 經常有香客問我,道長,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 63,009評論 1 312
  • 正文 為了忘掉前任,我火速辦了婚禮,結果婚禮上,老公的妹妹穿的比我還像新娘。我一直安慰自己,他們只是感情好,可當我...
    茶點故事閱讀 71,778評論 6 410
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著,像睡著了一般。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發上,一...
    開封第一講書人閱讀 55,218評論 1 324
  • 那天,我揣著相機與錄音,去河邊找鬼。 笑死,一個胖子當著我的面吹牛,可吹牛的內容都是我干的。 我是一名探鬼主播,決...
    沈念sama閱讀 43,281評論 3 441
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了?” 一聲冷哼從身側響起,我...
    開封第一講書人閱讀 42,436評論 0 288
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后,有當地人在樹林里發現了一具尸體,經...
    沈念sama閱讀 48,969評論 1 335
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...
    茶點故事閱讀 40,795評論 3 354
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發現自己被綠了。 大學時的朋友給我發了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 42,993評論 1 369
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖,靈堂內的尸體忽然破棺而出,到底是詐尸還是另有隱情,我是刑警寧澤,帶...
    沈念sama閱讀 38,537評論 5 359
  • 正文 年R本政府宣布,位于F島的核電站,受9級特大地震影響,放射性物質發生泄漏。R本人自食惡果不足惜,卻給世界環境...
    茶點故事閱讀 44,229評論 3 347
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧,春花似錦、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 34,659評論 0 26
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至,卻和暖如春,著一層夾襖步出監牢的瞬間,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 35,917評論 1 286
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人。 一個月前我還...
    沈念sama閱讀 51,687評論 3 392
  • 正文 我出身青樓,卻偏偏與公主長得像,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 47,990評論 2 374

推薦閱讀更多精彩內容

  • rljs by sennchi Timeline of History Part One The Cognitiv...
    sennchi閱讀 7,392評論 0 10
  • 我們小時候看世界,只有三種維度,更小的時候,連三種都沒有。它們是,長、寬、高。所以,我們叫他們大人,我們是小人。 ...
    周某某2019閱讀 174評論 0 0
  • 小白 男 我的“閨蜜” 小白是我給他取的專屬“愛稱” 認識小白 確切的說 是和他熟悉 是因為 散伙飯 那時候...
    _半夏閱讀 316評論 0 0
  • 我是林強,公眾號項目03號學員,七天的最后一天了,說說我的感受。 其實這心情很差,想想目前的情況妥妥的一個失敗者,...
    林強_a818閱讀 157評論 0 0
  • 公孫離已經在體驗服一段時間了 S9馬上就要結束了 今天已經1月16號了 傳說中S10將于1月25號開啟 而公孫離也...
    抖S貓叔閱讀 656評論 1 0