From af220ab388aeb3d0c643ecad5f2352ccbfab95ae Mon Sep 17 00:00:00 2001 From: David Minton Date: Mon, 23 Aug 2021 18:01:22 -0400 Subject: [PATCH] Updated initial conditions distribution to better match what Chambers did --- examples/symba_chambers_2013/init_cond.py | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/examples/symba_chambers_2013/init_cond.py b/examples/symba_chambers_2013/init_cond.py index 6a6dd8526..d02830a11 100755 --- a/examples/symba_chambers_2013/init_cond.py +++ b/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 = default_rng().uniform(0.0, 0.01, Nb) -evals = default_rng().uniform(0.0, 0.01, Ns) -incvalb = default_rng().uniform(0.0, 0.005 * 180 / np.pi, Nb) -incvals = default_rng().uniform(0.0, 0.005 * 180 / np.pi, 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) 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)