From 6347889801ac9642e15270647d9b513c39d9c9dd Mon Sep 17 00:00:00 2001 From: daminton Date: Thu, 9 Feb 2017 18:56:40 +0000 Subject: [PATCH] Overhauled crater formation routine. All raised rims are now part of the ejecta blanket --- src/ejecta/ejecta_distance_estimate.f90 | 4 +--- 1 file changed, 1 insertion(+), 3 deletions(-) diff --git a/src/ejecta/ejecta_distance_estimate.f90 b/src/ejecta/ejecta_distance_estimate.f90 index 5e2794bf..b04f2807 100644 --- 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