[yt-users] Fwd: cylindrical coordinates in Athena

Matthew Turk matthewturk at gmail.com
Thu Nov 15 14:50:50 PST 2012


Hi Kaitlin and John,

On Thu, Nov 15, 2012 at 5:49 PM, John ZuHone <jzuhone at gmail.com> wrote:
> Hi Kaitlin,
>
> At the moment yt supports 2D and 3D cartesian data for Athena (the former only very recently in the yt development tree).
>
> True cylindrical coordinate support (for all codes) is planned for yt 3.0. FLASH 2D cylindrical reads in fine, but this is mostly by accident since the data format in the file is the same as 2D cartesian. I'm guessing this means that the Athena vtk files for cylindrical coordinates are set up slightly differently from the cartesian coordinates.

As a quick followup, I *just now* blogged about 3.0 which includes
some discussion of non-cartesian coordinates:

http://blog.yt-project.org/post/WhatsUpWith30.html

I think John's idea is a good one, about using the existing code in
2.X as a stepping stone to swapping out the coordinate system in 3.0.

-Matt

>
> Are you using the stable or development version of yt? We might be able to get something working along the same lines as FLASH if you have a file you wouldn't mind providing, as a stepping-stone to true support for cylindrical coordinates.
>
> Best,
>
> John
>
> On Nov 15, 2012, at 5:36 PM, Kaitlin Kratter wrote:
>
>>
>>
>> Hi,
>>
>> I gather that cylindrical geometries aren't supported. Is it possible to force yt to read in cylindrical coordinates as though they were cartesian?
>>
>>
>> When I try to load an athena vtk in cylindrical, I get:
>>
>> In [32]: pf = load("cyl_disk.0001.vtk")
>> yt : [INFO     ] 2012-11-15 15:31:26,443 Temporarily setting domain_right_edge = -domain_left_edge. This will be corrected automatically if it is not the case.
>> yt : [INFO     ] 2012-11-15 15:31:26,445 Parameters: current_time              = 0.001080635
>> yt : [INFO     ] 2012-11-15 15:31:26,445 Parameters: domain_dimensions         = [        -37         128 -2147483648]
>> yt : [INFO     ] 2012-11-15 15:31:26,445 Parameters: domain_left_edge          = [ 1.       -1.570796  0.      ]
>> yt : [INFO     ] 2012-11-15 15:31:26,445 Parameters: domain_right_edge         = [-1.        1.570796 -0.      ]
>> yt : [INFO     ] 2012-11-15 15:31:26,445 Parameters: cosmological_simulation   = 0.0
>>
>>
>> This is a 256 by 128 grid, and the dimensions should be -pi/2, pi/2 and 1, 15.
>>
>>
>> Thanks,
>> Kaitlin
>>
>>
>>
>> _______________________________________________
>> 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



More information about the yt-users mailing list