Darn :).  I see what is wrong.  It doesn't seem that the 'unitary' unit is being set up correctly for Athena (it currently ends up being -infinity!).  I'll push a fix in a few minutes.

Sam

On Wed, Feb 27, 2013 at 10:11 AM, nick moeckel wrote:
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 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 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@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org