<div dir="ltr">Hi,<div><br></div><div>I'm looking to convert particle based quantities into a uniform grid based quantity (arepo data is input). So say I have a collection of gas elements with some density rho and x,y,z positions. I want to construct an 256^3 mesh representing the cells in a mesh.</div><div><br></div><div>x,y,z,rho --> cic mesh where the values at i,j,k is the smoothed density are created by the cic algorithm</div><div><br></div><div>Now, generalising this more, can I have rho to be any quantity? So now say I have star particles which similarly have x,y,z but I want to distribute their luminosity. Again, I want to do it for say the gas velocity field as well. Is there any general functionality for converting particle quantities to mesh quantites in yt with the deposit all_cic functionality? I have a basic code for just x,y,z particle data and distributing that based on density to a cic (see below)</div><div><br></div><div>ad = pf.all_data()</div><div>level = int(math.log(ndim,2))</div><div>cg = pf.covering_grid(level=level, left_edge=[0,0,0],dims=[ndim,ndim,ndim])</div><div>arr = cg['deposit', 'all_cic'],</div><div><br></div><div>but I'm not sure what to do for the other variables. Any help is appreciated in advance.</div><div><br></div><div>Thanks.</div><div>Brendan</div></div>