From 0c83693178d988d52a86a0a165ace4f1528095e2 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Fri, 10 Sep 2021 10:35:15 -0400 Subject: [PATCH] Fixed collision outcome initial minimum mass values to prevent runaway collisional cascade --- examples/symba_energy_momentum/param.disruption_headon.in | 8 ++++---- .../symba_energy_momentum/param.disruption_off_axis.in | 2 +- examples/symba_energy_momentum/param.hitandrun.in | 3 ++- .../param.supercatastrophic_headon.in | 2 +- .../param.supercatastrophic_off_axis.in | 2 +- 5 files changed, 9 insertions(+), 8 deletions(-) diff --git a/examples/symba_energy_momentum/param.disruption_headon.in b/examples/symba_energy_momentum/param.disruption_headon.in index 24e79be5d..a012c5505 100644 --- a/examples/symba_energy_momentum/param.disruption_headon.in +++ b/examples/symba_energy_momentum/param.disruption_headon.in @@ -1,6 +1,6 @@ T0 0.0e0 -TSTOP 0.000100 ! simulation length in seconds = 100 years -DT 0.0000001 ! stepsize in seconds +TSTOP 0.0001000 ! simulation length in seconds = 100 years +DT 0.0000001 ! stepsize in seconds CB_IN cb.in PL_IN disruption_headon.in TP_IN tp.in @@ -21,8 +21,8 @@ ENC_OUT enc.disruption_headon.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 -GMTINY 1.0e-16 -MIN_GMFRAG 1.0e-10 +GMTINY 1.0e-11 +MIN_GMFRAG 1.0e-11 FRAGMENTATION yes MU2KG 1.98908e30 TU2S 3.1556925e7 diff --git a/examples/symba_energy_momentum/param.disruption_off_axis.in b/examples/symba_energy_momentum/param.disruption_off_axis.in index b10ecee63..831f310d5 100644 --- a/examples/symba_energy_momentum/param.disruption_off_axis.in +++ b/examples/symba_energy_momentum/param.disruption_off_axis.in @@ -1,5 +1,5 @@ T0 0.0e0 -TSTOP 0.0000100 ! simulation length in seconds = 100 years +TSTOP 0.0001000 ! simulation length in seconds = 100 years DT 0.0000001 ! stepsize in seconds CB_IN cb.in PL_IN disruption_off_axis.in diff --git a/examples/symba_energy_momentum/param.hitandrun.in b/examples/symba_energy_momentum/param.hitandrun.in index c4c29a05b..1e2e2eea7 100644 --- a/examples/symba_energy_momentum/param.hitandrun.in +++ b/examples/symba_energy_momentum/param.hitandrun.in @@ -22,7 +22,8 @@ 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 +GMTINY 1.0e-11 +MIN_GMFRAG 1.0e-11 FRAGMENTATION yes MU2KG 1.98908e30 TU2S 3.1556925e7 diff --git a/examples/symba_energy_momentum/param.supercatastrophic_headon.in b/examples/symba_energy_momentum/param.supercatastrophic_headon.in index 1effdd15b..fd15c00dd 100644 --- a/examples/symba_energy_momentum/param.supercatastrophic_headon.in +++ b/examples/symba_energy_momentum/param.supercatastrophic_headon.in @@ -20,7 +20,7 @@ CHK_QMIN -1.0 ! ignore this check 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 +GMTINY 1.0e-11 MIN_GMFRAG 1.0e-11 FRAGMENTATION yes MU2KG 1.98908e30 diff --git a/examples/symba_energy_momentum/param.supercatastrophic_off_axis.in b/examples/symba_energy_momentum/param.supercatastrophic_off_axis.in index 47bce3775..51ab91e5d 100644 --- a/examples/symba_energy_momentum/param.supercatastrophic_off_axis.in +++ b/examples/symba_energy_momentum/param.supercatastrophic_off_axis.in @@ -21,7 +21,7 @@ ENC_OUT enc.supercatastrophic_off_axis.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 -GMTINY 1.0e-16 +GMTINY 1.0e-11 MIN_GMFRAG 1.0e-11 FRAGMENTATION yes MU2KG 1.98908e30