Project

General

Profile

Remap LonLat Grid to Lambert Conformal Projection with constant 5.5km spacing, but retain LatLon Information?

Added by Simon Hentschel 6 months ago

Hello,

I'm currently trying to remap ERA5 data to a Lambert Conformal Projection with 5.5km spacing, similar to that of CERRA, but moved to cover North America instead.
For this I'm using the "cdo remapbic" command with the grid-description from below.

While the reprojection itself works fine, for further processing I need a mapping from each grid cell to the coordinates it represents.
The converted file itself sadly only contains the x/y variables with the position encoded in 'meters' from the first grid point, instead of the lat/lon values.

Is there any way to retain the lat,lon coordinate map (similar to the original CERRA files)?

Thank you very much for your help!

Simon

#
  1. gridID similar to CERRA (moved to cover North America) #
    gridtype = projection
    gridsize = 1142761
    xsize = 1069
    ysize = 1069
    xname = x
    xunits = "m"
    yname = y
    yunits = "m"
    xfirst = 0
    xinc = 5500
    yfirst = 0
    yinc = 5500
    grid_mapping = Lambert_Conformal
    grid_mapping_name = lambert_conformal_conic
    standard_parallel = 50.
    earth_radius = 6371229.
    longitude_of_central_meridian = -94.514057
    latitude_of_projection_origin = 50.
    longitudeOfFirstGridPointInDegrees = -120
    latitudeOfFirstGridPointInDegrees = 20.292281
    longitudeOfSouthernPoleInDegrees = 0.
    latitudeOfSouthernPoleInDegrees = -90.

Replies (1)

RE: Remap LonLat Grid to Lambert Conformal Projection with constant 5.5km spacing, but retain LatLon Information? - Added by Uwe Schulzweida 6 months ago

Hello Simon,

You get lon/lat values if you add setgridtype,curvilinear after the interpolation:

cdo -setgridtype,curvilinear -remapbil,<targetgrid>  infile outfile
Cheers,
Uwe

    (1-1/1)