From 7e8b6fcde62ca5bf7e76421216a51bf2959a22a6 Mon Sep 17 00:00:00 2001 From: daminton Date: Mon, 5 Dec 2016 21:29:17 +0000 Subject: [PATCH] Ejecta modified to be more self-consistent with empirical data on ejecta thickness (calibrated for strengthless soft rock) --- src/ejecta/ejecta_thickness.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/ejecta/ejecta_thickness.f90 b/src/ejecta/ejecta_thickness.f90 index 10f42481..0928d47a 100644 --- a/src/ejecta/ejecta_thickness.f90 +++ b/src/ejecta/ejecta_thickness.f90 @@ -45,7 +45,7 @@ subroutine ejecta_thickness(user,crater,erad1,erad2,lrad1,lrad2,thick) if (crater%imp < user%basinimp) then ! compute ejecta thickness ce1sq = (2._DP / (CT * CT)) * mu**2 / (1._DP + 2 * mu + mu**2) - excav = abs(0.5_DP * ce1sq * PI * (erad2**3 - erad1**3)) + excav = 7.43_DP * (crater%fcrat)**(-0.26_DP) * abs(0.5_DP * ce1sq * PI * (erad2**3 - erad1**3)) else edDt = 1.150_DP ! Calibrated to Orientale using Fasset et al. 2011 excav = abs(0.5_DP * edDt * PI * (erad2**3 - erad1**3))