gadget4 reading
It turns out that the header in Gadget-4 is slightly different: it does not contain the attribute "NumPart_Total_HighWord", which used to be present in the headers of Gadget-2/3. And apparently glnemo2 expects to find this attribute, in order to recognize it as "Gadget3 HDF5".
-- The solutions were:
- A student of mine, Gustavo, made a small modification in the Gadget-4 source code. In the file gadget4/src/io/snap_io.cc (line 887) he added the line
write_vector_attribute(handle, "NumPart_Total_HighWord", header.npartTotal, H5T_NATIVE_UINT64, NTYPES);
Now the HDF5 saved by this modified Gadget-4 can be read succesfully by glnemo2.
- Alternatively, in glnemo2/plugins/gadgeth5/gadgeth5.cc I just commented out line 79, which was
header.NumPart_Total_HighWord = getAttribute("NumPart_Total_HighWord");