diff --git a/src/black_holes/EAGLE/black_holes.h b/src/black_holes/EAGLE/black_holes.h
index 468a72431d992e0890dbde152c2b4841bc3a94e4..8fdafadb9d53ce662c6b6d13c5db43d12c724c5e 100644
--- a/src/black_holes/EAGLE/black_holes.h
+++ b/src/black_holes/EAGLE/black_holes.h
@@ -23,6 +23,7 @@
 #include "black_holes_properties.h"
 #include "black_holes_struct.h"
 #include "cooling.h"
+#include "cooling/PS2020/cooling_tables.h"
 #include "cosmology.h"
 #include "dimension.h"
 #include "exp10.h"
diff --git a/src/cooling/EAGLE/cooling.c b/src/cooling/EAGLE/cooling.c
index 750934d70690b104cb76fcd60c860153b78cb5cd..90a977e4c1961e5cb2bbed7b9459cdededbeed2c 100644
--- a/src/cooling/EAGLE/cooling.c
+++ b/src/cooling/EAGLE/cooling.c
@@ -693,6 +693,31 @@ float cooling_get_temperature(
   return exp10(log_10_T);
 }
 
+/**
+ * @brief Compute the electron number density of a #part based on the cooling
+ * function.
+ *
+ * Does not exist in this model. We return 0.
+ *
+ * @param phys_const #phys_const data structure.
+ * @param hydro_props The properties of the hydro scheme.
+ * @param us The internal system of units.
+ * @param cosmo #cosmology data structure.
+ * @param cooling #cooling_function_data struct.
+ * @param p #part data.
+ * @param xp Pointer to the #xpart data.
+ */
+double cooling_get_electron_density(const struct phys_const *phys_const,
+                                    const struct hydro_props *hydro_props,
+                                    const struct unit_system *us,
+                                    const struct cosmology *cosmo,
+                                    const struct cooling_function_data *cooling,
+                                    const struct part *p,
+                                    const struct xpart *xp) {
+  return 0.;
+}
+
+
 /**
  * @brief Compute the electron pressure of a #part based on the cooling
  * function.
diff --git a/src/cooling/EAGLE/cooling.h b/src/cooling/EAGLE/cooling.h
index a44b219d9e65bc61c184f251d3b4ba2d34181402..c95964eb6dd5fca78891a3ed43ec4388ad321699 100644
--- a/src/cooling/EAGLE/cooling.h
+++ b/src/cooling/EAGLE/cooling.h
@@ -60,6 +60,14 @@ float cooling_timestep(const struct cooling_function_data *cooling,
                        const struct hydro_props *hydro_props,
                        const struct part *p, const struct xpart *xp);
 
+double cooling_get_electron_density(const struct phys_const *phys_const,
+                                    const struct hydro_props *hydro_props,
+                                    const struct unit_system *us,
+                                    const struct cosmology *cosmo,
+                                    const struct cooling_function_data *cooling,
+                                    const struct part *p,
+                                    const struct xpart *xp);
+
 double cooling_get_electron_pressure(
     const struct phys_const *phys_const, const struct hydro_props *hydro_props,
     const struct unit_system *us, const struct cosmology *cosmo,
@@ -170,7 +178,4 @@ void cooling_print_backend(const struct cooling_function_data *cooling);
 
 void cooling_clean(struct cooling_function_data *data);
 
-/*! Stub defined to let the BH model compile */
-#define colibre_cooling_N_elementtypes 1
-
 #endif /* SWIFT_COOLING_EAGLE_H */
diff --git a/src/extra_io/EAGLE/extra.h b/src/extra_io/EAGLE/extra.h
index 95698eca08b7380146f7d29ce177d6094a72e7a7..39839a6afd9f98953c4dd149483c6470f475ee76 100644
--- a/src/extra_io/EAGLE/extra.h
+++ b/src/extra_io/EAGLE/extra.h
@@ -23,6 +23,7 @@
 #include "cooling.h"
 #include "engine.h"
 #include "star_formation.h"
+#include "cooling/PS2020/cooling_tables.h"
 
 #define xray_table_date_string 20230110