Simplified the Fraggle position generation. The two largest bodies ar…

…e always placed approximately where the original two bodies (or their equivalent) were.
MintonGroup · Jan 31, 2023 · 1596c2c · 1596c2c
1 parent c43cbd2
commit 1596c2c
Showing 1 changed file with 5 additions and 12 deletions.
diff --git a/src/fraggle/fraggle_generate.f90 b/src/fraggle/fraggle_generate.f90
@@ -311,21 +311,14 @@ module subroutine fraggle_generate_pos_vec(collider, nbody_system, param)
                fragment_cloud_radius(:) = impactors%radius(:)
                fragment_cloud_center(:,1) = impactors%rc(:,1) 
                fragment_cloud_center(:,2) = impactors%rc(:,2) + rdistance * impactors%bounce_unit(:)
-            else if (lsupercat) then
-               fragment_cloud_center(:,1) = impactors%rc(:,1) 
-               fragment_cloud_center(:,2) = impactors%rc(:,2) 
-               fragment_cloud_radius(:) = cloud_size_scale_factor * rdistance 
             else
-               fragment_cloud_center(:,1) = impactors%rbimp(:) - impactors%radius(1) * impactors%y_unit(:)
-               fragment_cloud_center(:,2) = impactors%rbimp(:) + impactors%radius(2) * impactors%y_unit(:)
+               fragment_cloud_center(:,1) = impactors%rbimp(:) - 1.001_DP * max(fragments%radius(1),impactors%radius(1)) * impactors%y_unit(:)
+               fragment_cloud_center(:,2) = impactors%rbimp(:) + 1.001_DP * max(fragments%radius(2),impactors%radius(2)) * impactors%y_unit(:)
                fragment_cloud_radius(:) = cloud_size_scale_factor * rdistance
             end if
-            if (lsupercat) then
-               istart = 1
-            else
-               fragments%rc(:,1) = fragment_cloud_center(:,1)
-               istart = 2
-            end if
+            fragments%rc(:,1) = fragment_cloud_center(:,1)
+            fragments%rc(:,2) = fragment_cloud_center(:,2)
+            istart = 3
 
             do i = istart, nfrag
                if (loverlap(i)) then