From 055f18c3153cae6d1bcc4133515deaef9ba1bb3c Mon Sep 17 00:00:00 2001 From: David Minton Date: Thu, 17 Feb 2022 16:40:37 -0500 Subject: [PATCH] Fixed McGetchin equation as it was already supposed to be in units of meters --- src/crater/crater_dimensions.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/crater/crater_dimensions.f90 b/src/crater/crater_dimensions.f90 index 9a45e04f..f790cb3a 100644 --- a/src/crater/crater_dimensions.f90 +++ b/src/crater/crater_dimensions.f90 @@ -69,7 +69,7 @@ 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%ejrim = 0.14_DP * (crater%fcrat * 0.5_DP)**(0.74_DP) ! McGetchin et al. (1973) Thickness of ejecta at rim !find rim for counting purposes crater%frim = RIMFAC * crater%frad