From 7da0bf6193d7828d4d04d98ce503664d3318bf25 Mon Sep 17 00:00:00 2001 From: David Minton Date: Mon, 10 May 2021 12:29:15 -0400 Subject: [PATCH] Reduced fragment numbers back down to avoid blowing up the system with fragments --- src/symba/symba_casedisruption.f90 | 2 +- src/symba/symba_casehitandrun.f90 | 2 +- src/symba/symba_casesupercatastrophic.f90 | 2 +- 3 files changed, 3 insertions(+), 3 deletions(-) diff --git a/src/symba/symba_casedisruption.f90 b/src/symba/symba_casedisruption.f90 index 419c9f7bb..0559e1ed2 100644 --- a/src/symba/symba_casedisruption.f90 +++ b/src/symba/symba_casedisruption.f90 @@ -32,7 +32,7 @@ function symba_casedisruption (symba_plA, family, nmergeadd, mergeadd_list, x, v logical :: lmerge ! Collisional fragments will be uniformly distributed around the pre-impact barycenter - nfrag = 22 ! This value is set for testing. This needs to be updated such that it is calculated or set by the user + nfrag = 11 ! This value is set for testing. This needs to be updated such that it is calculated or set by the user allocate(m_frag(nfrag)) allocate(rad_frag(nfrag)) allocate(xb_frag(NDIM, nfrag)) diff --git a/src/symba/symba_casehitandrun.f90 b/src/symba/symba_casehitandrun.f90 index 30e300817..08f0682cb 100644 --- a/src/symba/symba_casehitandrun.f90 +++ b/src/symba/symba_casehitandrun.f90 @@ -51,7 +51,7 @@ function symba_casehitandrun(symba_plA, family, nmergeadd, mergeadd_list, id, x, nfrag = 0 lpure = .true. else ! Imperfect hit and run, so we'll keep the largest body and destroy the other - nfrag = 12 + nfrag = 11 lpure = .false. allocate(m_frag(nfrag)) allocate(id_frag(nfrag)) diff --git a/src/symba/symba_casesupercatastrophic.f90 b/src/symba/symba_casesupercatastrophic.f90 index f24ac637b..e14d8564e 100644 --- a/src/symba/symba_casesupercatastrophic.f90 +++ b/src/symba/symba_casesupercatastrophic.f90 @@ -32,7 +32,7 @@ function symba_casesupercatastrophic (symba_plA, family, nmergeadd, mergeadd_lis logical :: lmerge ! Collisional fragments will be uniformly distributed around the pre-impact barycenter - nfrag = 42 ! This value is set for testing. This needs to be updated such that it is calculated or set by the user + nfrag = 20 ! This value is set for testing. This needs to be updated such that it is calculated or set by the user allocate(m_frag(nfrag)) allocate(rad_frag(nfrag)) allocate(xb_frag(NDIM, nfrag))