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

Extended all the header arrays to have 7 entries in the SWIFT output plugin

This commit is contained in:
Matthieu Schaller 2021-02-28 17:49:24 +01:00
parent e3606275ad
commit e9ebdc9e5d

View file

@ -22,9 +22,9 @@
#include "HDF_IO.hh" #include "HDF_IO.hh"
template <typename T> template <typename T>
std::vector<T> from_6array(const std::array<T,6>& a) std::vector<T> from_7array(const std::array<T,7>& a)
{ {
return std::vector<T>{{a[0], a[1], a[2], a[3], a[4], a[5]}}; return std::vector<T>{{a[0], a[1], a[2], a[3], a[4], a[5], a[6]}};
} }
template <typename T> template <typename T>
@ -45,8 +45,8 @@ protected:
bool blongids_, bdobaryons_; bool blongids_, bdobaryons_;
std::string this_fname_; std::string this_fname_;
std::array<uint32_t,6> npart_, npartTotal_, npartTotalHighWord_; std::array<uint32_t,7> npart_, npartTotal_, npartTotalHighWord_;
std::array<double,6> mass_; std::array<double,7> mass_;
double time_; double time_;
public: public:
@ -119,10 +119,10 @@ public:
HDFWriteGroupAttribute(fname_, "Header", "Dimension", 3); HDFWriteGroupAttribute(fname_, "Header", "Dimension", 3);
HDFWriteGroupAttribute(fname_, "Header", "BoxSize", lunit_); // in Mpc, not Mpc/h HDFWriteGroupAttribute(fname_, "Header", "BoxSize", lunit_); // in Mpc, not Mpc/h
HDFWriteGroupAttribute(this_fname_, "Header", "NumPart_Total", from_6array<unsigned>(npartTotal_)); HDFWriteGroupAttribute(this_fname_, "Header", "NumPart_Total", from_7array<unsigned>(npartTotal_));
HDFWriteGroupAttribute(this_fname_, "Header", "NumPart_Total_HighWord", from_6array<unsigned>(npartTotalHighWord_)); HDFWriteGroupAttribute(this_fname_, "Header", "NumPart_Total_HighWord", from_7array<unsigned>(npartTotalHighWord_));
HDFWriteGroupAttribute(this_fname_, "Header", "NumPart_ThisFile", from_6array<unsigned>(npart_)); HDFWriteGroupAttribute(this_fname_, "Header", "NumPart_ThisFile", from_7array<unsigned>(npart_));
HDFWriteGroupAttribute(this_fname_, "Header", "MassTable", from_6array<double>(mass_)); HDFWriteGroupAttribute(this_fname_, "Header", "MassTable", from_7array<double>(mass_));
HDFWriteGroupAttribute(this_fname_, "Header", "Time", from_value<double>(time_)); HDFWriteGroupAttribute(this_fname_, "Header", "Time", from_value<double>(time_));
HDFWriteGroupAttribute(this_fname_, "Header", "Flag_Entropy_ICs", from_value<int>(0)); HDFWriteGroupAttribute(this_fname_, "Header", "Flag_Entropy_ICs", from_value<int>(0));