<div dir="ltr"><div><div>Hi Folks,<br><br></div>I'm trying to understand how the particle deposition fields work for grid-based codes. Specifically, how does yt handle particles that are close to the grid boundaries, so that the CIC interpolation stencil involves cells on grid(s) other than the one that contains the particle position?. Like, say I did:<br><br></div><div>import yt<br></div><div>ds = yt.load(fn)<br></div><div>g = ds.index.grids[0]<br></div><div>print g['deposit', 'all_cic']<br><br></div><div>How does yt handle the fact that particles other than those living on "g" can contribute to the total CIC density? Does it grab a layer of "ghost particles" somehow? I can't seem to find where that happens.<br><br></div><div>Sorry if I'm just missing something obvious here...<br><br></div><div>Best,<br>Andrew<br></div></div>