Re: [gembud] ECMWF grib2 files from WMO Essential no longer being decoded properly by NAWIPS or AWIPS

  • To: "Manousos, Peter C" <pmanousos@xxxxxxxxxxxxxxxxxxx>
  • Subject: Re: [gembud] ECMWF grib2 files from WMO Essential no longer being decoded properly by NAWIPS or AWIPS
  • From: Mike Zuranski <mzuranski@xxxxxxxx>
  • Date: Thu, 29 Jun 2023 12:33:41 -0500
Greetings,

Long story short, the ECMWF open data is now using a different compression
that g2clib can't seem to handle.  Luckily, the work around that's worked
for me so far is to pass it through wgrib2 to change the compression type,
after that other software can read it.

A good thread over on the GrADS email list explains what's happening in
more detail:
http://gradsusr.org/pipermail/gradsusr/2023-June/043431.html

As for AWIPS/EDEX, I'm going to reach out to the AWIPS team here and see if
we can do something to make that process easier.


Best,
-Mike


*Mike Zuranski*
Data Engineer II
Unidata Program Center
University Corporation for Atmospheric Research


On Thu, Jun 29, 2023 at 12:27 PM Manousos, Peter C via gembud <
gembud@xxxxxxxxxxxxxxxx> wrote:

> Greetings.  As of 06 28 2023 12z it seems the ECMWF’s freely available WMO
> Essential grib2 data has changed (here
> <https://data.ecmwf.int/forecasts/20230628/12z/0p4-beta/oper/>).   It is
> no longer fully processable by dcgrib2 or AWIPS EDEX
>
>
>
> We did not see any change notice on ECMWF web site or Google and were
> wondering if anyone else is experiencing the following (attached is a
> one-field <500 mb HGT> grib2 file), and if you did would you happen to know
> the cause and any workarounds?
>
>
>
> wgrib2 -v -V test.grb2  shows what looks like a viable grib2 file
>
>
>
> 1:0:vt=2023062803:500 mb:3 hour fcst:HGT Geopotential Height [gpm]:
>
>     ndata=405900:undef=0:mean=5562.91:min=4689.51:max=5955.01
>
>     grid_template=0:winds(N/S):
>
>       lat-lon grid:(900 x 451) units 1e-06 input WE:NS output WE:SN res 48
>
>       lat 90.000000 to -90.000000 by 0.400000
>
>       lon 180.000000 to 179.600000 by 0.400000 #points=405900
>
>       basic_ang=0 sub_angle=1000000 units=0.000001
>
>       unscaled lat=90000000 to -90000000 lon=180000000 to 179600000
>
>
>
> As additional info a wgrib2 -v –V on a grib2 file previous to June 28
> shows the same thing as above.
>
>
>
> However, NAWIPS 7.14.0.1   dcgrib2 -d temp.log -v 100 test.grd <
> test.grb2  gives the following error (in bold)
>
>
>
> [27248] 230628/1520[DCINIT 3] Version 7.14.0.1
>
> [27248] 230628/1520[DC 2]  read 102399/102399 bytes strt 0 newstrt 102399
>
> [27248] 230628/1520[DC 2]  read 1/1 bytes strt 102399 newstrt 102400
>
> [27248] 230628/1520[DC 2]  read 83204/102397 bytes strt 0 newstrt 83204
>
> *[27248] 230628/1520[CED_GNUM -1] No GID KXKY [900 451] LL [90.000000
> 180.000000] UR [-90.000000 179.600006]*
>
> *[27248] 230628/1520[CED_GNUM -1] No GID KXKY [900 451] LL [90.000000
> 180.000000] UR [-90.000000 179.600006]*
>
> [27248] 230628/1520[DC 2]  read 0/19196 bytes strt 83204 newstrt 83204
>
> [27248] 230628/1520[DCGRIB -9] End of Data
>
> [27248] 230628/1520[DC 5]  Normal termination.
>
> [27248] 230628/1520[DC 2]  Number of bulletins read and processed: 1
>
> [27248] 230628/1520[DC 6]  Shutting down.
>
>
>
> And nagrib2 also gives an error (in bold)
>
>
>
> getdrstemplate: DRS Template 5.42 not defined.
>
> *[NAGRIB2 -25]  Error unpacking GRIB2 metadata*
>
>          0 GRIB messages were read or scanned from the GRIB file:
>
>
> test.grb2
>
>          0 grids were written to the GEMPAK file:
>
>
> test.grd
>
> Parameters requested: GBFILE,GDOUTF,PROJ,GRDAREA,KXKY,MAXGRD,CPYFIL,GAREA,
>
> OUTPUT,G2TBLS,G2DIAG,OVERWR,PDSEXT.
>
>
>
>
>
> If you have AWIPS EDEX – submitting to EDEX shiows the following error (at
> the end in bold with preceding traceback messages)
>
>
>
> ERROR 2023-06-28 19:02:42,848 4308 [Ingest.GribDecode-6] grib: EDEX -
> Error processing file
> [/awips2/data_store/manual/grib/20230628/19/test.grb2]:
>
> Traceback (most recent call last):
>
>   File
> "/awips2/edex/lib/plugins/com.raytheon.edex.plugin.grib.jar/GribDecoder.py",
> line 184, in decode
>
>     gribDictList = grib2.decode(gribFile, self.startPosition,
> self.messageLength)
>
> *error: Failed to get grib field.*
>
>
>
> ------------------------------
>
> * The information contained in this message is intended only for the
> personal and confidential use of the recipient(s) named above. If the
> reader of this message is not the intended recipient or an agent
> responsible for delivering it to the intended recipient, you are hereby
> notified that you have received this document in error and that any review,
> dissemination, distribution, or copying of this message is strictly
> prohibited. If you have received this communication in error, please notify
> us immediately, and delete the original message. *
> _______________________________________________
> NOTE: All exchanges posted to Unidata maintained email lists are
> recorded in the Unidata inquiry tracking system and made publicly
> available through the web.  Users who post to any of the lists we
> maintain are reminded to remove any personal information that they
> do not want to be made public.
>
>
> gembud mailing list
> gembud@xxxxxxxxxxxxxxxx
> For list information or to unsubscribe,  visit:
> https://www.unidata.ucar.edu/mailing_lists/
>
  • 2023 messages navigation, sorted by:
    1. Thread
    2. Subject
    3. Author
    4. Date
    5. ↑ Table Of Contents
  • Search the gembud archives: