Changed order of when ejrim is called to ensure that it is set before…

… computing the ejecta radius
MintonGroup · Feb 17, 2022 · 3e4c1a5 · 3e4c1a5
1 parent 083ddb1
commit 3e4c1a5
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/crater/crater_dimensions.f90 b/src/crater/crater_dimensions.f90
@@ -68,8 +68,8 @@ subroutine crater_dimensions(user,crater,domain)
 
    ! Calculate the radius where the inner wall meets the original pre-existing surface
    ! This is used to demark the location where excavation transitions to deposition
-   crater%ejrad = max(crater_profile_find_r_inner_wall(user,crater) * crater%frad, crater%rad)
    crater%ejrim = 0.14_DP * (crater%fcrat * 1e-3 * 0.5_DP)**(0.74_DP) * 1e3_DP ! McGetchin et al. (1973) Thickness of ejecta at rim
+   crater%ejrad = max(crater_profile_find_r_inner_wall(user,crater) * crater%frad, crater%rad)
 
    !find rim for counting purposes
    crater%frim = RIMFAC * crater%frad