sp2gp fails with Tco199 grid
Added by Klaus Wyser 18 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 18 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 18 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 10 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 2 days 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