From 44bdf5865b6c3488b70a5a08770efdb0a7d59e11 Mon Sep 17 00:00:00 2001 From: Michael Krayer Date: Thu, 5 Aug 2021 01:08:37 +0200 Subject: [PATCH] needed another cell in masking bounding box --- particle.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/particle.py b/particle.py index c3af2a8..ce562b3 100644 --- a/particle.py +++ b/particle.py @@ -152,7 +152,7 @@ class Particles: rp = self.attr['r'][ipart]+padding for pos in self.position_with_duplicates(ipart,padding=padding): idxlo = np.array(fld.nearest_gridpoint(pos-rp,lower=True)) - sfdim = np.ceil(2*rp/fld.spacing+1).astype('int') + sfdim = np.ceil(2*rp/fld.spacing+2).astype('int') sf = fld.extract_subfield(idxlo,sfdim,deep=False,strict_bounds=False) xsf = sf.x().reshape((-1,1,1)) ysf = sf.y().reshape((1,-1,1))