Made sure to compute crater%ejrim before calling crater_profile for the first time

MintonGroup · Feb 11, 2022 · 64ea812 · 64ea812
1 parent 79c3f1a
commit 64ea812
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