Fixed bug in ejecta_table_define that was causing the root finding al…

…gorithm to fail to determine the ejection radious of the first element of the ejecta table
MintonGroup · Jan 27, 2022 · 338cd23 · 338cd23
1 parent 66816b7
commit 338cd23
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/ejecta/ejecta_table_define.f90 b/src/ejecta/ejecta_table_define.f90
@@ -48,7 +48,7 @@ subroutine ejecta_table_define(user,crater,domain,ejb,ejtble,melt)
                                                         ! We go out a factor of 3 to get the discontinuous ejecta thickness 
    domain%ejbres = (log(crater%ejdis) - log(crater%ejrad)) / EJBTABSIZE
    lrad = crater%ejrad 
-   erad = crater%ejrad / 2
+   erad = crater%rad 
    ejtble = EJBTABSIZE
    firstrun = .true.
    thick = 0._DP