Overhauled crater formation routine. All raised rims are now part of …

…the ejecta blanket
MintonGroup · Feb 9, 2017 · 6347889 · 6347889
1 parent 2cf4bd6
commit 6347889
Showing 1 changed file with 1 addition and 3 deletions.
diff --git a/src/ejecta/ejecta_distance_estimate.f90 b/src/ejecta/ejecta_distance_estimate.f90
@@ -61,9 +61,7 @@ subroutine ejecta_distance_estimate(user,crater,domain,ejdis_estimate)
    ! Extrapolate in logspace out to the minimum thickness
    slope = log(lrad2 / lrad1) / log(thick2 / thick1)
    ejdis_estimate = exp(slope * log(domain%small / thick1) + log(lrad1))
-   if (.not.user%discontinuous) then
-      ejdis_estimate = min(ejdis_estimate, 2.3_DP * crater%frad**(1.006_DP))
-   end if
+   ejdis_estimate = min(ejdis_estimate, crater%fcrat * user%ejecta_truncation)
 
    return
 end subroutine ejecta_distance_estimate