Radar-derived precipitation rates unrecognized (NEXRAD data, MRMS reanalysis dataset)

Disclaimer: I’m a beginner to python and wrangling climate datasets in general. I did my best to scour the web and the issues page to find a solution with no luck. I hope I’ve provided enough information!

Data source: http://edc.occ-data.org/nexrad/mosaic/#data-updates

• Using ‘MRMS_PrecipRate_00.00_20010101-120000.grib2’ which is the first file in the 2001 dataset

Problem: The ‘PrecipRate’ attribute I know to be in the grib is unrecognized resulting in na data values. I confirmed that these data are there using the pygrib library, but I need this to work with cgrib so I can use the ‘open_mfdataset’ function.

Code:

ds = xr.open_dataset('MRMS_PrecipRate_00.00_20010101-120000.grib2', engine = 'cfgrib')

print(ds)

<xarray.Dataset>
Dimensions:     (latitude: 3500, longitude: 7000)
Coordinates:
    time        datetime64[ns] 2001-01-01T12:00:00
    step        timedelta64[ns] 00:00:00
    surface     float64 0.0
  * latitude    (latitude) float64 54.99 54.98 54.98 54.97 ... 20.03 20.02 20.01
  * longitude   (longitude) float64 230.0 230.0 230.0 ... 300.0 300.0 300.0
    valid_time  datetime64[ns] 2001-01-01T12:00:00
Data variables:
    unknown     (latitude, longitude) float32 ...
Attributes:
    GRIB_edition:            2
    GRIB_centre:             161
    GRIB_centreDescription:  161
    GRIB_subCentre:          0
    Conventions:             CF-1.7
    institution:             161
    history:                 2021-05-04T13:11 GRIB to CDM+CF via cfgrib-0.9.9...

print(ds['unknown'].values)

[[nan nan nan ... nan nan nan]
 [nan nan nan ... nan nan nan]
 [nan nan nan ... nan nan nan]
 ...
 [nan nan nan ... nan nan nan]
 [nan nan nan ... nan nan nan]
 [nan nan nan ... nan nan nan]]

[YasharGIS]

thank you Shahram! it worked.

[shahramn]

The fields in these files use a non-standard (i.e., local) encoding. For example the discipline is 209 which is not a WMO standard value