From 61faacc8f217847ea6b1800decb770c415b46402 Mon Sep 17 00:00:00 2001 From: Chris Vogl <cvogl@uw.edu> Date: Tue, 14 Nov 2017 19:55:35 -0800 Subject: [PATCH] Added buffer in 3d/b4step3.f90 to avoid cells at edge of fault being assigned 0 slip --- 3d/b4step3.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/3d/b4step3.f90 b/3d/b4step3.f90 index 97f19bb..22fc61a 100644 --- a/3d/b4step3.f90 +++ b/3d/b4step3.f90 @@ -24,11 +24,11 @@ subroutine b4step3(mbc,mx,my,mz,meqn,q,xlower,ylower,zlower, & do j=1-mbc,my+mbc ycell = ylower + (j-0.5d0)*dy - if (ycb(1) <= ycell - 0.5d0*dy .and. ycell + 0.5d0*dy <= ycb(2)) then + if (ycb(1)-1.d-10 <= ycell - 0.5d0*dy .and. ycell + 0.5d0*dy <= ycb(2)+1.d-10) then do i=1-mbc,mx+mbc xcell = xlower + (i-0.5d0)*dx - if (xcb(1) <= xcell - 0.5d0*dx .and. xcell + 0.5d0*dx <= xcb(2)) then + if (xcb(1)-1.d-10 <= xcell - 0.5d0*dx .and. xcell + 0.5d0*dx <= xcb(2)+1.d-10) then do l=1,nsubfaults if (subfaults(l)%longitude*LAT2METER <= xcell .and. &