From 0853771c66b908a3aabb1e13acb18c6b4e412188 Mon Sep 17 00:00:00 2001 From: daminton Date: Thu, 23 Feb 2017 03:24:50 +0000 Subject: [PATCH] Fixed bug in crater_subpixel_diffusion where the kdiff array was not being calculated for the border cells properly --- src/crater/crater_subpixel_diffusion.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/crater/crater_subpixel_diffusion.f90 b/src/crater/crater_subpixel_diffusion.f90 index eca2f03d..a95e94be 100644 --- 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