Stand-alone FoF for non-cosmological snapshots
I want to run the stand-alone FoF on a planetary snapshot. The first, minor, issue is that the code currently requires several cosmological parameters to be in the .yml file.
The bigger problem is that although the printed output appears sensible in terms of the number of groups and their masses, the .dat catalogue file contains the group sizes and masses, but the max densities are all 0, max density index all -3, particle ID all -1, and Particle Density all blank. More importantly, a snapshot is also saved, but it only contains the various header information attributes like the units. So there is no actual output of which particles are in which fof groups.
Also, the snapshot name (after the base) was _-2147483648.hdf5
. I assume it should ideally add the new group IDs to the existing input snapshot rather than writing a new one?
Configured and run on cosma7 with ./configure --with-hydro=planetary --with-equation-of-state=planetary --with-gravity=with-potential --with-tbbmalloc --with-parmetis --enable-ipo --enable-stand-alone-fof
../../swiftsim/examples/fof --hydro -t 12 pE_T10_Th_T10_B45v100_1e6.yml
Parameter file: http://astro.dur.ac.uk/~cklv53/files/pE_T10_Th_T10_B45v100_1e6.yml
Example snapshot: http://astro.dur.ac.uk/~cklv53/files/pE_T10_Th_T10_B45v100_1e6_100000.hdf5