From 96034f76e6b46e6b1a6f8dfcc4b2e07f85f7053e Mon Sep 17 00:00:00 2001 From: daminton Date: Mon, 5 Dec 2016 21:29:29 +0000 Subject: [PATCH] Ejecta modified to be more self-consistent with empirical data on ejecta thickness (calibrated for strengthless soft rock) --- src/ejecta/ejecta_distance_estimate.f90 | 3 +++ 1 file changed, 3 insertions(+) diff --git a/src/ejecta/ejecta_distance_estimate.f90 b/src/ejecta/ejecta_distance_estimate.f90 index 65e87ce1..5e2794bf 100644 --- a/src/ejecta/ejecta_distance_estimate.f90 +++ b/src/ejecta/ejecta_distance_estimate.f90 @@ -61,6 +61,9 @@ 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 return end subroutine ejecta_distance_estimate