From 66a2b1a550d9f8b745f2d8dfe16c32b32fc89dd1 Mon Sep 17 00:00:00 2001 From: David Minton Date: Fri, 23 Jul 2021 14:35:35 -0400 Subject: [PATCH] Checked out from main branch --- 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 79c15e86..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 - 0.5_DP) - inneredge = crater%frad + 0.5_DP * domain%ejbres - 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