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

…cta thickness (calibrated for strengthless soft rock)
MintonGroup · Dec 5, 2016 · 96034f7 · 96034f7
1 parent 584cc82
commit 96034f7
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/src/ejecta/ejecta_distance_estimate.f90 b/src/ejecta/ejecta_distance_estimate.f90
@@ -61,6 +61,9 @@ subroutine ejecta_distance_estimate(user,crater,domain,ejdis_estimate)
    ! Extrapolate in logspace out to the minimum thickness
    slope = log(lrad2 / lrad1) / log(thick2 / thick1)
    ejdis_estimate = exp(slope * log(domain%small / thick1) + log(lrad1))
+   if (.not.user%discontinuous) then
+      ejdis_estimate = min(ejdis_estimate, 2.3_DP * crater%frad**(1.006_DP))
+   end if
 
    return
 end subroutine ejecta_distance_estimate