intlevelx3d - getting an error where once was none - CDO - Project Management Service

intlevelx3d - getting an error where once was none

Added by Brendan DeTracey almost 4 years ago

Hi all,

I am re-running the intlevelx3d command that worked in a previous cdo version (5 years ago ) but now emits errors with 1.9.8. Help! The input files have not changed.

$ cdo -v -delname,deptho -intlevelx3d,depth3D_CanESM2.nc thetao_Omon_CanESM2_historical_r1i1p1_198101.nc deptho_fx_CanESM2_historical_r0i0p0_filled.nc bottom_thetao_Omon_CanESM2_historical_r1i1p1_198101.nc
 OpenMP:  num_procs=8  max_threads=1
-delname,deptho -intlevelx3d,depth3D_CanESM2.nc thetao_Omon_CanESM2_historical_r1i1p1_198101.nc deptho_fx_CanESM2_historical_r0i0p0_filled.nc bottom_thetao_Omon_CanESM2_historical_r1i1p1_198101.nc 
cdo(1) intlevelx3d: Process started
cdo    delname: name 1 = deptho
Warning (cdfScanVarAttr): NetCDF: Variable not found - lat_bnds
Warning (cdfScanVarAttr): NetCDF: Variable not found - lev_bnds
Warning (cdfScanVarAttr): NetCDF: Variable not found - lon_bnds
cdo(1) intlevelx3d: 40 records input 3d vertical height
cdo(1) intlevelx3d: 1 records target 3d vertical height and gridsize 49152
cdo(1) intlevelx3d: Input vertical coordinate has no missing values.
cdo(1) intlevelx3d: Output vertical coordinate has no missing values.
Warning: Non monotonic zaxis!
cdo(1) intlevelx3d: lev1 0: 5
cdo(1) intlevelx3d: lev1 1: 15
cdo(1) intlevelx3d: lev1 2: 25.1
cdo(1) intlevelx3d: lev1 3: 35.3
cdo(1) intlevelx3d: lev1 4: 45.6
cdo(1) intlevelx3d: lev1 5: 56.1
cdo(1) intlevelx3d: lev1 6: 66.8
cdo(1) intlevelx3d: lev1 7: 77.7
cdo(1) intlevelx3d: lev1 8: 88.9
cdo(1) intlevelx3d: lev1 9: 100.5
cdo(1) intlevelx3d: lev1 10: 112.6
cdo(1) intlevelx3d: lev1 11: 125.5
cdo(1) intlevelx3d: lev1 12: 139.5
cdo(1) intlevelx3d: lev1 13: 154.9
cdo(1) intlevelx3d: lev1 14: 172.2
cdo(1) intlevelx3d: lev1 15: 192.1
cdo(1) intlevelx3d: lev1 16: 215.4
cdo(1) intlevelx3d: lev1 17: 243.1
cdo(1) intlevelx3d: lev1 18: 276.6
cdo(1) intlevelx3d: lev1 19: 317.7
cdo(1) intlevelx3d: lev1 20: 368.6
cdo(1) intlevelx3d: lev1 21: 431.9
cdo(1) intlevelx3d: lev1 22: 510.5
cdo(1) intlevelx3d: lev1 23: 607.7
cdo(1) intlevelx3d: lev1 24: 726.8
cdo(1) intlevelx3d: lev1 25: 870.8
cdo(1) intlevelx3d: lev1 26: 1042.1
cdo(1) intlevelx3d: lev1 27: 1242.1
cdo(1) intlevelx3d: lev1 28: 1471.1
cdo(1) intlevelx3d: lev1 29: 1728.1
cdo(1) intlevelx3d: lev1 30: 2010.9
cdo(1) intlevelx3d: lev1 31: 2316.6
cdo(1) intlevelx3d: lev1 32: 2642
cdo(1) intlevelx3d: lev1 33: 2983.8
cdo(1) intlevelx3d: lev1 34: 3338.9
cdo(1) intlevelx3d: lev1 35: 3704.6
cdo(1) intlevelx3d: lev1 36: 4078.6
cdo(1) intlevelx3d: lev1 37: 4459
cdo(1) intlevelx3d: lev1 38: 4844.3
cdo(1) intlevelx3d: lev1 39: 5233.3
cdo(1) intlevelx3d: lev1 40: 0
cdo(1) intlevelx3d: lev1 41: 0

cdo(1) intlevelx3d (Abort): Level 0 not found!

