Checked out from main branch

MintonGroup · Jul 23, 2021 · 66a2b1a · 66a2b1a
1 parent fbcd4da
commit 66a2b1a
Showing 1 changed file with 3 additions and 3 deletions.
diff --git 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