Project

General

Profile

Forums » Support »

proj_params

Added by Tony Payne almost 2 years ago

Hi

Can I check what the effect of grid_mapping and grid_mapping_name in a CDO grid description file is when I specify proj_params please? Is the projection in grid_mapping_name ignored in favour of the proj_params information? I tried omitting grid_mapping_name but got an error message suggesting that grid_mapping_name has to be specified.

gridtype = projection
gridsize = 168861
xsize = 301
ysize = 561
xunits = "m"
yunits = "m"
xfirst = -800000
xinc = 5000
yfirst = -3400000
yinc = 5000
grid_mapping = crs
grid_mapping_name = sinusoidal
proj_params='+proj=stere +lat_0=90 +ellps=WGS84 +x_0=0 +y_0=0 +lat_ts=71'

Sorry I am using version 1.9 because of compiler issues with more recent versions (hence cannot use grid_mapping_name = polar_stereographic).

Thanks
tony


Replies (1)

RE: proj_params - Added by Uwe Schulzweida almost 2 years ago

Hi Tony,

You are right, the grid_mapping_name should not be necessary if proj_params is defined. And this is also the case since at least CDO version 1.9.10. Since when this is exactly I can't say, because I don't have an older version installed.
You can check this with the operator topo and the option -v:

cdo -v -f nc topo,<myproj> outfile
This will create a field with topography data and the correct coordinates.

Cheers,
Uwe

    (1-1/1)