[yt-users] convert arbitrary particle based quantities to cloud-in-cell mesh

Brendan Griffen brendan.f.griffen at gmail.com
Wed Jun 24 04:28:02 PDT 2015


Hi,

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.

x,y,z,rho --> cic mesh where the values at i,j,k is the smoothed density
are created by the cic algorithm

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)

ad = pf.all_data()
level = int(math.log(ndim,2))
cg = pf.covering_grid(level=level, left_edge=[0,0,0],dims=[ndim,ndim,ndim])
arr = cg['deposit', 'all_cic'],

but I'm not sure what to do for the other variables. Any help is
appreciated in advance.

Thanks.
Brendan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20150624/f4a47769/attachment.htm>


More information about the yt-users mailing list