eccodes | ECMWF 's GRIB and BUFR | Bot library

Whenever I try importing cfgrib it gives me runtime error that it could not load ecCodes library

I'm trying to load GRIB2 files from DWD's ICON model using xarray and cfgrib. Most variables work fine, but for some (like for example CLCL), xarray.open_dataset raises the following error:

I am using ecCodes library in my project, and I have encountered an issue that memory is not freed between reading the files.

The minimal example representing the problem is this (and is basically a combination of those two library API usage examples [1](https://confluence.ecmwf.int/display/ECC/grib_get_keys) [2]:

I would like to read a GRIB file downloaded from server using ecCodes library in Rust. However, my current solution results in segmentation fault. The extracted example, replicating the problem, is below.

I download the file using reqwest crate and get the response as Bytes¹ using bytes(). To read the file with ecCodes I need to create a codes_handle using codes_grib_handle_new_from_file()², which as argument requires *FILE usually get from fopen(). However, I would like to skip IO operations. So I figured I could use libc::fmemopen() to get *FILE from Bytes. But when I pass the *mut FILE from fmemopen() to codes_grib_handle_new_from_file() segmentation fault occurs.

I suspect the issue is when I get from Bytes a *mut c_void required by fmemopen(). I figured I can do this like that:

I have accessed Anaconda Prompt for the environment that I will be working.

Once there, I have changed the directory to the place where I have the grib files that I will be working with, in order to install cfgrib, I have used the following command:

conda install -c conda-forge cfgrib

As the Python module depends on ECMWF ecCodes binary library, in the same Anaconda Prompt, I ran the following command:

conda install -c conda-forge eccodes

Then, opened the file in Spyder, and tried to run this simple script:

I have a netCDF file that consists of several meteorological variables (e.g. boundary layer height, surface temperature, surface pressure) measured once daily for a 10-year period over a region. As a result, there are over 4,000 time steps per variable.

I would like to use QGIS or Python to extract the value of each meteorological variable at several point locations for the entirety of the time series. The desired output is either (a) a new vector layer consisting of the point data and the extracted meteorological variable values over all time steps, or (b) a similar output, but in CSV/Excel format.

The "Add raster values to points" tool in QGIS performs the desired operation, but I can't apply this to my NC file since it is mesh data, not raster. Any help here would be appreciated. I'm also quite comfortable with Python is there's a Python package that could achieve this.

Here's the output of ncdump for my netCDF file: