Ejecta modified to be more self-consistent with empirical data on eje…

…cta thickness (calibrated for strengthless soft rock)
MintonGroup · Dec 5, 2016 · ac901fc · ac901fc
1 parent 7e8b6fc
commit ac901fc
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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.