Changed back to small e and i distributions for Chambers 2013

examples/symba_chambers_2013/init_cond.py

@@ -60,10 +60,10 @@
 # Define the initial orbital elements of the big and small bodies
 avalb = default_rng().uniform(0.3, 2.0, Nb)
 avals = default_rng().uniform(0.3, 2.0, Ns)
-evalb = np.abs(default_rng().normal(loc=0.0, scale=0.13, size=Nb))
-evals = np.abs(default_rng().normal(loc=0.0, scale=0.13, size=Ns))
-incvalb = default_rng().normal(loc=0.0, scale=1.0, size=Nb)
-incvals = default_rng().normal(loc=0.0, scale=1.0, size=Ns)
+evalb = default_rng().uniform(0.0, 0.01, Nb)
+evals = default_rng().uniform(0.0, 0.01, Ns)
+incvalb = default_rng().uniform(0.0, 0.5, Nb)
+incvals = default_rng().uniform(0.0, 0.5, Ns)
 capomvalb = default_rng().uniform(0.0, 360.0, Nb)
 capomvals = default_rng().uniform(0.0, 360.0, Ns)
 omegavalb = default_rng().uniform(0.0, 360.0, Nb)