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

Commit

Permalink
Merge branch 'OOPSymba' into OOPTides
Browse files Browse the repository at this point in the history
  • Loading branch information
daminton committed Jul 30, 2021
2 parents 3efff13 + 479cb33 commit dd3a63e
Show file tree
Hide file tree
Showing 87 changed files with 2,845 additions and 1,602 deletions.
2 changes: 1 addition & 1 deletion Makefile.Defines
Original file line number Diff line number Diff line change
Expand Up @@ -71,7 +71,7 @@ FORTRAN = ifort
#AR = xiar

#FORTRAN = gfortran
#FFLAGS = -ffree-line-length-none $(GDEBUG) #$(GMEM)
#FFLAGS = -ffree-line-length-none $(GDEBUG) $(GMEM)
AR = ar

# DO NOT include in CFLAGS the "-c" option to compile object only
Expand Down
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
0
0.00029591220819207774
39.476926408897626
0.004650467260962157
4.7535806948127355e-12
-2.2473967953572827e-18
51 changes: 51 additions & 0 deletions examples/helio_gr_test/init_cond.py
Original file line number Diff line number Diff line change
@@ -0,0 +1,51 @@
#!/usr/bin/env python3
import swiftest

sim = swiftest.Simulation()
sim.param['PL_IN'] = "pl.swiftest.in"
sim.param['TP_IN'] = "tp.swiftest.in"
sim.param['CB_IN'] = "cb.swiftest.in"
sim.param['BIN_OUT'] = "bin.swiftest.dat"
sim.param['ENC_OUT'] = "enc.swiftest.dat"

sim.param['MU2KG'] = swiftest.MSun
sim.param['TU2S'] = swiftest.YR2S
sim.param['DU2M'] = swiftest.AU2M
sim.param['T0'] = 0.0
sim.param['DT'] = 0.25 * swiftest.JD2S / swiftest.YR2S
sim.param['TSTOP'] = 1000.0
sim.param['ISTEP_OUT'] = 1461
sim.param['ISTEP_DUMP'] = 1461
sim.param['CHK_QMIN_COORD'] = "HELIO"
sim.param['CHK_QMIN'] = swiftest.RSun / swiftest.AU2M
sim.param['CHK_QMIN_RANGE'] = f"{swiftest.RSun / swiftest.AU2M} 1000.0"
sim.param['CHK_RMIN'] = swiftest.RSun / swiftest.AU2M
sim.param['CHK_RMAX'] = 1000.0
sim.param['CHK_EJECT'] = 1000.0
sim.param['OUT_FORM'] = "EL"
sim.param['OUT_STAT'] = "UNKNOWN"
sim.param['GR'] = 'YES'

bodyid = {
"Sun": 0,
"Mercury": 1,
"Venus": 2,
"Earth": 3,
"Mars": 4,
"Jupiter": 5,
"Saturn": 6,
"Uranus": 7,
"Neptune": 8,
}

for name, id in bodyid.items():
sim.add(name, idval=id)

sim.save("param.swiftest.in")
sim.param['PL_IN'] = "pl.swifter.in"
sim.param['TP_IN'] = "tp.swifter.in"
sim.param['BIN_OUT'] = "bin.swifter.dat"
sim.param['ENC_OUT'] = "enc.swifter.dat"
sim.save("param.swifter.in", codename="Swifter")


