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()