From 87be35c03b118355aa0a3fccc7f9ef7c00132f47 Mon Sep 17 00:00:00 2001 From: daminton Date: Thu, 9 Feb 2017 19:01:13 +0000 Subject: [PATCH] Added new user input that allows you to truncate the ejecta blanket at a specific radius --- src/ejecta/ejecta_distance_estimate.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/ejecta/ejecta_distance_estimate.f90 b/src/ejecta/ejecta_distance_estimate.f90 index b04f2807..6df703d1 100644 --- a/src/ejecta/ejecta_distance_estimate.f90 +++ b/src/ejecta/ejecta_distance_estimate.f90 @@ -61,7 +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)) - ejdis_estimate = min(ejdis_estimate, crater%fcrat * user%ejecta_truncation) + ejdis_estimate = min(ejdis_estimate, crater%frad * user%ejecta_truncation) return end subroutine ejecta_distance_estimate