From e45f4c452044b0777796d5e1f61c82a24dc86a6d Mon Sep 17 00:00:00 2001 From: daminton Date: Thu, 9 Feb 2017 19:01:23 +0000 Subject: [PATCH] Added new user input that allows you to truncate the ejecta blanket at a specific radius --- src/ejecta/ejecta_emplace.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/ejecta/ejecta_emplace.f90 b/src/ejecta/ejecta_emplace.f90 index bc01a30f..53c59b60 100644 --- a/src/ejecta/ejecta_emplace.f90 +++ b/src/ejecta/ejecta_emplace.f90 @@ -145,7 +145,7 @@ subroutine ejecta_emplace(user,surf,crater,domain,ejb,ejtble,deltaMtot) ! determine area to effect continuous = 2.348_DP * crater%frad**(1.006_DP) - inc = max(min(nint(crater%fcrat * user%ejecta_truncation / user%pix) + 1,PBCLIM*user%gridsize),1) + inc = max(min(nint(crater%frad * user%ejecta_truncation / user%pix) + 1,PBCLIM*user%gridsize),1) maxdistance = inc * user%pix ! Increase the box a bit to take into account possible ejecta pattern distortion due to topography inc = ceiling(inc * 1.5_DP)