Segmentation fault (core dumped) with MeteoFrance data

[patrickhoebeke]

Hello,

I am currently struggling to read grib files provided by MeteoFrance (forecast data).
These files are publicly available at:
http://dcpc-nwp.meteo.fr/services/PS_GetCache_DCPCPreviNum?token=__5yLVTdr-sGeHoPitnFc7TZ6MhBcJxuSsoZp6y0leVHU__&model=AROME&grid=0.01&package=HP1&time=00H&referencetime=2017-12-12T03:00:00Z&format=grib2
If the link does not work, simply change the date at the end of the url (the files are regularly updated)

Example of code:

import pygrib
fp='/data/AROME_0.01_HP1_00H_201712110600.grib2'
with pygrib.open(fp) as grbs:
grbs = pygrib.open(fp)
for grb in grbs:
    if (grb['level']==100) & (grb['shortName']==u'u'):
        break
x= grb.values
grbs.close()

running the code above in a script with the following command:
gdb -ex r --args /usr/bin/python wms.py

I got:

GNU gdb (Ubuntu 7.11.1-0ubuntu1~16.5) 7.11.1
Copyright (C) 2016 Free Software Foundation, Inc.
License GPLv3+: GNU GPL version 3 or later http://gnu.org/licenses/gpl.html
This is free software: you are free to change and redistribute it.
There is NO WARRANTY, to the extent permitted by law. Type "show copying"
and "show warranty" for details.
This GDB was configured as "x86_64-linux-gnu".
Type "show configuration" for configuration details.
For bug reporting instructions, please see:
http://www.gnu.org/software/gdb/bugs/.
Find the GDB manual and other documentation resources online at:
http://www.gnu.org/software/gdb/documentation/.
For help, type "help".
Type "apropos word" to search for commands related to "word"...
Reading symbols from /usr/bin/python...(no debugging symbols found)...done.
Starting program: /usr/bin/python wms.py
[Thread debugging using libthread_db enabled]
Using host libthread_db library "/lib/x86_64-linux-gnu/libthread_db.so.1".

Program received signal SIGSEGV, Segmentation fault.
0x00007ffff669cf45 in jpc_qmfb_join_row () from /usr/local/gribapi/lib/libgrib_api.so

I am working with version pygrib version 2.0.2

As it looks like related to gribapi, I've recompiled it using the latest version provided by ECMWF. But I still get the same issue.... :-(

Where is the bug then ?

Any help would be more than welcome.

Patrick

[poubelle83]

Hello Patrick,
the problem may be caused by Jasper. I have encountered the same (and for reading meteofrance files) and resolve it by using (from sources) : jasper 2.0.14 (installing in /usr/local/jasper_new and building as static) eccodes-2.7.3 (don't forget to ask for python interface and to give path to jasper lib and includes). This message is old but hope this helps.