In the latest dev version of yt, there is a ds.add_deposited_particle_field function attached to all dataset instances that automates <span></span>the workflow Britton suggests.<br><br>On Thursday, June 25, 2015, Britton Smith <<a href="mailto:brittonsmith@gmail.com">brittonsmith@gmail.com</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Brendan,<div><br></div><div>I would have a look in the yt source in yt/fields/particle_fields.py at the field functions like <font face="monospace, monospace">particle_mass</font> and <font face="monospace, monospace">particle_cic</font>.  They are good examples of just how the deposited fields are done.  In addition, there is a field function, <font face="monospace, monospace">_get_density_weighted_deposit_field</font>, for doing mass-weighted deposited fields.  I think this is exactly what you're looking for.  In fact, it seems that mass-weighted velocity fields already exist as something like ("deposit", "<PTYPE>_sum_velocity_x") or ("deposit", "<PTYPE>_nn_velocity_x").  If you have multiple particle types for stars and such, they should already be defined with these fields.  This should give you a template for implementing any other fields of this type.</div><div><br></div><div>Britton</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 24, 2015 at 12:28 PM, Brendan Griffen <span dir="ltr"><<a href="javascript:_e(%7B%7D,'cvml','brendan.f.griffen@gmail.com');" target="_blank">brendan.f.griffen@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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><span><font color="#888888"><div>Brendan</div></font></span></div>
<br>_______________________________________________<br>
yt-users mailing list<br>
<a href="javascript:_e(%7B%7D,'cvml','yt-users@lists.spacepope.org');" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>
</blockquote>