From cda6d2e35924e8f5bf3bceadda2256650085dddb Mon Sep 17 00:00:00 2001 From: daminton Date: Fri, 3 Feb 2017 16:05:10 +0000 Subject: [PATCH] Corrected ejecta interpolation algorithm for new ejecta table locations --- src/ejecta/ejecta_interpolate.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/ejecta/ejecta_interpolate.f90 b/src/ejecta/ejecta_interpolate.f90 index 79c15e86..8d53dd18 100644 --- a/src/ejecta/ejecta_interpolate.f90 +++ b/src/ejecta/ejecta_interpolate.f90 @@ -39,8 +39,8 @@ 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 + outeredge = crater%frad + domain%ejbres * EJBTABSIZE + inneredge = crater%frad k = max(min(1 + int((lrad - inneredge) / (outeredge - inneredge) * (EJBTABSIZE - 1.0_DP)),ejtble),1) loglrad = log(lrad) logtablerad = ejb(k)%lrad