CDO

New features:

• Proj 4 to 6 API Migration
• smooth/smooth9: Added support for gridtype PROJECTION [Feature #9202]
• Expr: Added function rand() and isMissing()
• Remap: Added support for Gaussian reduced grids
• trend, detrend: Added parameter equal=false for unequal timesteps
• Option --no_remap_weights: Switch off generation of remap weights

New operators:

• deltat: Difference between timesteps

Fixed bugs:

• ensavg: Wrong result if data contains missing values (same result as ensmean)
• Ydrunstat: Fix seg. fault

Currently I work on a redesign of the Python and Ruby interface to CDO towards a 2.0 release. This will include API changes that are not backward compatible to the current 1.5.x version.

I wrote a list of open points here and I'd like to invite everyone to join the discussion and come up with ideas on how to write down things in the future. Please create tickets at https://github.com/Try2Code/cdo-bindings/issues so that ideas and wishes can be collect in a transparent way.

• added option --worker <num>: Number of worker to decode/decompress GRIB records
• added option --pedantic: Warnings count as errors
• Yhourstat: added time bounds support
• expr: added support for ctimestep() in ternary conditional

• yearmaxidx: Yearly maximum indices
• yearminidx: Yearly minimum indices

• for: renamed to seq

• Build failed with GCC 9 (OpenMP data sharing) [Bug #9038]
• compile error: EXIT_FAILURE not declared in cdoDebugOutput.h [Bug #8899]
• eca_gsl: the 2nd input file was not closed [Bug #9033]
• ensrkhisttime/ensrkhistspace: don't work
• detrend: seg. fault if time series containts time constant fields
• inttime, intntime: handling of missing values is incorrect
• select: combination of some parameter (var, grid, zaxis) doesn't work
• expr:zonSTAT: wrong result
• expr::vertmean: fix wrong warning message about layer bounds
• mergetime: SKIP_SAME_TIME doesn't work in release 1.9.6

• Added support for polar stereographic projection
• Download ICON grids if necessary (http://icon-downloads.mpimet.mpg.de)
• Added global option --eccodes: Use ecCodes to decode/encode GRIB1 messages
• Operator chname: Added support to change coordinate names [Feature #8746]
• Operator diff: set exit status to 1 if inputs differ
• Renamed remapcon/gencon to remapscon/genscon
• Replaced remapcon/gencon by remapycon/genycon

• vars<stat> - Statistical values over all variables
• minc - Minimum of a field and a constant
• maxc - Maximum of a field and a constant

• Gradsdes.test fails [Bug #8614]
• Option --reduce_dim gives wrong result on time dimension [Bug #8615]
• Module Selbox: added support for grid cell area
• Operator setgridtype,regular: set nx=4*N+16 for octahedral reduced Gaussian grids
• Operator distgrid: seg. fault if last segment is larger than first segment
• Operator sellonlatbox: abort if grid coordinates missing
• Operator masklonlatbox: wrong result if lon1 > first lon || lon2 < last lon (bug introduce in 1.9.4) [Bug #8695]
• Operator maskindexbox: wrong result if idx1 > 1 || idx2 < nlon (bug introduce in 1.9.4) [Bug #8695]
• Absolute time axis (-a) returns wrong units in operator chain for NetCDF [Bug #8777]
• Relative time axis (-r) returns wrong first timestep in operator chain for NetCDF
• Wrong result with fldmean on zonal mean data (bug introduce in 1.9.5) [Bug #8834]

• Changed type of date from 32 to 64-bit integer to support years > 214748
• remapycon: optimized by changing cell search method
• expr: added support for zon<STAT> functions
• expr: added function sellevelrange() and sellevidxrange()
• gridfile: added extension ":N" to select grid number N from data file

• dhour<stat> - Multi-day hourly statistics

• seldate: segmentation fault (bug introduce in 1.9.4) [Bug #8499]
• select: wrong timestamp when combining select with selyear (bug introduce in 1.8.1) [Bug #8576]
• gradsdes: bug fix for rotated lon/lat grids
• silent option produces newlines [Bug #8538]
• remapnn/remapdis: wrong result with regular 2D source grids if nlat > nlon [Bug #8498]

New features:

• Large data support: changed type of gridsize from 32 to 64-bit integer
• remapbil, remapbic, remapnn, remapdis: optimized by changing point search method
• Fldstat, Vertstat: added option weight=false to disable weighting

Fixed bugs:

• option -r doesn't work [Bug #8334]
• enspctl: changed parameter type from int to float [Bug #8386]
• segfault with chained operators on timeseries data [Bug #8230]
• setattribute: added support for \n in text attributes
• expr: removed character [LlDd] from definition of float constants

• expr: added time coordinate function cdate(), ctime(), cdeltat(), ctimestep() ...

• not - logical NOT (1, if x equal 0; else 0)

• uvDestag: target grid undefined in output
• runpctl: fails since release 1.8.0
• read of reduced Gaussian grid description file failed [Bug #8146]
• read error on grid description file [Bug #8099]

• sign of grid size increment changes [Bug #7974]
• compilation fails on OpenBSD [Bug #7961]
• expr: nesting of ternary operator lost in cdo-1.9.1 [Bug #7992]
• rotuvb changed behavior in different versions [Bug #8084]
• select with start=end range aborts with 'Invalid character' [Bug #7976]

• Added support for NC_FORMAT_CDF5
• Extend option --reduce_dim to all dimension for all operators

• tee - Duplicate a data stream

• eof, eof3d: set default value of environment variable CDO_WEIGHT_MODE to off
• sinfo: Added time type
• ap2pl: added support for input data on half levels

• selindexbox: breaks uvRelativeToGrid flag [Bug #7901]
• expr: AND fall through OR
• --cmor option doesn't work for lon/lat bounds (introduced in 1.9.0)
• eof3d: weight array was allocated for only one level

• Code changed from ANSI C99 to ISO C++11
• Added configure option for ecCodes --with-eccodes=<yes|no|directory>
• Added range operator to all statistic modules (e.g. yearrange, zonrange)

• expr: improve ternary operator, no brackets needed anymore.
• expr: added support for clev in ternary operator.
• remapcon/remapycon produces wrong results for some grid combinations (introduced in 1.8.0) [Bug #7821]
• mergetime: wrong time information if first input file does not contain the first time step (bug introduced in 1.8.1) [Bug #7760]
• percentile: fix wrong result with method numpy (linear interpolation) and nist [Bug #7798]