add massive particles to SphericalCollapse problem

quokka-astro · Jan 18, 2024 · d484364 · d484364
1 parent 5128936
commit d484364
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 2 deletions.
diff --git a/src/SphericalCollapse/spherical_collapse.cpp b/src/SphericalCollapse/spherical_collapse.cpp
@@ -87,7 +87,14 @@ template <> void RadhydroSimulation<CollapseProblem>::setInitialConditionsOnGrid
 
 template <> void RadhydroSimulation<CollapseProblem>::createInitialParticles()
 {
-	// note: an implementation is only required if particles are used
+	// add particles at random positions in the box
+	bool generate_on_root_rank = true;
+	int iseed = 42;
+	int num_particles = 1000;
+	double particle_mass = 0.5 / static_cast<double>(num_particles);
+
+	quokka::CICParticleContainer::ParticleInitData pdata = {{particle_mass}, {0}, {0}, {0}}; // mass vx vy vz
+	CICParticles->InitRandom(num_particles, iseed, pdata, generate_on_root_rank);
 }
 
 template <> void RadhydroSimulation<CollapseProblem>::ErrorEst(int lev, amrex::TagBoxArray &tags, amrex::Real /*time*/, int /*ngrow*/)

diff --git a/tests/SphericalCollapse.in b/tests/SphericalCollapse.in
@@ -31,7 +31,8 @@ gravity.Gconst = 1.0    # gravitational constant
 
 do_reflux = 1
 do_subcycle = 0
-do_tracers = 1  # turn on tracer particles
+do_tracers = 0          # turn on tracer particles
+do_cic_particles = 1    # turns on CIC particles
 
 ascent_interval = 100
 plotfile_interval = 200