Added 0 value buffer to edges of ejecta diffusion zone

MintonGroup · Dec 5, 2016 · 9c9d21c · 9c9d21c
1 parent 2a9afeb
commit 9c9d21c
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/ejecta/ejecta_soften.f90 b/src/ejecta/ejecta_soften.f90
@@ -47,6 +47,8 @@ subroutine ejecta_soften(user,surf,N,indarray,cumulative_elchange)
 
    ! 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 
+   ! 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
    kdiff(:,1) = 0.0_DP
    kdiff(N,:) = 0.0_DP