To expand a little on what Matt said ealier, if you don't have the Dark_Matter_Density field, you should still be able to use the field particle_density_pyx.  This uses a yt implementation of the cic_deposit routine to take the particle data and make a grid field out of it.  Depending on the version of yt you're using, the field will either be called particle_density_pyx or just particle_density.<br>
<br>Britton<br><br><div class="gmail_quote">On Tue, Oct 26, 2010 at 12:07 PM, Matthew Turk <span dir="ltr"><<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Hi Elizabeth,<br>
<br>
Here's a script that does it; you can get this with "yt_lodgeit.py<br>
--download=1294".<br>
<br>
<a href="http://paste.enzotools.org/show/1294/" target="_blank">http://paste.enzotools.org/show/1294/</a><br>
<br>
Also, on re-reading your initial question, my initial suggestion of<br>
adding a profile using "ParticleRadius" and "ParticleMassMsun" *will*<br>
function, but it's not what you wanted.<br>
<br>
Best,<br>
<br>
Matt<br>
<div><div></div><div class="h5"><br>
On Tue, Oct 26, 2010 at 8:08 AM, Elizabeth Tasker <<a href="mailto:taskere@mcmaster.ca">taskere@mcmaster.ca</a>> wrote:<br>
> Hi Christine and Matt,<br>
><br>
> Unfortunately, I don't have a Dark_Matter_Density field, just the position,<br>
> velocities, type, mass and index. (From looking directly at the contents of<br>
> XXXX.cpu0000).<br>
><br>
> I'm afraid I also need my hand held more to use Matt's suggestion! Would it<br>
> be possible to give me an example of how to call those functions? I tried to<br>
> guess and failed miserably.<br>
><br>
> Thank you~<br>
><br>
> Elizabeth<br>
><br>
><br>
><br>
><br>
> Matthew Turk wrote:<br>
>><br>
>> Hi Elizabeth,<br>
>><br>
>> It should be.  You can either do the CIC deposited DM density using<br>
>> the field "particle_density_pyx" (<br>
>> <a href="http://yt.enzotools.org/doc/api/field_list.html#particle-density-pyx" target="_blank">http://yt.enzotools.org/doc/api/field_list.html#particle-density-pyx</a> )<br>
>> if you're on the 'stable' branch of the field "particle_density" if<br>
>> you're on the current development trunk.  That will convert the DM to<br>
>> be like a baryon quantity.<br>
>><br>
>> If you'd rather do raw particle binning, there is also functionality<br>
>> for that by using the x field "ParticleRadius" (<br>
>> <a href="http://yt.enzotools.org/doc/api/field_list.html#particleradius" target="_blank">http://yt.enzotools.org/doc/api/field_list.html#particleradius</a> ) or<br>
>> one of its derivatives, and the y field "ParticleMassMsun".<br>
>><br>
>> -Matt<br>
>><br>
>> On Mon, Oct 25, 2010 at 5:27 PM, Elizabeth Tasker <<a href="mailto:taskere@mcmaster.ca">taskere@mcmaster.ca</a>><br>
>> wrote:<br>
>><br>
>>><br>
>>> Hi,<br>
>>><br>
>>> I'd like to plot dark matter particle density as a function of radius for<br>
>>> my<br>
>>> particle-only simulation. Can I use add_profile_sphere, but for particles<br>
>>> rather than gas? I only have a single halo and I know the halo it's<br>
>>> centred<br>
>>> in the middle of my box (or it better be, since I put it there!).<br>
>>><br>
>>> Is this possible?<br>
>>><br>
>>> Thanks,<br>
>>><br>
>>> Elizabeth<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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
>>><br>
>>><br>
>><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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
>><br>
><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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
><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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</div></div></blockquote></div><br>