Project

General

Profile

Forums » Support »

Remapping with remapbil aborts because source grid rank !=2

Added by Rianne Giesen almost 3 years ago

Hello,

I am trying to remap an ECMWF grib file with a reduced Gaussian grid to a regular 0.25 degree lat,lon grid and save to netCDF.

I tried to achieve this with a bilinear interpolation in CDO: 'cdo remapbil,global_0.25 in.grib out.nc' but get the following messages:
cdo remapbil: Bilinear weights from gaussian_reduced (542080) to lonlat (1440x720) grid
cdo remapbil (Abort): Can't do bilinear interpolation when source grid rank != 2

Is there a simple way to preprocess the input file to adjust the rank of the source grid?

I have attached an example grib file with one field. Below is my CDO version information.

Thanks in advance!
Rianne

CDO version:
Climate Data Operators version 1.9.10 (https://mpimet.mpg.de/cdo)
System: x86_64-pc-linux-gnu
CXX Compiler: g++ -g -O2 -fopenmp -pthread
CXX version : g++ (GCC) 10.2.1 20201125 (Red Hat 10.2.1-9)
C Compiler: gcc -DHAVE_LIBNC_DAP -fopenmp -pthread -pthread
C version : gcc (GCC) 10.2.1 20201125 (Red Hat 10.2.1-9)
F77 Compiler: gfortran -fPIC
F77 version : GNU Fortran (GCC) 10.2.1 20201125 (Red Hat 10.2.1-9)
Features: 31GB 8threads C++14 OpenMP45 Fortran PTHREADS HDF5 NC4/HDF5 OPeNDAP SZ UDUNITS2 PROJ XML2 CURL FFTW3 SSE2
Libraries: HDF5/1.10.5 proj/6.3.2 xml2/2.9.10 curl/7.69.1
Filetypes: srv ext ieg grb1 grb2 nc1 nc2 nc4 nc4c nc5
CDI library version : 1.9.10
ecCodes library version : 2.20.0
NetCDF library version : 4.7.0 of Jan 29 2020 00:00:00 $
hdf5 library version : 1.10.5
exse library version : 1.4.2
FILE library version : 1.9.1

nwp_2019032_06_06_u10n.grib (795 KB) nwp_2019032_06_06_u10n.grib

Replies (4)

RE: Remapping with remapbil aborts because source grid rank !=2 - Added by Uwe Schulzweida almost 3 years ago

The native Gaussian reduced grid is not supported directly for remapping in CDO. This grid is handled like an unstructured grid. Therefor only the nearest neightbor (remapnn) or inverted distants method (remapdis) can be used.
For remapbil, the data must first be transformed to a regular grid:

cdo remapbil,global_0.25 -setgridtype,regular in.grib out.nc

RE: Remapping with remapbil aborts because source grid rank !=2 - Added by Rianne Giesen almost 3 years ago

Thanks Uwe, for your explanation and providing the options that work.

RE: Remapping with remapbil aborts because source grid rank !=2 - Added by Etienne Tourigny about 2 years ago

Are there any plans to add support for remapping to the reduced gaussian grids? What are the alternatives currently?

RE: Remapping with remapbil aborts because source grid rank !=2 - Added by Uwe Schulzweida about 2 years ago

For remapping to reduced gaussian grids, remapbil, remapbic, remapnn and remapdis can be used. For remapping between reduced gaussian grids only remapnn and remapdis can be used.

    (1-4/4)