[yt-dev] Dealing with overlapping grids in Athena data

Matthew Turk matthewturk at gmail.com
Thu Aug 22 07:03:41 PDT 2013


Hi John,

On Thu, Aug 22, 2013 at 9:02 AM, John ZuHone <jzuhone at gmail.com> wrote:
> Hi Matt,
>
>> Thanks for sending this.  My reading is that essentially the problem
>> is that Athena grids at level L=N can span multiple parent grids at
>> level L=N-1, but in the *union* of their parent grids they are
>> fully-contained.  If that's so, it mirrors how Orion/Boxlib operates,
>> and yt does support this behavior.
>
> Oh, I wasn't aware of this--so that's good.
>
>> I believe it may be as simple as converting Parent on the AMRGridPatch
>> object to be a list rather than a single object.  The Orion frontend
>> manages this, although in my "boxlib" bookmark I have cleaned up the
>> logic somewhat in yt/frontends/boxlib/ to make it easier to follow how
>> to do this.
>>
>> Does that sound workable?
>
> Yes, it does. I can have a look and work on it if you point me in the right direction. By the way, what about the parent/child relationship construction? Where is it done, if not in the frontend? Is there some kind of fallback somewhere?

Yup -- there's a function in the boxlib structure called
_reconstruct_parent_child:

http://paste.yt-project.org/show/3801/

which should probably be a base class method.

>
> Best,
>
> John
>
> _______________________________________________
> yt-dev mailing list
> yt-dev at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org



More information about the yt-dev mailing list