diff --git a/README.md b/README.md
index e91ea35..faaacd6 100644
--- a/README.md
+++ b/README.md
@@ -1,8 +1,9 @@
 # CSiBORG tools
 
 ## :scroll: Short-term TODO
-- [ ] Compare empirical $M_{500c}$ to the NFW expectation.
+- [x] Compare empirical $M_{500c}$ to the NFW expectation.
 - [ ] Calculate catalogues for all realisations.
+- [x] Add shortcut function for loading a catalogue
 
 
 ## :hourglass: Long-term TODO
diff --git a/scripts/utils.py b/scripts/utils.py
index 9131bde..874baa3 100644
--- a/scripts/utils.py
+++ b/scripts/utils.py
@@ -18,6 +18,7 @@ Notebook utility functions.
 
 
 import numpy
+from os.path import join
 from tqdm import trange
 from astropy.cosmology import FlatLambdaCDM
 
@@ -66,6 +67,27 @@ def load_mmains(N=None, verbose=True):
     return out
 
 
+def load_processed(Nsim, Nsnap):
+    simpath = csiborgtools.io.get_sim_path(Nsim)
+    outfname = join(
+        dumpdir, "ramses_out_{}_{}.npy".format(str(Nsim).zfill(5),
+                                               str(Nsnap).zfill(5)))
+    data = numpy.load(outfname)
+    # Add mmain
+    mmain = csiborgtools.io.read_mmain(Nsim, "/mnt/zfsusers/hdesmond/Mmain")
+    data = csiborgtools.io.merge_mmain_to_clumps(data, mmain)
+    # Cut on numbre of particles and finite m200
+    data = data[(data["npart"] > 100) & numpy.isfinite(data["m200"])]
+
+    # Do unit conversion
+    boxunits = csiborgtools.units.BoxUnits(Nsnap, simpath)
+    convert_cols = ["m200", "m500", "totpartmass", "mass_mmain",
+                    "r200", "r500", "Rs", "rho0", "peak_x", "peak_y", "peak_z"]
+    data = csiborgtools.units.convert_from_boxunits(
+        data, convert_cols, boxunits)
+    return data
+
+
 def load_planck2015(max_comdist=214):
     cosmo = FlatLambdaCDM(H0=70.5, Om0=0.307, Tcmb0=2.728)
     fpath = ("/mnt/zfsusers/rstiskalek/csiborgtools/"