<div dir="ltr"><div class="gmail_extra"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">​Hi Nathan,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">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..</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Jon</div><br><div class="gmail_quote">On Fri, Mar 10, 2017 at 4:05 PM,  <span dir="ltr"><<a href="mailto:yt-users-request@lists.spacepope.org" target="_blank">yt-users-request@lists.<wbr>spacepope.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id="m_5531975511922015978:191" class="m_5531975511922015978a3s m_5531975511922015978aXjCH m_5531975511922015978m15aba0a525517a43">Date: Fri, 10 Mar 2017 12:14:55 -0800<br>
From: Nathan Goldbaum <<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a>><br>
To: Discussion of the yt analysis package<br>
        <<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a>><br>
Subject: Re: [yt-users] cutting plane not working or changed behavior<br>
Message-ID:<br>
        <CAJXewOmWwY1CU4nrafG5_=GrAScG<wbr>Bnq=<a href="mailto:VMNxxh3s0OUzDzaGVw@mail.gmail.com" target="_blank">VMNxxh3s0OUzDzaGVw@mail.<wbr>gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi Jonathan,<br>
<br>
I'm not able to reproduce the behavior you're seeing with a test dataset<br>
(e.g. see <a href="http://paste.yt-project.org/show/7080/" rel="noreferrer" target="_blank">http://paste.yt-project.org/sh<wbr>ow/7080/</a>, this uses a dataset from<br>
<a href="http://yt-project.org/data" rel="noreferrer" target="_blank">yt-project.org/data</a>).<br>
<br>
Is there any chance you can share a dataset that triggers this?<br>
Alternatively, can you reproduce this using one of the test datasets on<br>
<a href="http://yt-project.org/data" rel="noreferrer" target="_blank">yt-project.org/data</a>?<br>
<br>
The yt curldrop might be useful to share a dataset that triggers the issue<br>
you're seeing, see <a href="http://docs.hub.yt" rel="noreferrer" target="_blank">docs.hub.yt</a>. If you're not comfortable sharing the<br>
dataset publicly, feel free to share it with me off-list.<br>
<br>
-Nathan<br>
<br>
On Fri, Mar 10, 2017 at 12:02 PM, Slavin, Jonathan <<a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a>><br>
wrote:<br>
<br>
> Hi all,<br>
><br>
> I have some code that uses a cutting plane to create a non-axis aligned<br>
> slice plot:<br>
> ds = yt.load(file)<br>
> ctr = YTArray([15.,15.,15.],'pc')<br>
> cut = ds.cutting(normal=(1.,-1.,0.), center=ctr, north_vector=(0.,0.,1.))<br>
> frb = cut.to_frb(width=(42.427,'pc')<wbr>,height=(30.,'pc'),resolution=<br>
> (512,724))<br>
> image = frb['density']/2.126E-24<br>
><br>
> I then view the image using the matplotlib imshow command.<br>
> The data is a 30 pc**3 cartesian FLASH AMR dataset with maximum resolution<br>
> of<br>
> ?512**3.  I used this code to produce an image, which I have, so I know it<br>
> used to work (last April).  Now however it does not seem to project in the<br>
> same direction and seems to be shaped incorrectly with a top portion of the<br>
> image filled with NaNs.  Has anyone else experienced this behavior?  I<br>
> tried various things, for example reducing the resolution or switching the<br>
> order of those values (that leads to a core dump!).  The goal here is to<br>
> create a slice along the plane aligned with the z axis and 45 degrees away<br>
> from the x and y axes (the data goes from 0 to 30 pc in each direction).<br>
> Any help would be appreciated.?<br>
><br>
> Jon<br>
><br>
> --<br>
> ______________________________<wbr>__________________________<br>
> Jonathan D. Slavin                 Harvard-Smithsonian CfA<br>
> <a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a>       60 Garden Street, MS 83<br>
> phone: <a href="tel:%28617%29%20496-7981" value="+16174967981" target="_blank">(617) 496-7981</a>       Cambridge, MA 02138-1516<br>
> cell: <a href="tel:%28781%29%20363-0035" value="+17813630035" target="_blank">(781) 363-0035</a>             USA<br>
> ______________________________<wbr>__________________________</div></blockquote></div><br><br><br clear="all"><div><br></div>-- <br><div class="m_5531975511922015978gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">______________________________<wbr>__________________________<br>Jonathan D. Slavin                 Harvard-Smithsonian CfA<br><a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a>       60 Garden Street, MS 83<br>phone: <a href="tel:(617)%20496-7981" value="+16174967981" target="_blank">(617) 496-7981</a>       Cambridge, MA 02138-1516<br>cell: <a href="tel:(781)%20363-0035" value="+17813630035" target="_blank">(781) 363-0035</a>             USA<br>______________________________<wbr>__________________________<br><br></div></div></div></div>
</div></div>