Hi Shankar,
This is something that Enzo actually does -- the Dark_Matter_Density
is generated via the cic_deposit routines, which are cloud-in-cell
interpolation. Where no particles contribute, it gets marked as zero,
which shows up as "bad values" in the log slice, and so it marks them
as white. Looks like everything in this image is working!
-Matt
On Mon, Mar 22, 2010 at 11:06 AM, Agarwal, Shankar
Hi Matt,
I just plotted a slice of Dark_Matter_Density. It looks ok except that there are some white patches (the density legend is VIBGYOR)
http://drop.io/Dark_Matter_Density#
It looks like those white patches are regions from where the density data was not read. Any thoughts?
Shankar
-----Original Message----- From: yt-users-bounces@lists.spacepope.org on behalf of Matthew Turk Sent: Mon 3/22/2010 11:58 AM To: Discussion of the yt analysis package Subject: Re: [yt-users] DM particles or baryon density field ?
Yup, you're right.
On Mon, Mar 22, 2010 at 9:57 AM, Agarwal, Shankar
wrote: Yes, "Dark_Matter_Density field projections and slices" is what I meant. I will try replacing "Density" with "Dark_Matter_Density".
Also, let me know if I am right...
p = pc.add_projection("Density", 0) ------- This is baryon density field.
p.modify["nparticles"](1.0) ------- But this is DM particles (since there is no baryon particles dataset ?).
Shankar
-----Original Message----- From: yt-users-bounces@lists.spacepope.org on behalf of Matthew Turk Sent: Mon 3/22/2010 11:32 AM To: Discussion of the yt analysis package Subject: Re: [yt-users] DM particles or baryon density field ?
Hi Shankar,
I guess the question is sort of ill-defined. Dark matter is discretized as particles; the idea of projecting dark matter is not the same as projecting a field quantity. The same goes for plotting slices.
However, you can make dark matter plots. You cna use the "add_particles" method on the PlotCollection, which accepts an axis and a width -- what it does is plot all the particles within that rectangular prism. This can be mocked up as a projection by setting the width equal to the domain width along that axis.
You can also try using "Dark_Matter_Density", which might be easier to visualize, as a projection or a slice.
If you run into any problems with the particle plot, let us know -- it's not as well-used as the other plots, so there might be subtle bugs. :)
-Matt
On Mon, Mar 22, 2010 at 9:27 AM, Agarwal, Shankar
wrote: Hi,
The projection, slice, velocity and overplot particles routines: They read the DM particles or gas ?
Looking at the density legend, it looks like this is baryon density field.
Is there a provision to plot projection and slice for DM ?
Also, I get this error when overplotting particles...
yt INFO 2010-03-22 12:17:44,051 Created plot collection with default plot-center = [0.5, 0.5, 0.5] yt INFO 2010-03-22 12:17:44,144 Getting the binary hierarchy yt INFO 2010-03-22 12:17:44,153 Finished with binary hierarchy reading yt INFO 2010-03-22 12:17:46,108 Added slice of Density at x = 0.5 with 'center' = [0.5, 0.5, 0.5] Traceback (most recent call last): File "slice.py", line 42, in <module> p.modify["nparticles"](0.005) File "/autohome/u118/agarwa23/yt-x86_64/src/yt-1.6/yt/raven/PlotTypes.py", line 52, in __getitem__ raise KeyError(item) KeyError: 'nparticles'
Shankar _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org