diff --git a/src/crater/crater_dimensions.f90 b/src/crater/crater_dimensions.f90 index 9a45e04f..fc6de8c9 100644 --- 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