CDO

Hi,

I am remapping a netcdf using the python interface of cdo. I would like to set the EPSG of the projection (EPSG:3035) in the target grid so that when I load the remapped netcdf in a GIS the projection is automatically recognized. In the grid file I have used the PROJ parameters for EPSG:3035, but it seems this is not enough to make the EPSG recognizable in a GIS. I was able to set the EPSG for the remapped netcdf with rioxarray, but I would like to know if there is a way to specify this in the grid file directly so I don't depend on rioxarray. Below is the grid file that I used, and in attachment you can find the complete workflow including the python script.

gridtype = projection
xsize = 14
ysize = 13
xunits = "meters"
yunits = "meters"
xfirst = 4077000
xinc = 3000
yfirst = 3019000
yinc = 3000
grid_mapping = crs
grid_mapping_name = lambert_azimuthal_equal_area
proj_params = "+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs"