From 2a9afeb95ff5ea4a3c1c43bf58fa955ea9c9a684 Mon Sep 17 00:00:00 2001
From: daminton <daminton@purdue.edu>
Date: Mon, 5 Dec 2016 22:18:54 +0000
Subject: [PATCH] Added 0 value buffer to edges of ejecta diffusion zone

---
 src/ejecta/ejecta_soften.f90 | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/src/ejecta/ejecta_soften.f90 b/src/ejecta/ejecta_soften.f90
index b570537e..f4b01e71 100644
--- a/src/ejecta/ejecta_soften.f90
+++ b/src/ejecta/ejecta_soften.f90
@@ -47,6 +47,11 @@ 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 
+   kdiff(1,:) = 0.0_DP
+   kdiff(:,1) = 0.0_DP
+   kdiff(N,:) = 0.0_DP
+   kdiff(:,N) = 0.0_DP
+
    !TESTING
    !kdiff = user%sf * user%pix**user%k * ebharr**user%p !**0.900