0

I'm trying to get data from a GFS APCP (Total Precipitation 'tp' accum), but only can access to time 2. Time 1 have the accumulated precipitation from RUN hour to the forecast hour. For example for 06Z , time 1 shows tha acumulated precipitation from 00 to 06. Time 2 have the accumulated fomr the sinoptic hour until the forecast hour. For example, 17Z have the precipitation from 12 to 17. With: backend_kwargs = {'filter_by_keys': {'stepType':'accum'},'indexpath': ''} ds = xr.open_dataset(file_path , engine="cfgrib" , backend_kwargs=backend_kwargs) only can be read time 2, not time 1.

As GFS APCP grib comes with 2 GRIB forecast observation time as Mixed Intervals Accumulation, I wonder to read both. Is that possible? Attached picture from Panoply with both time 1 and 2 to be selected in the same grib.Panoply APCP

Also I tried to merge as merge but the seems that the second grib doesn't exists when cfgrib open_dataset. So can't merge.

Thank you

0 Answers0