27 changes: 27 additions & 0 deletions examples/helio_gr_test/param.swifter.in
Original file line number Diff line number Diff line change
@@ -0,0 +1,27 @@
! VERSION Swifter parameter file converted from Swiftest
T0 0.0
TSTOP 1000.0
DT 0.0006844626967830253
ISTEP_OUT 1461
ISTEP_DUMP 1461
OUT_FORM EL
OUT_TYPE REAL8
OUT_STAT UNKNOWN
IN_TYPE ASCII
PL_IN pl.swifter.in
TP_IN tp.swifter.in
BIN_OUT bin.swifter.dat
ENC_OUT enc.swifter.dat
CHK_QMIN 0.004650467260962157
CHK_RMIN 0.004650467260962157
CHK_RMAX 1000.0
CHK_EJECT 1000.0
CHK_QMIN_COORD HELIO
CHK_QMIN_RANGE 0.004650467260962157 1000.0
EXTRA_FORCE NO
BIG_DISCARD NO
CHK_CLOSE YES
C 63241.07708426628
J2 4.7535806948127355e-12
J4 -2.2473967953572827e-18
RHILL_PRESENT YES
35 changes: 35 additions & 0 deletions examples/helio_gr_test/param.swiftest.in
Original file line number Diff line number Diff line change
@@ -0,0 +1,35 @@
! VERSION Swiftest parameter input
T0 0.0
TSTOP 1000.0
DT 0.0006844626967830253
ISTEP_OUT 1461
ISTEP_DUMP 1461
OUT_FORM EL
OUT_TYPE REAL8
OUT_STAT UNKNOWN
IN_TYPE ASCII
PL_IN pl.swiftest.in
TP_IN tp.swiftest.in
CB_IN cb.swiftest.in
BIN_OUT bin.swiftest.dat
ENC_OUT enc.swiftest.dat
CHK_QMIN 0.004650467260962157
CHK_RMIN 0.004650467260962157
CHK_RMAX 1000.0
CHK_EJECT 1000.0
CHK_QMIN_COORD HELIO
CHK_QMIN_RANGE 0.004650467260962157 1000.0
MU2KG 1.988409870698051e+30
TU2S 31557600.0
DU2M 149597870700.0
EXTRA_FORCE NO
BIG_DISCARD NO
CHK_CLOSE YES
FRAGMENTATION NO
ROTATION NO
TIDES NO
ENERGY NO
GR YES
YARKOVSKY NO
YORP NO
MTINY 0.0
36 changes: 36 additions & 0 deletions examples/helio_gr_test/pl.swifter.in
Original file line number Diff line number Diff line change
@@ -0,0 +1,36 @@
9
0 39.476926408897625196
0.0 0.0 0.0
0.0 0.0 0.0
1 6.5537098095653139645e-06 0.0014751234419554511911
1.6306381826061645943e-05
0.13267502226188271353 0.2786606257975073886 0.010601098875389479426
-11.331978934667442676 4.8184460126705647045 1.4332264599878684131
2 9.663313399581537916e-05 0.00675908960945781479
4.0453784346544178454e-05
-0.69398700025820403425 -0.19235393648106968723 0.03740673057980103272
1.9245789988923785786 -7.1528261190002948057 -0.20922405362759749996
3 0.000120026935827952453094 0.010044837538502923644
4.25875607065040958e-05
0.49463573470256239073 -0.8874896493821613497 4.051630875713834232e-05
5.386704768180099809 3.0357508899436080915 -0.00016218409216515533796
4 1.2739802010675941456e-05 0.0072467236860282326973
2.265740805092889601e-05
-1.5655322071100350456 0.56626121192188216824 0.050269397991054412533
-1.5477080637857006753 -4.370087697214287981 -0.05361768768801557225
5 0.037692251088985676735 0.35527094075555771578
0.00046732617030490929307
4.0891378954287338487 -2.9329188614380639066 -0.07930573161132697946
1.575024788882753283 2.3719591091996699917 -0.045089307261129988257
6 0.011285899820091272997 0.43765464106459166412
0.00038925687730393611812
6.3349788609660162564 -7.674600716671800882 -0.11868650931385750502
1.4598618704191345578 1.2948691245181617393 -0.080593167691228835176
7 0.0017236589478267730203 0.46956055286931676728
0.00016953449859497231466
14.832516206189200858 13.032608531076540714 -0.14378102535616668622
-0.9573374666934839659 1.014553546383260322 0.016118112341773867214
8 0.0020336100526728302319 0.7813163071687303693
0.000164587904124493665
29.561664938083289655 -4.6012285192418387325 -0.586585578731106283
0.17051705220469790965 1.1424784769020628332 -0.027423757798549895085
33 changes: 33 additions & 0 deletions examples/helio_gr_test/pl.swiftest.in
Original file line number Diff line number Diff line change
@@ -0,0 +1,33 @@
8
1 6.5537098095653139645e-06
1.6306381826061645943e-05
0.13267502226188271353 0.2786606257975073886 0.010601098875389479426
-11.331978934667442676 4.8184460126705647045 1.4332264599878684131
2 9.663313399581537916e-05
4.0453784346544178454e-05
-0.69398700025820403425 -0.19235393648106968723 0.03740673057980103272
1.9245789988923785786 -7.1528261190002948057 -0.20922405362759749996
3 0.000120026935827952453094
4.25875607065040958e-05
0.49463573470256239073 -0.8874896493821613497 4.051630875713834232e-05
5.386704768180099809 3.0357508899436080915 -0.00016218409216515533796
4 1.2739802010675941456e-05
2.265740805092889601e-05
-1.5655322071100350456 0.56626121192188216824 0.050269397991054412533
-1.5477080637857006753 -4.370087697214287981 -0.05361768768801557225
5 0.037692251088985676735
0.00046732617030490929307
4.0891378954287338487 -2.9329188614380639066 -0.07930573161132697946
1.575024788882753283 2.3719591091996699917 -0.045089307261129988257
6 0.011285899820091272997
0.00038925687730393611812
6.3349788609660162564 -7.674600716671800882 -0.11868650931385750502
1.4598618704191345578 1.2948691245181617393 -0.080593167691228835176
7 0.0017236589478267730203
0.00016953449859497231466
14.832516206189200858 13.032608531076540714 -0.14378102535616668622
-0.9573374666934839659 1.014553546383260322 0.016118112341773867214
8 0.0020336100526728302319
0.000164587904124493665
29.561664938083289655 -4.6012285192418387325 -0.586585578731106283
0.17051705220469790965 1.1424784769020628332 -0.027423757798549895085
192 changes: 192 additions & 0 deletions examples/helio_gr_test/swiftest_relativity.ipynb

