"""
Short example script to demonstrate how to read in FOF and Subhalo properties.
"""

import sim_tools as st       # Set of helper functions for reading in data
import numpy as np
from pdb import set_trace

# Specify directory of the simulation. Here, we use the Hydro version of CE-0
simdir = '/net/quasar/data3/Hydrangea/CE-0/HYDRO/'

# ... or, more easily, use a shortcut function:
# (returns exactly the same file path, but less effort!)
snapdir = st.form_files(simdir, isnap = 29, types = 'snap')


dm_coord = st.eagleread(snapdir, 'PartType1/Coordinates', astro = True)[0]