CDO学习1 CDO简介[通俗易懂]

CDO学习1 CDO简介[通俗易懂]参考自如下网站http://www.ceda.ac.uk/static/media/uploads/ncas-reading-2015/cdo.pdf介绍一个有几百种操作符的单独命令CDO受N

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

参考自如下网站 

http://www.ceda.ac.uk/static/media/uploads/ncas-reading-2015/cdo.pdf

介绍

  • 一个有几百种操作符的单独命令
  • CDO受NCO启发,通过命令行,提供了一系列的气候数据相关的操作
  • 主要是被设计用在操作netCDF3/4,GRIB1/2
  • 众多的函数可以被用在任何NetCDF/格点化数据
  • 非常适合与特定任务
  • 有效的调动内存

CDO:一个命令行工具

$ cdo <operator> [options] <files>

CDO操作的类别

  1. 文件信息和文件操作符
  2. 选择和比较
  3. 元数据的修改
  4. 算术操作符
  5. 统计分析
  6. 回归和内插
  7. 矢量和谱转换
  8. 格式化I/O
  9. 气候指数

众多操作符

CDO将其的功能分解成单独的操作符

当前(2015年)有超过650个操作符。如

showstdname Show standard names

sellonlatbox Select a longitude/latitude box

setmissval Set a new missing value

monadd Add monthly time series

zonstd Zonal standard deviation

eca_hd Heating degree days per time period

得到操作符的参考

如需参考,只需输入:cdo -h <operator>

 操作符:1.文件信息

 $ cdo infov ifile 

这是一个带有一个2D变量的数据集的示例结果,包括3个时间步长:

-1 : Date Time Varname Level Size Miss : Minimum Mean Maximum
1 : 1987?01?31 12:00:00 SST 0 2048 1361 : 232.77 266.65 305.31
2 : 1987?02?28 12:00:00 SST 0 2048 1361 : 233.64 267.11 307.15
3 : 1987?03?31 12:00:00 SST 0 2048 1361 : 225.31 267.52 307.67

 

$ cdo showtimestamp ggas2014121200_00-18.nc
2014-12-12T00:00:00 2014-12-12T06:00:00 2014-12-12T12:00:00 2014-12-12T18:00:00

 

$ cdo pardes ggas2014121200_00-18.nc
-1 CI Sea-ice cover [(0 - 1)]
-2 SSTK Sea surface temperature [K]
-3 MSL Mean sea-level pressure [Pa]
-4 TCC Total cloud cover [(0 - 1)]
-5 U10 10 metre U wind component [m s**-1]
-6 V10 10 metre V wind component [m s**-1]
-7 SKT Skin temperature [K]

 

很容易看到2个文件之间数值差异的概况:

$ cdo diff tas_rcp45_2055_ann_95p_change.nc tas_rcp45_2055_ann_05p_change.nc
Date Time Param Level Size Miss : S Z Max_Absdiff Max_Reldiff
1 : 2065-12-30 12:00:00 -1 0 64800 0 : T F 16.693 0.99997
1 of 1 records differ

 

操作符的所有文件列表包括

  • 赋值和连接
  • 合并场和时间
  • 按变量/层次/格点/时间分割场

一些例子

To copy a file and convert the output to NetCDF:

 $ cdo −f nc copy ifile ofile.nc 

To merge all files along the time axis:
 $ cdo mergetime ifile1 ifile2 ifile3 ofile 

可以用一些不同的方法选择数据包括:

  • 变量码,名称或者属性
  • 层次
  • 时间
  • 空间区域(纬度/经度)

 

操作符:2.选择

用id选择变量:“SSTK”和”CI”:

 $ cdo selname,SSTK,CI infile.nc outfile.nc 

To select a lat/lon bounding box:
 $ cdo sellonlatbox,120,-90,20,-20 infile.nc outfile.nc 

To select a date/time range:
 $ cdo seldate,201412-12T12:00:00, 201501-31T18:00:00 infile.nc outfile.nc 

 

操作符:3.修改

设置时间轴到198701-16 12:00, 时间增量一个月增量 使用

 $ cdo settaxis,19870116,12:00,1mon ifile ofile 

将一个2维场纬度从N到S反转成S到N

 $ cdo invertlat ifile ofile 

 

操作符:4.算术

计算所有场元素的平方根:

 $ cdo sqrt ifile ofile 

将所有输入场加入一个常数-273.15:

 $ cdo -addc,-273.15 ifile ofile 

 

操作符:5.统计

计算所有输入场的纬向平均:

 $ cdo zonmean ifile ofile 

假设输入数据集具有数年的月度平均值。 为了从月度计算季节性均值,必须跳过前两个月:Assume an input dataset has monthly means over several years. In order to compute seasonal means from monthly means the first two months must be skipped:

 $ cdo timselmean,3,2 ifile ofile 

多年逐日滑动百分值:

  • 为了计算一个百分值,必须知道最小和最大边界。

 $ cdo ydrunpctdl,p,nts infile minfile maxfile outfile 

该运算符写滑动的百分值对于一年中的每一天,从infile写到outfile。最小和最大边界用minfile和maxfile提供。

计算一个对年逐日的滑动百分率值,需要执行多步。

$ cdo ydrunmin,5 ifile minfile
$ cdo ydrunmax,5 ifile maxfile
$ cdo ydrunpctl,90,5 ifile minfile maxfile ofile

 

操作符:6 内插

