Matching haloes

[VictorForouhar]

Script is intended to be used to match haloes across simulations/halo finder runs using SOAP membership files.

[VictorForouhar]

I have basic matching done. I still need to add an initial sorting step in ascending particle ID, to handle different versions of the same simulation box (e.g. DMO vs HYDRO and HYDRO variations). I also want to make MPI rank sorting of subhaloes to match to be ordered and continuous in value.

[VictorForouhar]

This is the script I am using to test the implementation (HBT to VR central subgroups)

#!/bin/bash -l

module purge
module load python/3.12.4 gnu_comp/14.1.0 openmpi/5.0.3 parallel_hdf5/1.12.3
source /cosma8/data/dp004/dc-foro1/SOAP/openmpi-5.0.3-hdf5-1.12.3-env/bin/activate

membership_one="/cosma8/data/dp004/dc-foro1/halo_finder_comparison/FLAMINGO/L1000N0900/DMO_FIDUCIAL/SOAP_uncompressed/HBTplus/membership_0077/membership_0077.{file_nr}.hdf5"
membership_two="/cosma8/data/dp004/dc-foro1/halo_finder_comparison/FLAMINGO/L1000N0900/DMO_FIDUCIAL/SOAP_uncompressed/VR/membership_0077/membership_0077.{file_nr}.hdf5"
output_path="/cosma8/data/dp004/dc-foro1/halo_finder_comparison/matching_tests/output.hdf5"
                                                                                                                                                                                    
mpirun -np 8 -- python3 -m mpi4py match_halos.py $membership_one $membership_two $output_path --centrals_only --dmo --types 1