Large diffs are not rendered by default.

1 change: 1 addition & 0 deletions examples/helio_gr_test/tp.swifter.in
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
0
1 change: 1 addition & 0 deletions examples/helio_gr_test/tp.swiftest.in
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
0
Original file line number Diff line number Diff line change
@@ -0,0 +1,5 @@
0
0.00029591220819207774
0.004650467260962157
0.0
0.0
Original file line number Diff line number Diff line change
Expand Up @@ -25,8 +25,8 @@
sim.param['T0'] = 0.0
sim.param['DT'] = 1.0
sim.param['TSTOP'] = 365.25e1
sim.param['ISTEP_OUT'] = 11
sim.param['ISTEP_DUMP'] = 1
sim.param['ISTEP_OUT'] = 10
sim.param['ISTEP_DUMP'] = 10
sim.param['CHK_QMIN_COORD'] = "HELIO"
sim.param['CHK_QMIN'] = swiftest.RSun / swiftest.AU2M
sim.param['CHK_QMIN_RANGE'] = f"{swiftest.RSun / swiftest.AU2M} 1000.0"
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -2,8 +2,8 @@
T0 0.0
TSTOP 3652.5
DT 1.0
ISTEP_OUT 11
ISTEP_DUMP 1
ISTEP_OUT 10
ISTEP_DUMP 10
OUT_FORM XV
OUT_TYPE REAL8
OUT_STAT UNKNOWN
Expand All @@ -21,6 +21,6 @@ CHK_QMIN_RANGE 0.004650467260962157 1000.0
EXTRA_FORCE NO
BIG_DISCARD NO
CHK_CLOSE YES
J2 4.7535806948127355e-12
J4 -2.2473967953572827e-18
RHILL_PRESENT YES
J2 0.0
J4 0.0
Original file line number Diff line number Diff line change
Expand Up @@ -2,8 +2,8 @@
T0 0.0
TSTOP 3652.5
DT 1.0
ISTEP_OUT 11
ISTEP_DUMP 1
ISTEP_OUT 10
ISTEP_DUMP 10
OUT_FORM XV
OUT_TYPE REAL8
OUT_STAT UNKNOWN
Expand All @@ -25,6 +25,7 @@ DU2M 149597870700.0
EXTRA_FORCE NO
BIG_DISCARD NO
CHK_CLOSE YES
RHILL_PRESENT YES
FRAGMENTATION NO
ROTATION NO
TIDES NO
Expand Down
33 changes: 33 additions & 0 deletions examples/rmvs_swifter_comparison/8pl_16tp_encounters/pl.in
Original file line number Diff line number Diff line change
@@ -0,0 +1,33 @@
8
1 4.9125474498983623693e-11 0.0014751239400086721089
1.6306381826061645943e-05
-0.09861361766419070307 0.29750596935836171042 0.03335708456145129036
-0.032353632540864457612 -0.0078122718370876240157 0.0023293874953380202045
2 7.243452483873646905e-10 0.0067590794275223005208
4.0453784346544178454e-05
-0.6439817957564198947 -0.3248550380869373866 0.032702713983447248558
0.008969709495375973937 -0.018153139924556138673 -0.0007667345025597138231
3 8.9970113821660187435e-10 0.010044873080337524463
4.25875607065040958e-05
0.59421674333603324847 -0.82331253628773626296 3.7129329104855261984e-05
0.013670550280388280365 0.010004295439859960809 -5.226292361234363611e-07
4 9.549535102761465607e-11 0.0072467054748629370034
2.265740805092889601e-05
-1.592721551706784977 0.48166390206865000723 0.049163460846716633412
-0.0035287723306552309585 -0.01219974682608557931 -0.00016910795626524249315
5 2.825345908631354893e-07 0.35527074967975702942
0.00046732617030490929307
4.119089774477131094 -2.8872942462256898644 -0.080165336328135106125
0.004245402942744468111 0.0065414198811065849687 -0.00012215100047356211078
6 8.459715183006415395e-08 0.4376562090257202473
0.00038925687730393611812
6.3629100567525149756 -7.649727796147929304 -0.12023019299387090186
0.0039834472120812329868 0.0035613826786502411278 -0.00022039988214595340028
7 1.2920249163736673626e-08 0.4695793205674148502
0.00016953449859497231466
14.814154683311180349 13.052040295401360126 -0.14347198499748289868
-0.002625101393275708784 0.0027742356008832688187 4.416821810149910185e-05
8 1.5243589003230834323e-08 0.7813388398513013378
0.000164587904124493665
29.564924658285640646 -4.579331535234244299 -0.5871109926822926095
0.00046449847307956888343 0.003128345390031967918 -7.5036135696161668576e-05
36 changes: 36 additions & 0 deletions examples/rmvs_swifter_comparison/8pl_16tp_encounters/pl.swifter.in
Original file line number Diff line number Diff line change
@@ -0,0 +1,36 @@
9
0 0.00029591220819207775568
0.0 0.0 0.0
0.0 0.0 0.0
1 4.9125474498983623693e-11 0.0014751239400086721089
1.6306381826061645943e-05
-0.09861361766419070307 0.29750596935836171042 0.03335708456145129036
-0.032353632540864457612 -0.0078122718370876240157 0.0023293874953380202045
2 7.243452483873646905e-10 0.0067590794275223005208
4.0453784346544178454e-05
-0.6439817957564198947 -0.3248550380869373866 0.032702713983447248558
0.008969709495375973937 -0.018153139924556138673 -0.0007667345025597138231
3 8.9970113821660187435e-10 0.010044873080337524463
4.25875607065040958e-05
0.59421674333603324847 -0.82331253628773626296 3.7129329104855261984e-05
0.013670550280388280365 0.010004295439859960809 -5.226292361234363611e-07
4 9.549535102761465607e-11 0.0072467054748629370034
2.265740805092889601e-05
-1.592721551706784977 0.48166390206865000723 0.049163460846716633412
-0.0035287723306552309585 -0.01219974682608557931 -0.00016910795626524249315
5 2.825345908631354893e-07 0.35527074967975702942
0.00046732617030490929307
4.119089774477131094 -2.8872942462256898644 -0.080165336328135106125
0.004245402942744468111 0.0065414198811065849687 -0.00012215100047356211078
6 8.459715183006415395e-08 0.4376562090257202473
0.00038925687730393611812
6.3629100567525149756 -7.649727796147929304 -0.12023019299387090186
0.0039834472120812329868 0.0035613826786502411278 -0.00022039988214595340028
7 1.2920249163736673626e-08 0.4695793205674148502
0.00016953449859497231466
14.814154683311180349 13.052040295401360126 -0.14347198499748289868
-0.002625101393275708784 0.0027742356008832688187 4.416821810149910185e-05
8 1.5243589003230834323e-08 0.7813388398513013378
0.000164587904124493665
29.564924658285640646 -4.579331535234244299 -0.5871109926822926095
0.00046449847307956888343 0.003128345390031967918 -7.5036135696161668576e-05
Original file line number Diff line number Diff line change
@@ -0,0 +1,33 @@
8
1 4.9125474498983623693e-11 0.0014751239400086721089
1.6306381826061645943e-05
-0.09861361766419070307 0.29750596935836171042 0.03335708456145129036
-0.032353632540864457612 -0.0078122718370876240157 0.0023293874953380202045
2 7.243452483873646905e-10 0.0067590794275223005208
4.0453784346544178454e-05
-0.6439817957564198947 -0.3248550380869373866 0.032702713983447248558
0.008969709495375973937 -0.018153139924556138673 -0.0007667345025597138231
3 8.9970113821660187435e-10 0.010044873080337524463
4.25875607065040958e-05
0.59421674333603324847 -0.82331253628773626296 3.7129329104855261984e-05
0.013670550280388280365 0.010004295439859960809 -5.226292361234363611e-07
4 9.549535102761465607e-11 0.0072467054748629370034
2.265740805092889601e-05
-1.592721551706784977 0.48166390206865000723 0.049163460846716633412
-0.0035287723306552309585 -0.01219974682608557931 -0.00016910795626524249315
5 2.825345908631354893e-07 0.35527074967975702942
0.00046732617030490929307
4.119089774477131094 -2.8872942462256898644 -0.080165336328135106125
0.004245402942744468111 0.0065414198811065849687 -0.00012215100047356211078
6 8.459715183006415395e-08 0.4376562090257202473
0.00038925687730393611812
6.3629100567525149756 -7.649727796147929304 -0.12023019299387090186
0.0039834472120812329868 0.0035613826786502411278 -0.00022039988214595340028
7 1.2920249163736673626e-08 0.4695793205674148502
0.00016953449859497231466
14.814154683311180349 13.052040295401360126 -0.14347198499748289868
-0.002625101393275708784 0.0027742356008832688187 4.416821810149910185e-05
8 1.5243589003230834323e-08 0.7813388398513013378
0.000164587904124493665
29.564924658285640646 -4.579331535234244299 -0.5871109926822926095
0.00046449847307956888343 0.003128345390031967918 -7.5036135696161668576e-05
Loading

0 comments on commit dd3a63e

Please sign in to comment.