[yt-users] phaseplot set zlim

[Nathan Goldbaum]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20170328/fbd786b1/attachment.html>