Problem with bottomvalue operator

Added by Brendan DeTracey over 3 years ago

Hi everyone!
I was lost in space but now I return for some help.

$ cdo -V
Climate Data Operators version 1.9.9 (https://mpimet.mpg.de/cdo)
System: x86_64-conda-linux-gnu
CXX Compiler: /home/conda/feedstock_root/build_artifacts/cdo_1604652667339/_build_env/bin/x86_64-conda-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_1604652667339/work=/usr/local/src/conda/cdo-1.9.9 -fdebug-prefix-map=/home/brendan/miniconda3/envs/cmip6=/usr/local/src/conda-prefix -fopenmp -pthread
CXX version : unknown
C Compiler: /home/conda/feedstock_root/build_artifacts/cdo_1604652667339/_build_env/bin/x86_64-conda-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_1604652667339/work=/usr/local/src/conda/cdo-1.9.9 -fdebug-prefix-map=/home/brendan/miniconda3/envs/cmip6=/usr/local/src/conda-prefix -fopenmp -pthread -pthread
C version : unknown
F77 Compiler: /home/conda/feedstock_root/build_artifacts/cdo_1604652667339/_build_env/bin/x86_64-conda-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_1604652667339/work=/usr/local/src/conda/cdo-1.9.9 -fdebug-prefix-map=/home/brendan/miniconda3/envs/cmip6=/usr/local/src/conda-prefix
F77 version : GNU Fortran (crosstool-NG 1.24.0.133_b0863d8_dirty) 9.3.0
Features: 62GB 24threads C++11 OpenMP45 Fortran PTHREADS HDF5 NC4/HDF5/threadsafe OPeNDAP UDUNITS2 PROJ XML2 CURL FFTW3 SSE3
Libraries: HDF5/1.10.6 proj/7.2.0 xml2/2.9.10 curl/7.71.1
Filetypes: srv ext ieg grb1 grb2 nc1 nc2 nc4 nc4c nc5
     CDI library version : 1.9.9
 cgribex library version : 1.9.5
 ecCodes library version : 2.19.1
  NetCDF library version : 4.7.4 of Oct 21 2020 19:55:48 $
    hdf5 library version : 1.10.6 threadsafe
    exse library version : 1.4.1
    FILE library version : 1.9.1

Command line:

cdo bottomvalue thetao_Omon_MRI-ESM2-0_historical_r1i1p1f1_gn_195001_1-60.nc thetao_Omon_MRI-ESM2-0_historical_r1i1p1f1_gn_195001_bottom.nc

The result is... strange. This is my last little twiddle for a paper. And I have properly cited cdo because cdo is such an awesome tool.
Here is the screenshot. I know something is wrong with the input file, but I can't find it...

Ich liebe cdo! Ich liebe dich!
Brendan
EDIT: oops.. forgot the file attachents.. I am out of practise.

Download all files

Screenshot 2021-09-21 115015.jpg (133 KB) Screenshot 2021-09-21 115015.jpg
thetao_Omon_MRI-ESM2-0_historical_r1i1p1f1_gn_195001_bottom.nc (10.5 MB) thetao_Omon_MRI-ESM2-0_historical_r1i1p1f1_gn_195001_bottom.nc
thetao_Omon_MRI-ESM2-0_historical_r1i1p1f1_gn_195001_1-60.nc (41.4 MB) thetao_Omon_MRI-ESM2-0_historical_r1i1p1f1_gn_195001_1-60.nc

Replies (3)

RE: Problem with bottomvalue operator - Added by Ralf Mueller over 3 years ago

hi Brendan (ich liebe dich auch ;-))

something is really off with your input. I cannot change the missval to a constant and get an internal error when using setmisstonn.

Maybe Uwe can find out more

cheers
ralf

RE: Problem with bottomvalue operator - Added by Uwe Schulzweida over 3 years ago

The missing values in the file are NaN and the corresponding attribute has the value 1.e+20f. Since CDO always uses the missing_value attribute, all calculations with NaN lead to unpredictable results. You can either set all NaNs to missing_value or set the missing_value to NaN:

cdo bottomvalue -setctomiss,NaN  infile  outfile
cdo bottomvalue -setmissval,NaN  infile  outfile

Cheers,
Uwe

RE: Problem with bottomvalue operator - Added by Brendan DeTracey over 3 years ago

Thanks Guys... My project is finally wrapping up! Feels good man! Real good! Press

!{EJECT}!     !!{HARDER}!! .....   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!{....................}!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!    ......     !^∞{.^∞}!^∞

(blip)

(1-3/3)

Project

General

Profile

CDO