clusters/halo.py: Add property "volume" (unit: [cm^3])

Signed-off-by: Aaron LI <aly@aaronly.me>

1 files changed, 9 insertions, 0 deletions

diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index 44d1e1a..d7b3ec4 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -163,6 +163,15 @@ class RadioHalo:
         return r_halo
 
     @property
+    def volume(self):
+        """
+        The halo volume, calculated from the above radius.
+        Unit: [cm^3]
+        """
+        r_cm = self.radius * AUC.kpc2cm
+        return (4*np.pi/3) * r_cm**3
+
+    @property
     def magnetic_field(self):
         """
         The magnetic field strength at the simulated observation