[yt-users] Cutting region confusion

[Matthew Turk]

Hi Stephanie,

What do they look like if you project along all three axes and use these
objects as sources?

Matt
On Aug 4, 2014 1:34 PM, "Stephanie Tonnesen" <stonnes at gmail.com> wrote:

> Hi John,
>
> When I change the tracerp25 selections to all have <= or >= symbols, the
> value basically doesn't change.  With the gt/lt and equal-to "Ones" is
> definitely different:
>
> tracerp25 is 1099285
> tdp25 is 2028955
>
> So the selection criteria is definitely more stringent on
> tracerp25 = alld.cut_region(["grid['
> specific_scalar[0]'] > 0.25","grid['z'] <= 0.1925","grid['z'] >=
> -0.1925","grid['cyl_Rcode'] <= 1.2"])
>
> -Stephanie
>
>
>
> On Mon, Aug 4, 2014 at 2:11 PM, John ZuHone <jzuhone at gmail.com> wrote:
>
>> Hi Stephanie,
>>
>> I suspect the problem is that these are not *exactly* the same regions.
>> Have you tried checking the result of:
>>
>> tdp25["Ones"].sum()
>>
>> tracerp25["Ones"].sum()
>>
>> which would be a quick way to determine the number of cells in each
>> region.
>>
>> What might be going on (but I have yet to verify) is that the disk
>> selection container object is using >= and <= for the zone selection,
>> whereas you used > and < in the definition of your cut region. You might
>> try swapping those out and seeing if that changes things.
>>
>> Best,
>>
>> John Z
>>
>> On Aug 4, 2014, at 1:15 PM, Stephanie Tonnesen <stonnes at gmail.com> wrote:
>>
>> > Hi all,
>> >
>> > I am trying to get the total mass in a disk region of my simulation.  I
>> am using athena data and yt 2.7-dev.  I thought I was defining the region
>> in the same way using two methods:
>> >
>> > alld = pf.h.all_data()
>> > tracerp25 = alld.cut_region(["grid['specific_scalar[0]'] >
>> 0.25","grid['z'] < 0.1925","grid['z'] > -0.1925","grid['cyl_Rcode'] < 1.2"])
>> >
>> > disk = pf.h.disk([0.0,0.0,0.0],[0.0,0.0,1.0],1.2,0.385)
>> > tdp25 = disk.cut_region(["grid['specific_scalar[0]'] > 0.25"])
>> >
>> > but when I print the total mass I get answers that differ by 10%
>> >
>> > print tracerp25.quantities["TotalQuantity"]("CellMassMsun")
>> > print tdp25.quantities["TotalQuantity"]("CellMassMsun")
>> >
>> > >>[8019047231.5192223]
>> > >>[8886968784.4710732]
>> >
>> > Can someone point out where my difference is?
>> >
>> > Thanks!
>> >
>> > Stephanie
>> > _______________________________________________
>> > 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/20140804/cef0ac73/attachment.htm>