Fixed issue where big ejecta blankets were not being truncated proper…

…ly based on the target body radius.
MintonGroup · Nov 22, 2022 · 0ffc41a · 0ffc41a
1 parent 87db67b
commit 0ffc41a
Show file tree

Hide file tree

Showing 2 changed files with 1 addition and 9 deletions.
diff --git a/src/ejecta/CTEM_debug.code-workspace b/src/ejecta/CTEM_debug.code-workspace
diff --git a/src/ejecta/ejecta_emplace.f90 b/src/ejecta/ejecta_emplace.f90
@@ -144,7 +144,7 @@ subroutine ejecta_emplace(user,surf,crater,domain,ejb,ejtble,deltaMtot,age,age_r
    inc = ceiling(inc * 1.5_DP)
    krad = user%ejecta_truncation * crater%frad
    dradsq = int(krad / user%pix) + 3
-   inc = max(inc,dradsq)
+   inc = min(nint(PI * user%trad / user%pix),max(inc,dradsq)) ! Ensure that the ejecta doesn't get any bigger than the surface can accomodate
    dradsq = dradsq**2
 
    if (user%dosoftening) kdiffmax = crater_degradation_function(user,crater%frad)