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