Skip to content

Missing plotting support for gridType="lambert" #122

@uandrae

Description

@uandrae

Is your feature request related to a problem? Please describe.

The plotting, as in quickplot() works well for gridtype lambert_lam whereas the pure lambert projection, https://codes.ecmwf.int/grib/format/grib2/templates/3/30/, returns the following error

File .../lib/python3.10/site-packages/earthkit/data/core/fieldlist.py:336, in Field.to_points(self, flatten, dtype, index)
    334     return dict(x=lon, y=lat)
    335 else:
--> 336     raise ValueError("to_points(): geographical coordinates in original CRS are not available")
    338 # convert values to array format
    339 assert r

ValueError: to_points(): geographical coordinates in original CRS are not available

Applying grib_set -s gridType=lambert_lam on the same file makes the plotting work.

Describe the solution you'd like

No response

Describe alternatives you've considered

No response

Additional context

No response

Organisation

No response

Metadata

Metadata

Assignees

Labels

enhancementNew feature or request

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions