Missing values when using remapbil
Added by Fredrik Boberg over 3 years ago
Dear cdo experts,
I have (ERA5) data for a region around Greenland (=input) that I want to regrid using remapbil to a smaller region (=grid). I end up with the right output grid but where the upper corners have missing values (=output). Any suggestions on how I can end up with an output file without missing values is greatly appreciated.
The attached jpgs are very simple plots. Sorry about that!
Best regards,
Fredrik
input.jpg (573 KB) input.jpg | land sea mask from input file | ||
output.jpg (423 KB) output.jpg | land sea mask from output file | ||
grid.nc (2.77 MB) grid.nc | target grid file | ||
input.nc (612 KB) input.nc | input file | ||
output.nc (3.7 MB) output.nc | output file |
Replies (6)
RE: Missing values when using remapbil - Added by Fredrik Boberg over 3 years ago
More info:
I use cdo 1.9.10 and the command
cdo remapbil,grid.nc input.nc output.nc
/Fredrik
RE: Missing values when using remapbil - Added by Karin Meier-Fleischer over 3 years ago
Hi Fredrik,
you should use remapnn:
cdo -remapnn,grid.nc input.nc output.nc
output.png (156 KB) output.png |
RE: Missing values when using remapbil - Added by Fredrik Boberg over 3 years ago
Thanks Karin,
If you look at the sftlf variable you will see that remapnn is not working either (at least not for me). It will just extrapolate the values outside the "missing value corners".
Best,
Fredrik
RE: Missing values when using remapbil - Added by Karin Meier-Fleischer over 3 years ago
Let's have a look at your files.
The grid.nc file upper left and upper right corner is outside the input.nc grid.
During the interpolation to the grid.nc grid this cells of missing data will be filled with _FillValue since not present.
1. grid.nc
2. input.nc
3. output.nc
output.png (132 KB) output.png | |||
grid.png (184 KB) grid.png | |||
input.png (156 KB) input.png |
RE: Missing values when using remapbil - Added by Fredrik Boberg over 3 years ago
Dear Karin,
I'm so sorry for this! There was a typo in my extraction script. I have now put the correct longitude limits in. Many thanks for helping out a confused old man.
Best,
Fredrik
RE: Missing values when using remapbil - Added by Karin Meier-Fleischer over 3 years ago
You're welcome!