1
0
Fork 0
mirror of https://github.com/cosmo-sims/monofonIC.git synced 2024-09-19 17:03:45 +02:00

Fix constant used in the SWIFT i/o plugin to match exactly the value of the physical constants used internally by SWIFT

This commit is contained in:
Matthieu Schaller 2021-03-14 18:50:46 +01:00
parent 5207492b8e
commit 098dafa1ef

View file

@ -63,7 +63,7 @@ public:
num_ranks_ = 1;
real_t astart = 1.0 / (1.0 + cf_.get_value<double>("setup", "zstart"));
const double rhoc = 27.7519737; // in h^2 1e10 M_sol / Mpc^3
const double rhoc = 27.7536609198; // in h^2 1e10 M_sol / Mpc^3 assuming SWIFT's internal constants
hubble_param_ = pcc->cosmo_param_["h"];
zstart_ = cf_.get_value<double>("setup","zstart");
@ -297,6 +297,7 @@ public:
// now each node writes its own chunk in a round-robin fashion, appending at the end of the currently existing data
for (int rank = 0; rank < num_ranks_; ++rank) {
// wait until the initialisation or the previous rank in the loop is done
MPI_Barrier(MPI_COMM_WORLD);
if (rank == this_rank_) {