Changed solar impact to be a solar system run in units of AU-yr inste…

…ad of a Mars run
MintonGroup · Jun 21, 2021 · 882231e · 882231e
1 parent d371a72
commit 882231e
Show file tree

Hide file tree

Showing 3 changed files with 44 additions and 37 deletions.
diff --git 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
@@ -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
@@ -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