$ cdo -V
Climate Data Operators version 1.9.8 (https://mpimet.mpg.de/cdo)
System: x86_64-conda_cos6-linux-gnu
CXX Compiler: /home/conda/feedstock_root/build_artifacts/cdo_1582618008272/_build_env/bin/x86_64-conda_cos6-linux-gnu-c++ -fPIC -DPIC -g -O2 -std=c++11 -fopenmp -DACCEPT_USE_OF_DEPRECATED_PROJ_API_H=1  -march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -ffunction-sections -pipe -isystem /home/brendan/miniconda3/envs/cmip6/include -fdebug-prefix-map=/home/conda/feedstock_root/build_artifacts/cdo_1582618008272/work=/usr/local/src/conda/cdo-1.9.8 -fdebug-prefix-map=/home/brendan/miniconda3/envs/cmip6=/usr/local/src/conda-prefix -fopenmp 
CXX version : unknown
C Compiler: /home/conda/feedstock_root/build_artifacts/cdo_1582618008272/_build_env/bin/x86_64-conda_cos6-linux-gnu-cc -DACCEPT_USE_OF_DEPRECATED_PROJ_API_H=1  -march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -ffunction-sections -pipe -isystem /home/brendan/miniconda3/envs/cmip6/include -fdebug-prefix-map=/home/conda/feedstock_root/build_artifacts/cdo_1582618008272/work=/usr/local/src/conda/cdo-1.9.8 -fdebug-prefix-map=/home/brendan/miniconda3/envs/cmip6=/usr/local/src/conda-prefix -fopenmp  
C version : unknown
F77 Compiler: /home/conda/feedstock_root/build_artifacts/cdo_1582618008272/_build_env/bin/x86_64-conda_cos6-linux-gnu-gfortran -fopenmp -march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -ffunction-sections -pipe -isystem /home/brendan/miniconda3/envs/cmip6/include -fdebug-prefix-map=/home/conda/feedstock_root/build_artifacts/cdo_1582618008272/work=/usr/local/src/conda/cdo-1.9.8 -fdebug-prefix-map=/home/brendan/miniconda3/envs/cmip6=/usr/local/src/conda-prefix
F77 version : unknown
Features: 11GB 8threads C++11 Fortran DATA PTHREADS OpenMP HDF5 NC4/HDF5/threadsafe OPeNDAP UDUNITS2 PROJ XML2 CURL FFTW3 SSE3
Libraries: HDF5/1.10.5 proj/6.3.1 xml2/2.9.10 curl/7.68.0
Filetypes: srv ext ieg grb1 grb2 nc1 nc2 nc4 nc4c nc5 
     CDI library version : 1.9.8
 cgribex library version : 1.9.4
 ecCodes library version : 2.17.0
  NetCDF library version : 4.7.3 of Jan 23 2020 20:07:17 $
    hdf5 library version : 1.10.5 threadsafe
    exse library version : 1.4.1
    FILE library version : 1.8.3

Download all files

deptho_fx_CanESM2_historical_r0i0p0_filled.nc (207 KB) deptho_fx_CanESM2_historical_r0i0p0_filled.nc
depth3D_CanESM2.nc (15 MB) depth3D_CanESM2.nc
thetao_Omon_CanESM2_historical_r1i1p1_198101.nc (22.5 MB) thetao_Omon_CanESM2_historical_r1i1p1_198101.nc

Replies (5)

RE: intlevelx3d - getting an error where once was none - Added by Ralf Mueller almost 4 years ago

Hi Brendan!

This error does not happen with the development version, so it will be fixed in the upcoming release. the development tarball (https://code.mpimet.mpg.de/attachments/22243/cdo-1.9.9rc4.tar.gz) should have all you need.

There is also a development package for cdo on conda forge, that has this fix. you can install it with

conda install --channel conda-forge/label/dev cdo

or with a new environment:

conda create --name cdo-dev --channel conda-forge/label/dev cdo

This binary is for shipping cdo with CMOR.

hth
ralf

RE: intlevelx3d - getting an error where once was none - Added by Brendan DeTracey almost 4 years ago

Thanks ralf! I'll give this a try. Do you have a guesstimate on the date of the next cdo release?

RE: intlevelx3d - getting an error where once was none - Added by Uwe Schulzweida almost 4 years ago

The next release is planned for October.

RE: intlevelx3d - getting an error where once was none - Added by Brendan DeTracey almost 4 years ago

Thanks Uwe

RE: intlevelx3d - getting an error where once was none - Added by Brendan DeTracey almost 4 years ago

A intlevelx3d bug has resurfaced. Reported here: https://code.mpimet.mpg.de/issues/9822

(1-5/5)

Project

General

Profile

CDO