From 091c9dee757db707fe7232088f50c67b7e78b9ef Mon Sep 17 00:00:00 2001
From: Michael Buehlmann <buehlmann.michi@gmail.com>
Date: Mon, 18 Jul 2022 14:37:40 -0500
Subject: [PATCH] fix HACC masses for HACCHydro, single species

---
 src/plugins/output_genericio.cc | 3 +--
 1 file changed, 1 insertion(+), 2 deletions(-)

diff --git a/src/plugins/output_genericio.cc b/src/plugins/output_genericio.cc
index 9e4a850..773ca0b 100644
--- a/src/plugins/output_genericio.cc
+++ b/src/plugins/output_genericio.cc
@@ -74,8 +74,7 @@ public:
 
     void write_particle_data(const particle::container &pc, const cosmo_species &s, double Omega_species)
     {
-        double rhoc = 27.7519737; // in h^2 1e10 M_sol / Mpc^3
-        double boxmass = Omega_species * rhoc * std::pow(cf_.get_value<double>("setup", "BoxLength"), 3.);
+        double boxmass = Omega_species * munit_;
         double particle_mass = boxmass / pc.get_global_num_particles();
 
         size_t npart = pc.get_local_num_particles();