Undid some damage due to previous refactoring as this subroutine actu…

…ally wants an mtiny value not Gmtiny
MintonGroup · Aug 6, 2021 · 8d5455c · 8d5455c
1 parent 9214b2a
commit 8d5455c
Showing 1 changed file with 2 additions and 3 deletions.
diff --git a/src/fragmentation/fragmentation.f90 b/src/fragmentation/fragmentation.f90
@@ -835,8 +835,7 @@ end subroutine restructure_failed_fragments
    end subroutine fragmentation_initialize
 
 
-
-   module subroutine fragmentation_regime(Mcb, m1, m2, rad1, rad2, xh1, xh2, vb1, vb2, den1, den2, regime, Mlr, Mslr, Gmtiny, Qloss)
+   module subroutine fragmentation_regime(Mcb, m1, m2, rad1, rad2, xh1, xh2, vb1, vb2, den1, den2, regime, Mlr, Mslr, mtiny, Qloss)
       !! Author: Jennifer L.L. Pouplin, Carlisle A. Wishard, and David A. Minton
       !!
       !! Determine the collisional regime of two colliding bodies. 
@@ -931,7 +930,7 @@ module subroutine fragmentation_regime(Mcb, m1, m2, rad1, rad2, xh1, xh2, vb1, v
       Qloss = 0.0_DP
       U_binding = (3.0_DP * Mtot) / (5.0_DP * Rp) ! LS12 eq. 27
 
-      if ((m1 < Gmtiny).or.(m2 < Gmtiny)) then 
+      if ((m1 < mtiny).or.(m2 < mtiny)) then 
          regime = COLLRESOLVE_REGIME_MERGE !perfect merging regime
          Mlr = Mtot
          Mslr = 0.0_DP