From 64ea812f24002e57c114067f648a71ed53f3f072 Mon Sep 17 00:00:00 2001 From: David Minton Date: Fri, 11 Feb 2022 13:03:35 -0500 Subject: [PATCH] Made sure to compute crater%ejrim before calling crater_profile for the first time --- 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..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