[yt-users] athena region problem

[nick moeckel]

in my haste to get an example my zeroes didn't have decimals there (shame
on me), but that doesn't fix it.

In [37]: lc = [0.0, 0.0, 0.4]
In [38]: rc = [1.0, 1.0, 0.6]
In [39]: region = pf.h.region(center, lc, rc)
In [40]: region.volume()
Out[40]: -inf


On Wed, Feb 27, 2013 at 6:06 PM, Sam Skillman <samskillman at gmail.com> wrote:

> Hi Nick,
>
> It looks like lc only has 2 values in it.  If fixing that to lc = [0.0,
> 0.0, 0.4] works (or still breaks), let us know.
>
> Cheers,
> Sam
>
>
> On Wed, Feb 27, 2013 at 10:02 AM, nick moeckel <nickolas1 at gmail.com>wrote:
>
>> Hi there,
>>
>> I'm dealing with some athena data right now, and running into some issues
>> with defining regions.
>>
>> an example, with domain_left_edge = [0, 0, 0], and right_edge = [1, 1, 1]
>>
>> In [18]: center=[0.5,0.5,0.5]
>> In [19]: lc = [0,0,0.4]
>> In [20]: rc=[1.0,1.0,0.6]
>> In [21]: region = pf.h.region(center, lc, rc)
>> In [22]: region.volume()
>> Out[22]: -inf
>>
>> I'm able to get at region.quantities, but trying to make a projection
>> plot fails.
>>
>> Interestingly, or maybe not, if I try to do something similar with data
>> that has its domain centered on the origin, the volume is finite and
>> everything I've tried with plotting works fine. I guess athena support is
>> still a bit preliminary, but I thought I'd throw this out there.
>>
>> best,
>> Nick
>>
>> _______________________________________________
>> 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/20130227/d2f12055/attachment.html>