Set initial guess for erad to transient radius crater%rad, instead of…

… just inward of the innermost lrad value
MintonGroup · Feb 1, 2022 · 0cc2b1b · 0cc2b1b
1 parent 988a780
commit 0cc2b1b
Showing 1 changed file with 1 addition and 2 deletions.
diff --git 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