From 0ffc41af5981461167766087aa7d0184f5e20e1c Mon Sep 17 00:00:00 2001 From: David A Minton Date: Tue, 22 Nov 2022 16:34:12 -0500 Subject: [PATCH] Fixed issue where big ejecta blankets were not being truncated properly based on the target body radius. --- src/ejecta/CTEM_debug.code-workspace | 8 -------- src/ejecta/ejecta_emplace.f90 | 2 +- 2 files changed, 1 insertion(+), 9 deletions(-) delete mode 100644 src/ejecta/CTEM_debug.code-workspace diff --git a/src/ejecta/CTEM_debug.code-workspace b/src/ejecta/CTEM_debug.code-workspace deleted file mode 100644 index 7cc6eb86..00000000 --- a/src/ejecta/CTEM_debug.code-workspace +++ /dev/null @@ -1,8 +0,0 @@ -{ - "folders": [ - { - "path": "../.." - } - ], - "settings": {} -} \ No newline at end of file diff --git a/src/ejecta/ejecta_emplace.f90 b/src/ejecta/ejecta_emplace.f90 index 1539bba7..4635e64f 100644 --- 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)