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.
Here's a script that does it; you can get this with "yt_lodgeit.py
Also, on re-reading your initial question, my initial suggestion of
adding a profile using "ParticleRadius" and "ParticleMassMsun" *will*
function, but it's not what you wanted.
On Tue, Oct 26, 2010 at 8:08 AM, Elizabeth Tasker <firstname.lastname@example.org> wrote:
> Hi Christine and Matt,
> Unfortunately, I don't have a Dark_Matter_Density field, just the position,
> velocities, type, mass and index. (From looking directly at the contents of
> I'm afraid I also need my hand held more to use Matt's suggestion! Would it
> be possible to give me an example of how to call those functions? I tried to
> guess and failed miserably.
> Thank you~
> Matthew Turk wrote:
>> Hi Elizabeth,
>> It should be. You can either do the CIC deposited DM density using
>> the field "particle_density_pyx" (
>> http://yt.enzotools.org/doc/api/field_list.html#particle-density-pyx )
>> if you're on the 'stable' branch of the field "particle_density" if
>> you're on the current development trunk. That will convert the DM to
>> be like a baryon quantity.
>> If you'd rather do raw particle binning, there is also functionality
>> for that by using the x field "ParticleRadius" (
>> http://yt.enzotools.org/doc/api/field_list.html#particleradius ) or
>> one of its derivatives, and the y field "ParticleMassMsun".
>> On Mon, Oct 25, 2010 at 5:27 PM, Elizabeth Tasker <email@example.com>
>>> I'd like to plot dark matter particle density as a function of radius for
>>> particle-only simulation. Can I use add_profile_sphere, but for particles
>>> rather than gas? I only have a single halo and I know the halo it's
>>> in the middle of my box (or it better be, since I put it there!).
>>> Is this possible?
>>> yt-users mailing list
>> yt-users mailing list
> yt-users mailing list
yt-users mailing list