Project

General

Profile

Interpolate a regional dataset

Added by Noelia Otero almost 10 years ago

Hi,

I am trying to interpolate a 1º grid into a coarser and regular grid 2.5º. The problem I have is that the original file only contains only 39 lat and 48 lon.
What I did so far (I don't know if it makes sense, but I had always used the cdo for interpolating global datasets, not one particular region):

cdo sellonlatbox,347,34,70,35 -remapbil,r144x73 file.nc newgrid.nc

This is the info of the newgrid.nc file:

File format: netCDF
-1 : Institut Source Param Time Typ Grid Size Num Levels Num
1 : unknown unknown -1 var F64 285 1 1 1
Horizontal grids :
1 : lonlat > size : dim = 285 nlon = 19 nlat = 15
lon : first = -12.5 last = 32.5 inc = 2.5 degrees_east
lat : first = 35 last = 70 inc = 2.5 degrees_north
Vertical grids :
1 : surface : 0
Time axis : 3653 steps
RefTime = 2000-01-01 00:00:00 Units = days Calendar = STANDARD
YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss
2000-01-01 00:00:00 2000-01-02 00:00:00 2000-01-03 00:00:00 2000-01-04 00:00:00

When I try to look into the data the file contains NA in all grids, so it is not working...
what could be the proper way to get the interpolation? Any suggestion?

Many thanks!

cheers,

Noelia


Replies (2)

RE: Interpolate a regional dataset - Added by Jaison-Thomas Ambadan almost 10 years ago

create a template of the region and use that for remap, like,

cdo -f nc -sellonlatbox,34,74,70,35 -remapcon,r144x73 -topo template.nc
cdo -f nc -remapbil,template.nc ifile.nc newgrid.nc

Cheers,
J

RE: Interpolate a regional dataset - Added by Noelia Otero almost 10 years ago

Thanks J.

It worked perfectly!

Cheers,

Noelia

    (1-2/2)