Question for you NOAA model experts out there (boosts welcome):

what is an easy (*) way to plot the RAP 13km domain data (https://rapidrefresh.noaa.gov/
) on a map using python/xarray/cartopy?

The high res domains are pretty straight forward, but the outer domain with its unusual grid does not work out of the box.

I can use wgrib2 to extract the record I want, but it ends up flattening the data into 1D and has no coordinates.
wgrib2 to netcdf segfaults, even though they seem to have added some support for the 'weird NCEP grid'.

I've found this approach: https://www.ssec.wisc.edu/realearth/solved-north-american-domain-rap-13-5km-rotated-pole-grib2-for-gis/

but I don't really want to work with a) that toolchain b) TIFFs - if I can avoid it.

Any ideas welcome.

(*) easy as in the least number of tools possible and only common ones. If it could be done with just wgrib2 and some python packages, that would be ideal.

Edit: added hashtags
#NOAA #NCEP #WGRIB #WRF #GRIB2

Satisfaction comes from wherever you can get it. An issue I filed on Github in 2021 and which users gradually piled onto to say, oh yeah, we need that? It's resolved!

New release of my software coming in the next couple days.

#netcdf #grib2