empty grids when regridding
Added by Sentia Goursaud almost 3 years ago
Hello,
I am trying to regrid a netcdf file (see the attached file infile.nc) provided with a 1800x900 lonlat grid (grid.dat), using remapbil (eg t63 or r96x73), but the resulting output shows a file with empty grids (see the file oroguk-20210623.png).
Does anyone have an idea for this bad regridding?
Best regards,
Sentia
| oroguk-20210623.png (74.5 KB) oroguk-20210623.png | |||
| grid.dat (5.8 MB) grid.dat | |||
| infile.nc (12.4 MB) infile.nc |
Replies (2)
RE: empty grids when regridding - Added by Karin Meier-Fleischer almost 3 years ago
Hi Sentia,
I would say the interpolation of your data with remapbil is correct. The data file infile.nc has 1434411 grid cells with missing values for each variable (Topo, lsm). If you interpolate the data from r1800x900 to a very coarse grid r96x73 the number of missing values are 6416 of 7008, only 592 cells with data.
Instead of remapbil you can use the operator remapdis https://code.mpimet.mpg.de/projects/cdo/embedded/cdo.pdf#subsection.2.12.4. The default number of nearest neighbours is 4. With the default value 4 the result is similar to remapbil. But you can use any number of neighbours.
cdo remapdis,r96x73,20 infile.nc outfile.nc
-Karin
RE: empty grids when regridding - Added by Sentia Goursaud almost 3 years ago
Thank you Karin. I did not know this command but it works !
Sentia