Fixed up encounter test cases and verified that they work

MintonGroup · Oct 8, 2021 · b958476 · b958476
1 parent 70e6147
commit b958476
Show file tree

Hide file tree

Showing 16 changed files with 136,159 additions and 72 deletions.
diff --git a/examples/rmvs_swifter_comparison/8pl_16tp_encounters/param.swiftest.in b/examples/rmvs_swifter_comparison/8pl_16tp_encounters/param.swiftest.in
@@ -33,3 +33,4 @@ ROTATION         NO
 TIDES            NO
 ENERGY           NO
 GR               NO
+ENCOUNTER_CHECK SORTSWEEP
diff --git a/examples/rmvs_swifter_comparison/8pl_16tp_encounters/swiftest_rmvs_vs_swifter_rmvs.ipynb b/examples/rmvs_swifter_comparison/8pl_16tp_encounters/swiftest_rmvs_vs_swifter_rmvs.ipynb
diff --git a/examples/symba_chambers_2013/param.in b/examples/symba_chambers_2013/param.in
@@ -36,3 +36,5 @@ GR               NO
 ENERGY_OUT       energy.dat
 GMTINY           1.1053539394491336e-08
 MIN_GMFRAG       1.1053539394491336e-10
+ENCOUNTER_CHECK  ADAPTIVE
+INTERACTION_LOOPS ADAPTIVE
diff --git a/examples/symba_clement_2018/param.in b/examples/symba_clement_2018/param.in
@@ -1,9 +1,9 @@
 ! VERSION        Swiftest parameter input from file param.in
 T0               0.0
-TSTOP            100.0
+TSTOP            1.000e2
 DT               0.016
-ISTEP_OUT        625
-ISTEP_DUMP       6250
+ISTEP_OUT        100
+ISTEP_DUMP       1000
 OUT_FORM         XVEL
 OUT_TYPE         NETCDF_DOUBLE
 OUT_STAT         REPLACE
@@ -35,3 +35,5 @@ ENERGY           YES
 GR               NO
 GMTINY           2.96336237e-09
 MIN_GMFRAG       2.96336237e-11
+ENCOUNTER_CHECK  ADAPTIVE
+INTERACTION_LOOPS FLAT
diff --git a/examples/symba_mars_disk/param.in b/examples/symba_mars_disk/param.in
@@ -1,6 +1,6 @@
 !Parameter file for the SyMBA-RINGMOONS test
 T0             0.0 
-TSTOP          1.0e12 
+TSTOP          6000.00
 DT             600.0
 CB_IN          cb.in
 PL_IN          mars.in
@@ -22,7 +22,7 @@ CHK_QMIN_COORD HELIO
 CHK_QMIN_RANGE 3389500.0 338950000000.0
 EXTRA_FORCE    no
 RHILL_PRESENT  yes
-GMTINY         1.4e4
+GMTINY         1.4e5
 MIN_GMFRAG     1.0e4
 ENERGY         yes
 FRAGMENTATION  yes 
@@ -31,4 +31,5 @@ MU2KG          1.0
 DU2M           1.0
 TU2S           1.0
 SEED           2 3080983 2220830
-INTERACTION_LOOPS ADAPTIVE
+INTERACTION_LOOPS FLAT
+ENCOUNTER_CHECK ADAPTIVE
diff --git a/examples/symba_pl02k/Jul27_2k_fully.in b/examples/symba_pl02k/Jul27_2k_fully.in
diff --git a/examples/symba_pl02k/cb.in b/examples/symba_pl02k/cb.in
@@ -0,0 +1,5 @@
+Sun
+39.476926408897626
+0.005
+0.0
+0.0
diff --git a/examples/symba_pl02k/param.in b/examples/symba_pl02k/param.in
@@ -0,0 +1,35 @@
+!
+!Parameter file for the Entire Solar System for the Terrestrial Disk
+!
+T0              0.0e0                           !starting time in years
+TSTOP           0.100                      !simulation length in years
+DT              0.001                           !step size in years
+CB_IN           cb.in
+PL_IN           Jul27_2k_fully.in
+TP_IN           tp.in
+IN_TYPE         ASCII                           !format of input file
+IN_FORM         XV                              !format of input file
+ISTEP_OUT       100                       !output cadence 
+ISTEP_DUMP      100                      !keep same as ISTEP_OUT
+BIN_OUT         bin.nc                          !output file
+OUT_TYPE        NETCDF_DOUBLE                   !double precision real output
+OUT_FORM        XVEL                            !output type in XV OR EL
+OUT_STAT        REPLACE
+CHK_CLOSE       yes                             !check for planetary close encounters
+CHK_RMAX        10000.0                         !drop things from the simulation in AU
+CHK_EJECT       -1.0                            !turns this off
+CHK_QMIN        -1.0                            !turns this off
+EXTRA_FORCE     no                              !this would be user defined forces
+BIG_DISCARD     no                              !output all planets if anything is discarded
+RHILL_PRESENT   yes                             !Rhill is definted in the input file
+FRAGMENTATION   no                              !Turns fragmentation model on
+ROTATION        no
+TIDES           no
+GR              no
+MU2KG           5.038e30                        !Conversion from mass unit to grams (M_sun*G --> kg where G=4pi^2 AU^3/(y^2*M_sun))
+DU2M            1.5e11                     		!Conversion from distance unit to centimeters (AU --> m)
+TU2S            3.16e7                          !Conversion from time unit to seconds (years --> seconds)
+GMTINY           1e-20						!(G*M)
+ENERGY          no
+INTERACTION_LOOPS ADAPTIVE
+ENCOUNTER_CHECK ADAPTIVE
diff --git a/examples/symba_pl02k/runjob.sh b/examples/symba_pl02k/runjob.sh
@@ -0,0 +1,28 @@
+#!/bin/zsh -l
+#SBATCH -A standby
+#SBATCH --time=4:00:00
+#SBATCH --nodes=1
+#SBATCH --ntasks=1
+#SBATCH --ntasks-per-node=1
+#SBATCH --job-name=swiftest-2k
+#SBATCH --cpus-per-task=128
+#SBATCH --mem-per-cpu=MaxMemPerCPU
+#SBATCH -o %x.out
+#SBATCH -e %x.err
+
+cd $SLURM_SUBMIT_DIR
+filename="${SLURM_SUBMIT_DIR}/swiftest-2k-timehist.log"
+rm ${filename}
+echo "N cores,wall time(s)" > ${filename}
+
+count=1
+while [ $count -lt 129 ]
+do
+  export OMP_NUM_THREADS=${count}
+  ./swiftest_driver symba param.in > term${count}.out
+  grep Integration term${count}.out | tail -n1 | awk -v NTHREADS=$OMP_NUM_THREADS '{print NTHREADS,",",$6}' | sed 's/.$//' >> ${filename} 
+  mv interaction_timer.log interaction_timer${count}.log
+  mv encounter_check_timer.log encounter_check_timer${count}.log
+  count=`expr $count + 1`
+done
+
diff --git a/examples/symba_pl02k/tp.in b/examples/symba_pl02k/tp.in
@@ -0,0 +1 @@
+0