Updated ejecta soften with better data on degradation

MintonGroup · Dec 6, 2016 · b02e6af · b02e6af
1 parent a623e45
commit b02e6af
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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