enforceDensityFloor on setInitialConditions

Author: chongchonghe

src/problems/StarCluster/testStarCluster.cpp

@@ -259,9 +259,6 @@ auto problem_main() -> int
 	// initialize
 	sim.setInitialConditions();
 
-	// evolve
-	sim.evolve();
-
 	// read output variables
 	auto [position, values] = fextract(sim.state_new_cc_[0], sim.Geom(0), 2, 0.0, true); // z direction
 	const int nz = static_cast<int>(position.size());
@@ -312,6 +309,9 @@ auto problem_main() -> int
 	matplotlibcpp::save("./star_cluster_density_floor_z.pdf");
 #endif // HAVE_PYTHON
 
+	// evolve
+	sim.evolve();
+
 	amrex::Print() << "Finished." << '\n';
 	return status;
 }

src/simulation.hpp

@@ -972,6 +972,8 @@ template <typename problem_t> void AMRSimulation<problem_t>::setInitialCondition
 		ReadCheckpointFile();
 	}
 
+	enforceDensityFloor(0, state_new_cc_[0]); // enforce density floor at level 0 on the initial conditions
+
 	// Ensure consistency between particle radiation settings and luminosity data table configuration
 	if constexpr (Physics_Traits<problem_t>::is_radiation_enabled) {
 		if (particleRegister_.HasRadiatingParticles()) {