From 0cc2b1beffc157b2c1aa2817294ace0e89ee6ae9 Mon Sep 17 00:00:00 2001
From: David A Minton <daminton@purdue.edu>
Date: Tue, 1 Feb 2022 14:05:31 -0500
Subject: [PATCH] Set initial guess for erad to transient radius crater%rad,
 instead of just inward of the innermost lrad value

---
 src/ejecta/ejecta_table_define.f90 | 3 +--
 1 file changed, 1 insertion(+), 2 deletions(-)

diff --git a/src/ejecta/ejecta_table_define.f90 b/src/ejecta/ejecta_table_define.f90
index 04451325..1a91af58 100644
--- a/src/ejecta/ejecta_table_define.f90
+++ b/src/ejecta/ejecta_table_define.f90
@@ -49,7 +49,7 @@ subroutine ejecta_table_define(user,crater,domain,ejb,ejtble,melt)
    domain%ejbres = (log(crater%ejdis) - log(crater%ejrad)) / EJBTABSIZE
    lrad = crater%ejrad
    lrad = exp(log(lrad) + domain%ejbres)
-   erad = crater%ejrad 
+   erad = crater%rad
    ejtble = EJBTABSIZE
    firstrun = .true.
    thick = 0._DP
@@ -64,7 +64,6 @@ subroutine ejecta_table_define(user,crater,domain,ejb,ejtble,melt)
          vimp = crater%impvel
       end if 
    
-
       call regolith_melt_zone(user,crater,dimp,vimp,rmelt,depthb)
    end if