为了将所有场线性插值到一个T42高斯格点。

 $ cdo remapbil,t42grid ifile ofile 

为了将混合模式层数据插值到气压层925,850,500和200 hPa

 $ cdo ml2pl,92500,85000,50000,20000 ifile ofile 

 

操作符:9 气候指数

得到一个逐日降水总量时间序列找那个连续干日的最大数目。

 $ cdo eca_cdd rrfile ofile 

其中rrfile是逐日降水总量RR的时间序列,然后统计了连续天数的最大值,当RR小于1 mm时。

为了得到热带夜晚一个时间序列逐日最小温度

 $ cdo eca_tr tnfile ofile 

其中tnfile是逐日最低温度TN的时间序列,随后统计的是当TN>T时候的天数。T是一个可选的参数,其缺省值T=20℃。

 

除此之外CDO还能做什么? CDO有许多通用的命令行开关,可以与(大多数)运算符一起使用。 这是一个选择:

-a 生成绝对时间轴Generate an absolute time axis
-f <format> 指定格式Specify format (“grb”, “nc”, “nc4” etc)
-m <val> 设置默认缺测值Set the default missing value
-Q 对netCDF变量名排序Sort netCDF variable names
-r 生成相对时间轴Generate a relative time axis
-s 静默模式Silent mode
-V Version of CDO
-v Verbose – print extra details.
-z zip Deflate compression of netCDF4 vars.

 

组合运算符

所有的运算符都有一个固定的输入流和一个输出流,可以直接将结果管道给其它的操作符。这个操作符必须以“-”开始,为了和其它相结合。可以通过一下,来促进性能:

  • 减少不必要的磁盘I/O
  • 并行执行

例如,我们可以把

$ cdo timavg ifile1 tmp1
$ cdo dayavg ifile2 tmp2
$ cdo sub tmp2 tmp1 ofile
$ rm tmp1 tmp2

替换成

$ cdo sub -dayavg ifile2 -timavg ifile1 ofile

 

组合运算符:例子2

以下三个命令是计算一个多年逐日滑动百分率,需要多步:

$ cdo ydrunmin,5 ifile minfile
$ cdo ydrunmax,5 ifile maxfile
$ cdo ydrunpctl,90,5 ifile minfile maxfile ofile

可以通过一个简单的命令替换

$ cdo ydrunpctl,90,5 ifile -ydrunmin ifile -ydrunmax ifile ofile

 

更多信息 

CDO User Guide: https://code.mpimet.mpg.de/projects/cdo/embedded/cdo.pdf

CDO Home page:  https://code.zmaw.de/projects/cdo

CDO Tutorial:  https://code.zmaw.de/projects/cdo/wiki/Tutorial

CDO Documentation:  https://code.zmaw.de/projects/cdo/wiki#Documentation

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/166771.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • eclipse导入maven工程出现的问题「建议收藏」

    eclipse导入maven工程出现的问题「建议收藏」最近自己找了一个ssm框架想学习一下,但是用eclipse导入maven工程后出现了以下问题:error1:DescriptionResourcePathLocationTypeProjectconfigurationisnotup-to-datewithpom.xml.RunMaven4MyEclipse->UpdateProject

  • Nodejs安装教程

    Nodejs安装教程nodejs和npm安装详细教程

  • SpringCloud-Alibaba-Sentinel-服务降级-热点限流-服务熔断

    SpringCloud-Alibaba-Sentinel-服务降级-热点限流-服务熔断前言:除了流量控制以外,对调用链路中不稳定的资源进行熔断降级也是保障高可用的重要措施之一。一个服务常常会调用别的模块,可能是另外的一个远程服务、数据库,或者第三方API等。例如,支付的时候,可能需要远程调用银联提供的API;查询某个商品的价格,可能需要进行数据库查询。然而,这个被依赖服务的稳定性是不能保证的。如果依赖的服务出现了不稳定的情况,请求的响应时间变长,那么调用服务的方法的响应时间也会变长,线程会产生堆积,最终可能耗尽业务自身的线程池,服务本身也变得不可用熔断策略Sentinel提供

  • 金蝶K3-航天信息税控发票开票软件接口程序

    金蝶K3-航天信息税控发票开票软件接口程序##金蝶K3-航天信息税控发票开票软件接口程序用于从航天信息[上海爱信诺航天信息有限公司]税控发票开票软件(金税盘版)文本接口导出TXT文件中提取数据。###开发环境SharpDevelop5.0###调用库NPOI(用于EXCEL文件生成)###背景公司之前使用一个金蝶开发的…

  • mysql smalldatetime_datetime与smalldatetime之间的区别

    mysql smalldatetime_datetime与smalldatetime之间的区别1、一直以为smalldatetime和datetime的差别只是在于时间范围:smalldatetime的有效时间范围1900/1/1~2079/6/6datetime的有效时间范围1753/1/1~9999/12/31所以我判断如果该值不用到太远的日期范围,就会使用smalldatetime。2、但我忽略了更关键的差别,那就是smalldatetime只精准到分,而datetime则可精准到3…

  • Anaconda 的安装教程(图文)「建议收藏」

    Anaconda的介绍Anaconda指的是一个开源的Python发行版本,其包含了conda、Python等180多个科学包及其依赖项。因为包含了大量的科学包,Anaconda的下载文件比较大。这么说可能有点抽象,大家可以直接把Anaconda理解为一个python的傻瓜捆绑包就行了Anaconda下载下载地址:https://www.anaconda.com/download/…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号