Fixed bug in crater_subpixel_diffusion where the kdiff array was not …

…being calculated for the border cells properly
MintonGroup · Feb 23, 2017 · 0853771 · 0853771
1 parent c72f600
commit 0853771
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/crater/crater_subpixel_diffusion.f90 b/src/crater/crater_subpixel_diffusion.f90
@@ -70,8 +70,8 @@ subroutine crater_subpixel_diffusion(user,surf,prod,nflux,domain,finterval,kdiff
 
    end do
 
-   do j = 1,user%gridsize
-      do i = 1,user%gridsize
+   do j = 0,user%gridsize + 1
+      do i = 0,user%gridsize + 1
          do n = 1,ntot
             dburial = EXFAC * 0.5_DP * nflux(1,n)
             if (surf(i,j)%ejcov > dburial) then