From b02e6aff3039fdef46d5d1aaeb27e4ea9b826a71 Mon Sep 17 00:00:00 2001 From: daminton Date: Tue, 6 Dec 2016 14:32:39 +0000 Subject: [PATCH] Updated ejecta soften with better data on degradation --- src/ejecta/ejecta_soften.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/ejecta/ejecta_soften.f90 b/src/ejecta/ejecta_soften.f90 index 8ea3d65f..f26f7d2c 100644 --- a/src/ejecta/ejecta_soften.f90 +++ b/src/ejecta/ejecta_soften.f90 @@ -32,7 +32,7 @@ subroutine ejecta_soften(user,surf,N,indarray,cumulative_elchange) integer(I4B),dimension(2,N,N),intent(in) :: indarray real(DP),dimension(N,N),intent(inout) :: cumulative_elchange !real(DP),parameter :: SOFTEN_FACTOR = 2.50_DP ! Constant in topographic diffusion term for ejecta blanket softening - real(DP),parameter :: EJECTA_SOFTEN_FACTOR = 2.00_DP ! Constant in topographic diffusion term for ejecta blanket softening + real(DP),parameter :: EJECTA_SOFTEN_FACTOR = 1.50_DP ! Constant in topographic diffusion term for ejecta blanket softening ! Internal variables integer(I4B) :: maxhits @@ -46,7 +46,7 @@ subroutine ejecta_soften(user,surf,N,indarray,cumulative_elchange) maxhits = 1 ! Diffusion constant for 1 time unit was found to be proportional to ejecta thickness times the pixel size - kdiff = EJECTA_SOFTEN_FACTOR * user%pix * ebharr + kdiff = EJECTA_SOFTEN_FACTOR * ebharr**2 ! This prevents an infinite hole bug by making a 0 value buffer along the ! edges of the diffusion constant matrix kdiff(1,:) = 0.0_DP