Merge branch 'debug'

MintonGroup · Aug 23, 2021 · 24458dd · 24458dd
2 parents be2d23b + 01597dd
commit 24458dd
Show file tree

Hide file tree

Showing 3 changed files with 467 additions and 467 deletions.
diff --git a/examples/symba_chambers_2013/init_cond.py b/examples/symba_chambers_2013/init_cond.py
@@ -30,6 +30,7 @@
 sim.param['ENERGY_OUT'] = "energy.dat"
 sim.param['PL_IN'] = "pl_chambers_2013.in"
 sim.param['CB_IN'] = "sun_MsunAUYR.in"
+sim.param['ENC_OUT'] = ""
 
 # Simulation parameters
 sim.param['FRAGMENTATION'] = "YES"
@@ -59,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)

diff --git a/examples/symba_chambers_2013/param.in b/examples/symba_chambers_2013/param.in
@@ -22,7 +22,6 @@ MU2KG            1.988409870698051e+30
 TU2S             31557600.0
 DU2M             149597870700.0
 IN_FORM          EL
-ENC_OUT          enc.dat
 EXTRA_FORCE      NO
 DISCARD_OUT      discard.out
 PARTICLE_OUT     particle.dat