<div dir="ltr">And maybe you need to clarify what you mean by "loading the data using my tools". Are you ever creating a yt Dataset object? If not, then Matt's proviso about using yt's internal APIs applies. If you're creating a fake in memory yt dataset using e.g. load_particles, then you should be able to use 'add_deposited_particle_field` as I suggested in my e-mail above.</div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Aug 19, 2015 at 1:08 PM, Matthew Turk <span dir="ltr"><<a href="mailto:matthewturk@gmail.com" target="_blank">matthewturk@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 Brendan,<br><div><br></div><div>The CIC deposition routines *themselves* are subject to change and not necessarily API stable. You can try using them and the raw routines, but they may be subject to changes.</div><div><br></div><div>-Matt</div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Aug 19, 2015 at 1:02 PM, Brendan Griffen <span dir="ltr"><<a href="mailto: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 Nathan,<div><br></div><div>Is there an example of how one might do this for say an array which contains particle based x,y,z (for an arepo cell) and some arbitrary other quantity of the same length (e.g. HI density, velocity magnitudes etc.)? Seems like you can only do it if you've loaded the data through yt. I wanted to load in the quantities using my tools as numpy arrays then feed that into the deposition.</div><div><br></div><div>Thanks.</div><span><font color="#888888"><div>Brendan</div></font></span></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 25, 2015 at 11:31 AM, Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">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.<div><div><br><br>On Thursday, June 25, 2015, Britton Smith <<a href="mailto:brittonsmith@gmail.com" target="_blank">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>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>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>
</div></div><br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto: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>
</div></div><br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto: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>
</div></div><br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">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>