issues: 96185559

This data as json

id node_id number title user state locked assignee milestone comments created_at updated_at closed_at author_association active_lock_reason draft pull_request body reactions performed_via_github_app state_reason repo type
96185559 MDU6SXNzdWU5NjE4NTU1OQ== 484 segfault with hdf4 file 1197350 closed 0     5 2015-07-20T23:15:06Z 2015-07-21T02:34:16Z 2015-07-21T02:34:16Z MEMBER      

I am trying to read data from the NASA MERRA reanalysis. An example file is: ftp://goldsmr3.sci.gsfc.nasa.gov/data/s4pa/MERRA/MAI3CPASM.5.2.0/2014/01/MERRA300.prod.assim.inst3_3d_asm_Cp.20140101.hdf The file format is hdf4 (NOT hdf5). (full file specification)

This file can be read by netCDF4.Dataset

python from netCDF4 import Dataset fname = 'MERRA300.prod.assim.inst3_3d_asm_Cp.20140101.hdf' nc = Dataset(fname) nc.variables['SLP'][0]

No errors

However, with xray

python import xray ds = xray.open_dataset(fname)

I get a segfault.

Is this behavior unique to my system? Or is this a reproducible bug?

Note: I am not using anaconda's netCDF package, because it does not have hdf4 file support. I had my sysadmin build us a custom netcdf and netCDF4 python.

 
{
    "url": "https://api.github.com/repos/pydata/xarray/issues/484/reactions",
    "total_count": 0,
    "+1": 0,
    "-1": 0,
    "laugh": 0,
    "hooray": 0,
    "confused": 0,
    "heart": 0,
    "rocket": 0,
    "eyes": 0
}
   completed 13221727 issue

Links from other tables