From e4fda58c04ac4d1b8881da4cc326cf9db111d241 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Sun, 12 Feb 2023 12:09:13 -0500 Subject: [PATCH] Updated the Chambers example with new fragmentation parameters --- examples/Chambers2013/init_cond.py | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/examples/Chambers2013/init_cond.py b/examples/Chambers2013/init_cond.py index 18582769e..1b2437b6c 100755 --- a/examples/Chambers2013/init_cond.py +++ b/examples/Chambers2013/init_cond.py @@ -39,10 +39,11 @@ Ns = 140 Mb = 2.8e-7 * 14 / Nb Ms = 2.8e-8 * 140 / Ns -dens = 3000.0 / (sim.MU2KG / sim.DU2M**3) +dens = 3000.0 * sim.KG2MU / sim.M2DU**3 mtiny = 1e-2 * Ms -mininum_fragment_mass = 1e-4 * Ms +minimum_fragment_mass = 1e-5 * Ms +nfrag_reduction = 10.0 rng = default_rng(seed=3031179) runname = "Chambers (2013)" @@ -128,7 +129,7 @@ def sample(r_inner, r_break, r_outer, slope1, slope2): sim.add_solar_system_body(["Sun","Jupiter","Saturn","Uranus","Neptune"]) sim.add_body(name=nameb, a=avalb, e=evalb, inc=incvalb, capom=capomvalb, omega=omegavalb, capm=capmvalb, mass=Mvalb, radius=Rvalb, rot=rotvalb, Ip=Ipvalb) sim.add_body(name=names, a=avals, e=evals, inc=incvals, capom=capomvals, omega=omegavals, capm=capmvals, mass=Mvals, radius=Rvals, rot=rotvals, Ip=Ipvals) -sim.set_parameter(mtiny=mtiny, minimum_fragment_mass=mininum_fragment_mass) +sim.set_parameter(mtiny=mtiny, minimum_fragment_mass=minimum_fragment_mass, nfrag_reduction=nfrag_reduction) sim.set_parameter(tstop=3e8, dt=6.0875/365.25, istep_out=60000, dump_cadence=10) sim.clean()