CDO

The scripting language bindings for Ruby have been enlarged to Python! Please have a look at cdo.{rb,py} for installation instructions and some examples.

For reading all unit tests, please register yourself and have a look at source:trunk/cdo/contrib/python/test/test_cdo.py for Python and source:trunk/cdo/contrib/ruby/test/test_cdo.rb for Ruby.

This is the first release of cdo.rb, a Ruby module to call CDO.

Installation¶

cdo.rb is released as a ruby gem, Ruby's own package format. Install cdo.rb with

gem install cdo

Documentation¶

Please have a look at the automatic documentation from rubygem: http://rubydoc.info/gems/cdo/frames

Examples¶

CDO's operators are implemented as module methods:

Cdo.info(:in => "input.nc")

Input and output stream arguments are handled special with hash keys :in and :out. Same is true for the other options of CDO like "-f" or "-P":

Cdo.topo(:out => "topography.grb")
Cdo.copy(:in => "topography.grb",:out => "topography.nc4",:options => "-f nc4")

Operator options can be set as standard method arguments:

Cdo.stdatm(0,10,20,50,100,200,500,1000,:out => "standard_atm.nc",:options => "-f nc")

More example can be found in the unit tests: source:trunk/cdo/contrib/ruby/test/test_cdo.rb.

• Variable input parameter for ECA operators: eca_cdd, eca_cwd, eca_rr1, eca_sdii

• deflate compression with netCDF4 doesn't work (option: -z zip)
• sellonlatbox: correct lon bounds if necessary
• ifthen, ifthenelse: uses only the first time step of the first input file
• module Monarith (monadd, monsub, monmul, mondiv): wrong result for 3D variables

• replace: added support to replace single levels

• remapeta: Changed minimum pressure level for condensation from 1000Pa to 0Pa.
  Use the environment variable REMAPETA_PTOP to set the minimum pressure level for condensation.
  Above this level the humidity is set to the constant 1.E-6.

• invertlat: bug fix for CURVILINEAR grids
• ymon<stat>: preserve time axis attributes (type and calendar)
• import_binary: added support for OPTION ZREV
• expr/exprf: wrong result for expression 'constant-field' and 'constant/field' (e.g. 1-field)
  This bug was introduced in CDO version 1.5.1.
• eof, eoftime, eofspatial, eof3d - Empirical Orthogonal Functions:
  There was a bug in the calculation of the Frobenius norm, which has only been triggered in some cases
  when using a low precision. The normalization has been changed thus that the eigenvectors are not weighted
  and their absolute is 1. The default settings for convergence have been changed to be more conservative:
  CDO_SVD_MODE=jacobi
  MAX_JACOBI_ITER=12
  FNORM_PRECISION=1.e-12

• Added support for netCDF4 classic format; option -f nc4c
• import_binary: Added support for 1 and 2 byte integer

• intlevel3d: vertical interpolation to/from 3d vertical coordinates
• ensrkhistspace: Ranked Histogram averaged over time
• ensrkhisttime: Ranked Histogram averaged over space
• ensroc: Ensemble Receiver Operating characteristics
• enscrps: Ensemble CRPS and decomposition
• ensbrs: Ensemble Brier score

• Exprf: wrong result for missing values != (double) -9.e33
• detrend: added support for time bounds
• Filter: added support for time bounds
• eofspatial: integer overflow; wrong result for grid size > 46340
• eca_*: use the input calendar for the output streams

• netcdf-4.1.2
• hdf5-1.8.7
• grib_api-1.9.9

• Compile grib_api + jasper with PIC support when using make libs4cdo.

Download it here

• GRIB2 support via ECMWF GRIB_API
• Added support for netCDF level bounds
• Added option -O to overwrite existing output file (only for ens<STAT>, merge, mergetime)

• selparam: Select parameters by identifier
• delparam: Delete parameters by identifier
• splitparam: Split parameter identifiers
• chparam: Change parameter identifier

• expr: added functions abs(), int(), nint(), sqr()
• sinfo: changed output of table and code number to parameter identifier
• info: changed output of code number to parameter identifier
• diff: changed output of code number to parameter identifier

• sellonlatbox: does not work as expected when selecting all longitutes
• sellonlatbox: initialization missing for unstructured grids
• ml2pl and ml2hl: wrong result if input file contains full and half level data
• trend and detrend: integer overflow; wrong result for ntimesteps > 46340
• spcut: only correct results with continuous wave numbers starting at 1
• remaplib: fixed data race in calculation of bin_addr (OpenMP)

• hdf5 1.8.6
• grib_api 1.9.8
• proj4 4.7.0

• Add grib_api to default build on unix (with make libs4cdo) for CDO's upcomming GRIB2 support
• Remove netcdf-4.0.1

Downloads

• improved support for netCDF output from WRF model (import time and grid variables)

• ydayadd, ydaysub, ydaymul, ydaydiv: Multi-year daily arithmetic
• eca_pd: Precipitation days index per time period
• dv2ps: Divergence and vorticity to velocity potential and stream function

• import_cmsaf: added time information also for time constant fields

• eof: fix memory access violation
• eofcoeff: fix memory access violation
• fldmean: gives wrong result for grid units [radian]
• Yseasstat: bug fix for datasets with time constant fields
• sellevel: fix problem with hybrid model levels and netCDF output
• sellonlatbox: fix rounding error of the last lon index
• Settime: bug fix for time independent variables in combination with other operators (pipes)

• Using libtool for linking (rpath)
• Changed predefined gaussian grid names from t<RES>grid to n<N>
  Use n80 instead of t106grid to define a Gaussian N80 grid
• Changed percentile parameter type from integer to float

• bandpass: Bandpass filtering
• lowpass: Lowpass filtering
• highpass: Highpass filtering

• eca_gsl: adjust implementation to fit definition by ECA
• expr, exprf: added missing values support

• sellevel: copy zaxis meta data name and units
• seldate: open output file only when time steps found