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
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@gmail.com> wrote:
I've posted a PR here: https://bitbucket.org/yt_analysis/yt/pull-request/448/bugfix-for-units-setup... When it is accepted it should do the trick. Thanks for pointing this out! Cheers, Sam On Wed, Feb 27, 2013 at 10:21 AM, Sam Skillman <samskillman@gmail.com>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 <samskillman@gmail.com> wrote:
I've posted a PR here: https://bitbucket.org/yt_analysis/yt/pull-request/448/bugfix-for-units-setup... When it is accepted it should do the trick. Thanks for pointing this out! Cheers, Sam On Wed, Feb 27, 2013 at 10:21 AM, Sam Skillman <samskillman@gmail.com>wrote:
participants (2)
-
nick moeckel -
Sam Skillman