sp2gp fails with Tco199 grid
Added by Klaus Wyser 17 days 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 17 days 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 17 days 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 9 days 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 1 day 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