From 882231ed1f0603b059cf6a2f350068a3406b816b Mon Sep 17 00:00:00 2001 From: David A Minton Date: Mon, 21 Jun 2021 16:22:59 -0400 Subject: [PATCH] Changed solar impact to be a solar system run in units of AU-yr instead of a Mars run --- examples/symba_energy_momentum/escape.in | 8 ++-- examples/symba_energy_momentum/param.sun.in | 43 +++++++++++---------- examples/symba_energy_momentum/sun.in | 30 ++++++++------ 3 files changed, 44 insertions(+), 37 deletions(-) diff --git a/examples/symba_energy_momentum/escape.in b/examples/symba_energy_momentum/escape.in index 41a44b9a6..1cbcc07db 100644 --- a/examples/symba_energy_momentum/escape.in +++ b/examples/symba_energy_momentum/escape.in @@ -1,9 +1,9 @@ 3 1 39.47841760435743 -0.0 0.0 0.0 -0.0 0.0 0.0 -0.4 0.4 0.4 !Ip -0.0 0.0 0.0 !rot +0.0 0.0 0.0 ! x y z +0.0 0.0 0.0 ! vx vy vz +0.0 0.0 0.07 ! ip +11.2093063 -38.75937204 82.25088158 ! rot (radian / year) 2 1e-07 0.0009 7e-06 99.9 0.0 0.0 diff --git a/examples/symba_energy_momentum/param.sun.in b/examples/symba_energy_momentum/param.sun.in index 25a2562ef..0e26d328f 100644 --- a/examples/symba_energy_momentum/param.sun.in +++ b/examples/symba_energy_momentum/param.sun.in @@ -1,7 +1,7 @@ !Parameter file for the SyMBA-RINGMOONS test T0 0.0 -TSTOP 1000.0 -DT 50.0 +TSTOP 1.0 +DT 1e-3 PL_IN sun.in TP_IN tp.in IN_TYPE ASCII @@ -10,25 +10,26 @@ ISTEP_DUMP 1 BIN_OUT bin.sun.dat PARTICLE_FILE particle.sun.dat OUT_TYPE REAL8 -OUT_FORM EL +OUT_FORM XV ! osculating element output OUT_STAT REPLACE -J2 0.0 -J4 0.0 -CHK_CLOSE yes -CHK_RMIN 3389500.0 -CHK_RMAX 3389500000.0 -CHK_EJECT 3389500000.0 -CHK_QMIN 3389500.0 -CHK_QMIN_COORD HELIO -CHK_QMIN_RANGE 3389500.0 338950000000.0 -ENC_OUT /dev/null -EXTRA_FORCE no -BIG_DISCARD no -RHILL_PRESENT yes -MTINY 1000.0 -ENERGY yes +ISTEP_DUMP 1 ! system dump cadence +J2 0.0 ! no J2 term +J4 0.0 ! no J4 term +CHK_CLOSE yes ! check for planetary close encounters +CHK_RMIN 0.005 +CHK_RMAX 1e2 +CHK_EJECT -1.0 ! ignore this check +CHK_QMIN -1.0 ! ignore this check +!CHK_QMIN_COORD HELIO ! commented out here +!CHK_QMIN_RANGE 1.0 1000.0 ! commented out here +ENC_OUT enc.escape.dat +EXTRA_FORCE no ! no extra user-defined forces +BIG_DISCARD no ! output all planets if anything discarded +RHILL_PRESENT no ! Hill's sphere radii in input file +MTINY 1.0e-16 FRAGMENTATION yes +MU2KG 1.98908e30 +TU2S 3.1556925e7 +DU2M 1.49598e11 +ENERGY yes ROTATION yes -MU2KG 1.0 -DU2M 1.0 -TU2S 1.0 diff --git a/examples/symba_energy_momentum/sun.in b/examples/symba_energy_momentum/sun.in index aa5d8d077..a84e3b125 100644 --- a/examples/symba_energy_momentum/sun.in +++ b/examples/symba_energy_momentum/sun.in @@ -1,12 +1,18 @@ - 2 - 1 42838866200000.0 - 0.0 0.0 0.0 - 0.0 0.0 0.0 - 0.400000000000000 0.400000000000000 0.400000000000000 - 0.000000000000000E+000 0.000000000000000E+000 0.000000000000000E+000 - 2099 21798.5890113035 2312.73971815119 - 3002.96021708737 - -3746045.07840497 -412560.437993027 325192.799507880 - 1689.06790223103 -3033.03138547776 -586.746139832494 - 0.400000000000000 0.400000000000000 0.400000000000000 - -1.898375172053203E-003 -3.848932263142081E-004 -8.383582091179525E-004 +3 +1 39.47841760435743 +0.0 0.0 0.0 ! x y z +0.0 0.0 0.0 ! vx vy vz +0.0 0.0 0.07 ! ip +11.2093063 -38.75937204 82.25088158 ! rot (radian / year) +2 2e-08 +3e-06 +5e-2 0.0 0.0 +0.00 10.00 0.0 +0.4 0.4 0.4 !Ip +0.0 0.0 2300.0 !rot +3 2e-08 +3e-06 +1.0 0.00E-05 0.0 +0.00 6.28 0.0 +0.4 0.4 0.4 !Ip +0.0 0.0 2300.0 !rot