sp2gp fails with Tco199 grid
Added by Klaus Wyser 10 months ago
We are using a Tco199 grid from ECMWF and I want transform a variable from spectral to gridpoint. This fails with
cdo sinfo -sp2gp,cubic spec-orog-tco199.grb cdo(1) sp2gp: Process started FFT does not work with len=816 (n=17)! cdo(1) sp2gp (Abort): FFT error!
sp2gp,cubic works fine for some other Tco grids that I have tested with, but not Tco199.
Can this be fixed somehow? Testfile is attached.
The CDO version that I'm using is pretty fresh:
Climate Data Operators version 2.3.0 (https://mpimet.mpg.de/cdo)
System: x86_64-pc-linux-gnu
CXX Compiler: mpiicpc -g -O2 -fopenmp -pthread
CXX version : icpc (ICC) 2021.9.0 20230302
C Compiler: mpiicc -g -O2 -fopenmp -pthread -pthread
C version : icc (ICC) 2021.9.0 20230302
F77 Compiler: f77 -g -O2
F77 version : GNU Fortran (GCC) 12.3.0
Features: 755GB 32threads c++17 OpenMP45 pthreads HDF5 NC4/HDF5 OPeNDAP udunits2 proj cmor sse2
Libraries: yac/3.0.1 NetCDF/4.9.2 HDF5/1.12.2 proj/7.2.0 cmor/3.4.0
CDI data types: SizeType=size_t
CDI file types: srv ext ieg grb1 grb2 nc1 nc2 nc4 nc4c nc5 nczarr
CDI library version : 2.3.0
cgribex library version : 2.1.1
ecCodes library version : 2.32.0
NetCDF library version : 4.9.2 of Feb 14 2024 15:24:19 $
exse library version : 1.5.0
FILE library version : 1.9.1
| spec-orog-tco199.grb (79.5 KB) spec-orog-tco199.grb |
Replies (4)
RE: sp2gp fails with Tco199 grid - Added by Klaus Wyser 10 months ago
Addendum: sp2gp also fails with a similar error with Tco159 or Tco399, but works fine with Tco95 or Tco319.
RE: sp2gp fails with Tco199 grid - Added by Klaus Wyser 10 months ago
And another intersting piece of information: sp2gp,cubic works fine with the Tco199 grid on my Windows laptop with cygwin, but not on our new HPC ;-)
cdo -V on my PC yields
Features: 15GB 8threads C++14 OpenMP45 Fortran PTHREADS HDF5 NC4/HDF5 OPeNDAP SZ UDUNITS2 PROJ XML2 CURL FFTW3 SSE2
and on the HPC
Features: 755GB 32threads c++17 OpenMP45 pthreads HDF5 NC4/HDF5 OPeNDAP udunits2 proj cmor sse2
The interesting thing is the FFTW3 library on the PC, could that make the difference?
RE: sp2gp fails with Tco199 grid - Added by Klaus Wyser 10 months ago
It was the missing FFTW3 library on our HPC. After rebuilding cdo with fftw everything works fine.
Issue can be closed.
RE: sp2gp fails with Tco199 grid - Added by Uwe Schulzweida 10 months ago
Hi Klaus,
Thank you for the information!
We will improve the error message in the next CDO version 2.4.1 and point out that FFTW3 is missing.
Cheers,
Uwe