Skip to content
This repository was archived by the owner on Aug 28, 2024. It is now read-only.

Commit

Permalink
added hit and run case to symba_energy_momentum example directory
Browse files Browse the repository at this point in the history
  • Loading branch information
cwishard committed Sep 9, 2021
1 parent f428be3 commit 7d61dd4
Show file tree
Hide file tree
Showing 2 changed files with 45 additions and 0 deletions.
13 changes: 13 additions & 0 deletions examples/symba_energy_momentum/hitandrun.in
Original file line number Diff line number Diff line change
@@ -0,0 +1,13 @@
2
1 1e-07 0.0009
7e-06
1.0 -4.20E-05 0.0
0.00 6.28 0.0
0.4 0.4 0.4 !Ip
0.0 0.0 6.0e4 !rot
2 7e-10 0.0004
3.25e-06
1.0 4.20E-05 0.0
-1.50 -6.28 0.0
0.4 0.4 0.4 !Ip
0.0 0.0 1.0e5 !rot
32 changes: 32 additions & 0 deletions examples/symba_energy_momentum/param.hitandrun.in
Original file line number Diff line number Diff line change
@@ -0,0 +1,32 @@
T0 0.0e0
TSTOP 0.000100 ! simulation length in seconds = 100 years
DT 0.0000001 ! stepsize in seconds
CB_IN cb.in
PL_IN hitandrun.in
TP_IN tp.in
IN_TYPE ASCII
ISTEP_OUT 1 ! output cadence every year
BIN_OUT bin.hitandrun.nc
PARTICLE_OUT particle.hitandrun.dat
OUT_TYPE NETCDF_DOUBLE ! double precision real output
OUT_FORM XVEL ! osculating element output
OUT_STAT REPLACE
ISTEP_DUMP 1 ! system dump cadence
CHK_CLOSE yes ! check for planetary close encounters
CHK_RMIN 0.005
CHK_RMAX 1e6
CHK_EJECT -1.0 ! ignore this check
CHK_QMIN -1.0 ! ignore this check
ENC_OUT enc.hitandrun.dat
DISCARD_OUT discard.hitandrun.out
EXTRA_FORCE no ! no extra user-defined forces
BIG_DISCARD no ! output all planets if anything discarded
RHILL_PRESENT yes ! Hill's sphere radii in input file
GMTINY 1.0e-16
FRAGMENTATION yes
MU2KG 1.98908e30
TU2S 3.1556925e7
DU2M 1.49598e11
ENERGY yes
ROTATION yes
SEED 8 933097 -220886113 -118730874 233084005 32111237 -823335422 524551114 -61162322

0 comments on commit 7d61dd4

Please sign in to comment.