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

Commit

Permalink
Added Chambers 2013 initial conditions to examples folder
Browse files Browse the repository at this point in the history
  • Loading branch information
daminton committed May 28, 2021
1 parent 5aa3dd2 commit 51b7f69
Show file tree
Hide file tree
Showing 6 changed files with 985 additions and 0 deletions.
3 changes: 3 additions & 0 deletions examples/symba_chambers2013/.gitignore
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
!*.in
!*.sh
!*.ipynb
942 changes: 942 additions & 0 deletions examples/symba_chambers2013/chambers_2013.in

Large diffs are not rendered by default.

2 changes: 2 additions & 0 deletions examples/symba_chambers2013/enc.in
Original file line number Diff line number Diff line change
@@ -0,0 +1,2 @@
param.in
10
36 changes: 36 additions & 0 deletions examples/symba_chambers2013/param.in
Original file line number Diff line number Diff line change
@@ -0,0 +1,36 @@
!
! Parameter file for Chambers 2013 in units of Msun, AU, year
!
T0 0.0e0
TSTOP 1e8 ! simulation length in years
DT 0.016 ! stepsize in years
PL_IN chambers_2013.in
TP_IN tp.in
IN_TYPE ASCII
ISTEP_OUT 6250 ! output cadence
ISTEP_DUMP 6250 ! system dump cadence
BIN_OUT bin.dat
PARTICLE_FILE particle.dat
OUT_TYPE REAL8 ! double precision real output
OUT_FORM EL ! osculating element output
OUT_STAT REPLACE
J2 0.0 ! no J2 term
J4 0.0 ! no J4 term
CHK_CLOSE yes ! check for planetary close encounters
CHK_RMIN 0.00465047 ! check for close solar encounters in AU
CHK_RMAX 100000.0 ! discard outside of
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.dat
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
FRAGMENTATION yes
MU2KG 1.98847e30 ! (M_sun-> kg)
DU2M 1.495979e11 ! distance unit to meters (AU --> m)
TU2S 3.15569259747e7 ! time unit to seconds (years --> seconds)
MTINY 1e-10 ! (G*M)
ENERGY yes
ROTATION yes
1 change: 1 addition & 0 deletions examples/symba_chambers2013/start.in
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
param.in
1 change: 1 addition & 0 deletions examples/symba_chambers2013/tp.in
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
0

0 comments on commit 51b7f69

Please sign in to comment.