From ac901fcfe7aef6298a81771124a6fff8f7e7f0b8 Mon Sep 17 00:00:00 2001 From: daminton Date: Mon, 5 Dec 2016 21:29:22 +0000 Subject: [PATCH] Ejecta modified to be more self-consistent with empirical data on ejecta thickness (calibrated for strengthless soft rock) --- src/ejecta/ejecta_blanket.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/ejecta/ejecta_blanket.f90 b/src/ejecta/ejecta_blanket.f90 index b0ac702b..232e1ebc 100644 --- a/src/ejecta/ejecta_blanket.f90 +++ b/src/ejecta/ejecta_blanket.f90 @@ -65,7 +65,7 @@ subroutine ejecta_blanket(user,crater,domain,erad,lrad,vejsq,ejang,firstrun) ytrans2 = PItwo * (PIfour**(-THIRD)) ytrans = ytrans1 * (ytrans2**(2._DP / (mu + 2._DP))) ce2asq = ce1sq * abs((user%gaccel * crater%grad * trho) / (tstr + ytrans + vsmall)) - p = 2.0_DP / mu + p = 6.0_DP / mu ce2sq = ce2asq * (crater%grad / crater%rad)**p yterm = ce2sq * (tstr / trho) firstrun = .false.