[yt-users] Projection of 2D cylindrical FLASH data

[Suoqing JI]

One way (imperfect but definitely doable) you might want to try is to map 2D cylindrical data into a 3D numpy array, which could be loaded by yt.load_uniform_grid and passed into the volume rendering module.

See here: http://yt-project.org/doc/examining/generic_array_data.html

Best wishes,
--
Suoqing JI
Ph.D Student
Department of Physics
University of California, Santa Barbara
http://web.physics.ucsb.edu/~suoqing

> On May 14, 2015, at 6:01 PM, Nathan Goldbaum <nathan12343 at gmail.com> wrote:
> 
> 
> 
> On Thu, May 14, 2015 at 5:35 PM, Ken Shen <kenshen at astro.berkeley.edu <mailto:kenshen at astro.berkeley.edu>> wrote:
> I'm really sorry, I can't seem to explain this very clearly in words!
> 
> I have 2D axisymmetric data.  Let's suppose that I map this to a 3D Cartesian dataset.  I'd like to make a 2D plot from this 3D data that consists of integrating the density along rays from the observer through the dataset onto a plane on the other side of the data.  Ideally, I'd like to do this with the observer at an arbitrary position with respect to the axes.
> 
> I haven't done it, but it seems like this is built in to yt's functionality for 3D datasets.  But I haven't found how to do it for a 2D axisymmetric dataset.
> 
> 
> In yt we would call this a volume rendering using rays that are aligned with the cartesian axes and a projection transfer function.
> 
> I don't think this has been implemented yet, but volume renderings of data with non-cartesian geometries or with unstructured meshes is something we would like to do.
>  
> 
> Sorry for all the messages!
> Ken
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> On Thu, May 14, 2015 at 5:28 PM, Michael Zingale <michael.zingale at stonybrook.edu <mailto:michael.zingale at stonybrook.edu>> wrote:
> to be clear here, what Flash calls cylindrical is 2-d axisymmetric, r-z.  If I understand what you are asking, you basically want to still have a 2-d plot where each cell is the integral over an annular region as if that cell were rotated about the symmetry axis -- is that right?  
> 
> On Thu, May 14, 2015 at 8:21 PM, Ken Shen <kenshen at astro.berkeley.edu <mailto:kenshen at astro.berkeley.edu>> wrote:
> Hi, sorry I wasn't very clear before.  I'd like to make 2D column density plots of my 2D cylindrical data as if my data were 3D cylindrical data with no phi-dependence.
> 
> (The ProjectionPlot fails for axis = 1 or 3.  For axis = 2, it essentially recreates the regular density slice plot.)
> 
> 
> Thanks,
> Ken
> 
> 
> 
> 
> On Thu, May 14, 2015 at 5:00 PM, Matthew Turk <matthewturk at gmail.com <mailto:matthewturk at gmail.com>> wrote:
> Hi Ken,
> 
> (I think Suoqing's answer may work as well.)  If you want to project
> along r, z, or theta, you can do a ProjectionPlot (or "ds.proj") and
> specify "r", "z", or "theta", and the result will be along that axis
> with the correct path lengths.  I think that the "r" direction may be
> what you're looking for here, i.e., a soup-can label.
> 
> -Matt
> 
> On Thu, May 14, 2015 at 6:31 PM, Ken Shen <kenshen at astro.berkeley.edu <mailto:kenshen at astro.berkeley.edu>> wrote:
> > Hi all, I'd like to make column density plots of my 2D cylindrical FLASH
> > data.  Ideally, these would be at arbitrary angles to the z-axis.  I'm new
> > to yt, so I could very well have missed it, but I couldn't seem to find this
> > in the documentation.  Any hints?
> >
> >
> > Thanks,
> > Ken
> >
> > _______________________________________________
> > yt-users mailing list
> > yt-users at lists.spacepope.org <mailto:yt-users at lists.spacepope.org>
> > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org <http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org>
> >
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org <mailto:yt-users at lists.spacepope.org>
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org <http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org>
> 
> 
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org <mailto:yt-users at lists.spacepope.org>
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org <http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org>
> 
> 
> 
> 
> -- 
> Michael Zingale
> Associate Professor
> 
> Dept. of Physics & Astronomy • Stony Brook University • Stony Brook, NY 11794-3800
> phone:  631-632-8225 <tel:631-632-8225>
> e-mail: Michael.Zingale at stonybrook.edu <mailto:Michael.Zingale at stonybrook.edu>
> web: http://www.astro.sunysb.edu/mzingale <http://www.astro.sunysb.edu/mzingale>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org <mailto:yt-users at lists.spacepope.org>
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org <http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org>
> 
> 
> 
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org <mailto:yt-users at lists.spacepope.org>
> http://lists.spacepope.org/listinfo.cgi/yt-users-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/20150514/09ccc5f4/attachment.htm>
-------------- next part --------------
_______________________________________________
yt-users mailing list
yt-users at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org