[yt-users] cutting plane not working or changed

[Slavin, Jonathan]

​Hi Nathan,

It turns out to be a bit different from what I said.  I tried with somewhat
different parameters so as to make the frb square and it worked.  So the
issue seems to be with making non-square cuts, or maybe in making the frb
from the​ cut.  I wonder if this is a recurrence of issue #1023, but it
doesn't look like it.  I don't have time at the moment to post my dataset,
but I will soon.  It is fairly large so I hope that's not a problem..

Jon

On Fri, Mar 10, 2017 at 4:05 PM, <yt-users-request at lists.spacepope.org>
wrote:

> Date: Fri, 10 Mar 2017 12:14:55 -0800
> From: Nathan Goldbaum <nathan12343 at gmail.com>
> To: Discussion of the yt analysis package
>         <yt-users at lists.spacepope.org>
> Subject: Re: [yt-users] cutting plane not working or changed behavior
> Message-ID:
>         <CAJXewOmWwY1CU4nrafG5_=GrAScGBnq=VMNxxh3s0OUzDzaGVw at mail.
> gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hi Jonathan,
>
> I'm not able to reproduce the behavior you're seeing with a test dataset
> (e.g. see http://paste.yt-project.org/show/7080/, this uses a dataset from
> yt-project.org/data).
>
> Is there any chance you can share a dataset that triggers this?
> Alternatively, can you reproduce this using one of the test datasets on
> yt-project.org/data?
>
> The yt curldrop might be useful to share a dataset that triggers the issue
> you're seeing, see docs.hub.yt. If you're not comfortable sharing the
> dataset publicly, feel free to share it with me off-list.
>
> -Nathan
>
> On Fri, Mar 10, 2017 at 12:02 PM, Slavin, Jonathan <
> jslavin at cfa.harvard.edu>
> wrote:
>
> > Hi all,
> >
> > I have some code that uses a cutting plane to create a non-axis aligned
> > slice plot:
> > ds = yt.load(file)
> > ctr = YTArray([15.,15.,15.],'pc')
> > cut = ds.cutting(normal=(1.,-1.,0.), center=ctr, north_vector=(0.,0.,1.))
> > frb = cut.to_frb(width=(42.427,'pc'),height=(30.,'pc'),resolution=
> > (512,724))
> > image = frb['density']/2.126E-24
> >
> > I then view the image using the matplotlib imshow command.
> > The data is a 30 pc**3 cartesian FLASH AMR dataset with maximum
> resolution
> > of
> > ?512**3.  I used this code to produce an image, which I have, so I know
> it
> > used to work (last April).  Now however it does not seem to project in
> the
> > same direction and seems to be shaped incorrectly with a top portion of
> the
> > image filled with NaNs.  Has anyone else experienced this behavior?  I
> > tried various things, for example reducing the resolution or switching
> the
> > order of those values (that leads to a core dump!).  The goal here is to
> > create a slice along the plane aligned with the z axis and 45 degrees
> away
> > from the x and y axes (the data goes from 0 to 30 pc in each direction).
> > Any help would be appreciated.?
> >
> > Jon
> >
> > --
> > ________________________________________________________
> > Jonathan D. Slavin                 Harvard-Smithsonian CfA
> > jslavin at cfa.harvard.edu       60 Garden Street, MS 83
> > phone: (617) 496-7981       Cambridge, MA 02138-1516
> > cell: (781) 363-0035             USA
> > ________________________________________________________
>




-- 
________________________________________________________
Jonathan D. Slavin                 Harvard-Smithsonian CfA
jslavin at cfa.harvard.edu       60 Garden Street, MS 83
phone: (617) 496-7981       Cambridge, MA 02138-1516
cell: (781) 363-0035             USA
________________________________________________________
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20170310/b0bbf372/attachment-0001.htm>