From cbece6cb76628758a4f6cfb109eacec038e075d2 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Mon, 20 Mar 2023 11:05:07 -0400 Subject: [PATCH] Updated the sundiver test with more massive body to test for angular momentum conservation --- examples/solar_impact/sundiver.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) mode change 100644 => 100755 examples/solar_impact/sundiver.py diff --git a/examples/solar_impact/sundiver.py b/examples/solar_impact/sundiver.py old mode 100644 new mode 100755 index 02e320f69..c279b4ec4 --- a/examples/solar_impact/sundiver.py +++ b/examples/solar_impact/sundiver.py @@ -46,9 +46,9 @@ q = 0.9 * swiftest.RSun * sim.M2DU a = 0.1 e = 1.0 - q / a -M = 1e-4 * swiftest.MEarth * sim.KG2MU +M = 2e0 * swiftest.MEarth * sim.KG2MU R = (3 * M / (4 * np.pi * density)) ** (1.0 / 3.0) -rot = 2 * sim.init_cond.sel(name="Earth")['rot'] +rot = 4 * sim.init_cond.sel(name="Earth")['rot'] sim.add_body(name="Sundiver", a=a, e=e, inc=0.0, capom=0.0, omega=0.0, capm=180.0, mass=M, radius=R, Ip=[0.4,0.4,0.4], rot=rot) sim.get_parameter()