Commit ec104141 authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

For now, do not remove any particles in end_force().

parent 5f0e8673
......@@ -2050,9 +2050,10 @@ void runner_do_end_force(struct runner *r, struct cell *c, int timer) {
/* Finish the force loop */
hydro_end_force(p, cosmo);
// MATTHIEU: Temporary star-formation law
if (p->rho > 1.5e7 && e->step > 2) {
message("Removing particle id=%lld rho=%e", p->id, p->rho);
cell_convert_part_to_gpart(e, c, p, xp);
// message("Removing particle id=%lld rho=%e", p->id, p->rho);
// cell_convert_part_to_gpart(e, c, p, xp);
// cell_remove_part(e,c,p,xp);
}
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment