[yt-users] phaseplot set zlim

Tue Mar 28 09:53:04 PDT 2017

Hmm, that is strange, I didn't notice that earlier. It might be a bug? I'd
need to take a closer look at it. Is there any chance you can share the
script you used to make those two plots?

On Tue, Mar 28, 2017 at 11:40 AM, Stephanie Tonnesen <stonnes at gmail.com>
wrote:

> 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
>>
>>
>
> _______________________________________________
> 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/14afa3df/attachment-0001.htm>