<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Mar 10, 2017 at 1:50 PM, Slavin, Jonathan <span dir="ltr"><<a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div style="font-family:arial,helvetica,sans-serif;font-size:small">Hi Nathan,</div><div style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div 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></div></blockquote><div><br></div><div>No worries, I think the limit on the curldrop is 10 GB.</div><div><br></div><div>In my test script I was making a non-square cutting plane and wasn't seeing NaNs, but it's also possible that I wasn't triggering the failure mode you're seeing. Unfortunately due to the nature of yt working with opaque data on-disk usually it's sometimes very difficult to figure out how to reproduce an issue without a test dataset to work with (even if the issue isn't specific to that dataset).</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div 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.spacep<wbr>ope.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_4685417845547668254m_5531975511922015978:191" class="m_4685417845547668254m_5531975511922015978a3s m_4685417845547668254m_5531975511922015978aXjCH m_4685417845547668254m_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.gm<wbr>ail.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><span class="HOEnZb"><font color="#888888"><br><br><br clear="all"><div><br></div>-- <br><div class="m_4685417845547668254m_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>
</font></span></div></div>
<br>______________________________<wbr>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
<br></blockquote></div><br></div></div>