diff --git a/src/plugins/output_swift.cc b/src/plugins/output_swift.cc index e66ee2c..8008155 100644 --- a/src/plugins/output_swift.cc +++ b/src/plugins/output_swift.cc @@ -152,7 +152,11 @@ public: std::vector data( npart_[0], ceint ); HDFWriteDataset(this_fname_, "PartType0/InternalEnergy", data); - data.assign( npart_[0], boxsize_ / cf_.get_value("setup","GridRes") ); + const double h = boxsize_ / hubble_param_ / cf_.get_value("setup","GridRes"); + + music::ilog.Print("Swift : set initial smoothing length to mean inter-part separation: %.2f Mpc", h); + + data.assign( npart_[0], h); HDFWriteDataset(this_fname_, "PartType0/SmoothingLength", data); }