From 602d8fe4b6b567dc232fe96a0b237d9b0dfb6201 Mon Sep 17 00:00:00 2001 From: daminton Date: Thu, 9 Feb 2017 18:56:53 +0000 Subject: [PATCH] Overhauled crater formation routine. All raised rims are now part of the ejecta blanket --- src/ejecta/ejecta_interpolate.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/ejecta/ejecta_interpolate.f90 b/src/ejecta/ejecta_interpolate.f90 index 8d53dd18..a8ba7931 100644 --- a/src/ejecta/ejecta_interpolate.f90 +++ b/src/ejecta/ejecta_interpolate.f90 @@ -39,9 +39,9 @@ subroutine ejecta_interpolate(crater,domain,lrad,ejb,ejtble,ebh,vsq,theta,erad,m ! Executable code ! Locate ourselves in the table - outeredge = crater%frad + domain%ejbres * EJBTABSIZE - inneredge = crater%frad - k = max(min(1 + int((lrad - inneredge) / (outeredge - inneredge) * (EJBTABSIZE - 1.0_DP)),ejtble),1) + inneredge = crater%rad + outeredge = crater%rad * exp(domain%ejbres * EJBTABSIZE) + k = max(min(1 + int((log(lrad) - log(inneredge)) / (log(outeredge) - log(inneredge)) * (EJBTABSIZE - 1.0_DP)),ejtble),1) loglrad = log(lrad) logtablerad = ejb(k)%lrad