[yt-users] athena region problem

[nick moeckel]

Good deal, thanks!


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

> I've posted a PR here:
>
> https://bitbucket.org/yt_analysis/yt/pull-request/448/bugfix-for-units-setup-in-athena/diff
>
> 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 at gmail.com>wrote:
>
>> 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 <nickolas1 at 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 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
>>>>
>>>>
>>>
>>> _______________________________________________
>>> 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/ff2d9519/attachment.html>