[yt-users] phaseplot set zlim

Stephanie Tonnesen stonnes at gmail.com
Tue Mar 28 09:40:27 PDT 2017


Awesome!  This looks like exactly what I want to do, thanks!

But here was my problem with the histograms--the gap made sense, just like
you said.  What was really strange was that more of the histogram bins were
yellow in the "mcut" histogram--which has less mass to work with!  So why
did cutting out the low mass gas make more of the histogram have higher
masses (the yellow)?  This isn't dire, I am just wondering so I can
understand better how the code works.

Thanks,
Stephanie


--
Dr. Stephanie Tonnesen
Alvin E. Nashman Postdoctoral Fellow
Carnegie Observatories, Pasadena, CA
stonnes at gmail.com

On Tue, Mar 28, 2017 at 6:39 AM, Nathan Goldbaum <nathan12343 at gmail.com>
wrote:

> Ah sorry, I gave you bad advice in my initial e-mail.
>
> The cut_region selects individual cells in your simulation that have
> masses greater than the cutoff you chose. However that's not what the
> profile is giving you. It's showing you regions in density-temperature
> parameter space, with each box in this parameter space getting
> contributions from several zones in your simulation. The reason your
> PhasePlot made from the cut_region has "holes" you weren't expecting is
> that those regions are the ones preferentially filled by low-mass zones in
> your simulation.
>
> To get the plot you want you're going to need to mask out the 2D profile
> PhasePlot uses to create the 2D histogram plot. The most straightforward
> way to do this is to create the profile object manually using
> create_profile, tweak it by zeroing out histogram bins that don't contain
> enough mass, and then create a PhasePlot from the profile object using
> PhasePlot.from_profile.
>
> Here's a worked example showing how to do this:
>
> https://gist.github.com/3a2afd006be4d1918f28d98ef2a1e077
>
> Again, apologies for giving you bad advice on this issue initially, hope
> this example helps.
>
> -Nathan
>
> On Tue, Mar 28, 2017 at 1:32 AM Stephanie Tonnesen <stonnes at gmail.com>
> wrote:
>
> Just to be clear, I *do* have weight_field=none set.
>
>
>
> --
> Dr. Stephanie Tonnesen
> Alvin E. Nashman Postdoctoral Fellow
> Carnegie Observatories, Pasadena, CA
> stonnes at gmail.com
>
> On Mon, Mar 27, 2017 at 11:29 PM, Stephanie Tonnesen <stonnes at gmail.com>
> wrote:
>
> Hi Nathan and yt-users,
>
> So I used the cut_region to only select cells with cell_mass > 1e36, and
> got a surprising plot (I think).  I am attaching both 2d profiles here for
> comparison.  I am not surprised that there are unfilled regions of the
> figure--those may have been filled with cells that had less than 1e36 g,
> although the positions are a bit strange.  What I find really surprising is
> that so much more of the figure that has the mass cut is yellow--shoudn't
> there be equal or less mass than when I am not making a cut?
>
> Thanks,
> Stephanie
>
>
> --
> Dr. Stephanie Tonnesen
> Alvin E. Nashman Postdoctoral Fellow
> Carnegie Observatories, Pasadena, CA
> stonnes at gmail.com
>
> On Sun, Mar 26, 2017 at 6:38 PM, Stephanie Tonnesen <stonnes at gmail.com>
> wrote:
>
> Hi Nathan,
>
> Thanks!  I was hoping there was a way to do this within the plotting
> commands, but that will definitely work.
>
> Stephanie
>
> On Sun, Mar 26, 2017 at 4:21 PM Nathan Goldbaum <nathan12343 at gmail.com>
> wrote:
>
> Hi Stephanie,
>
> One way to do this would be to create the PhasePlot using a cut_region
> data object. If you're already using another data object to construct your
> phase plot, the cut_region could be passed to the other data object via the
> data_source keyword argument.
>
> Nathan
>
> On Sun, Mar 26, 2017 at 4:24 PM Stephanie Tonnesen <stonnes at gmail.com>
> wrote:
>
> Hi yt-users,
>
> I am trying to make some nice looking phaseplots, and want to set the mass
> range that I show in the figure.  I use set_zlim and instead of only
> showing the cells that call into that mass range, it just changes the color
> coding.  You can see the attached for what I am talking about.  Is there a
> way to have phaseplot only show me the cells that have "cell_mass" larger
> that 1e36 (or whatever I happen to set)?
>
> Thanks!
>
> Stephanie
>
> --
> Dr. Stephanie Tonnesen
> Alvin E. Nashman Postdoctoral Fellow
> Carnegie Observatories, Pasadena, CA
> stonnes at gmail.com
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
> --
> --
> Dr. Stephanie Tonnesen
> Alvin E. Nashman Postdoctoral Fellow
> Carnegie Observatories, Pasadena, CA
> stonnes at gmail.com
>
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20170328/b95b1e0d/attachment.html>


More information about the yt-users mailing list