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