<div dir="ltr">I've issued a pull request which I hope will fix the issues that were brought up in this thread:<div><br></div><div><a href="https://github.com/yt-project/yt/pull/1409">https://github.com/yt-project/yt/pull/1409</a></div><div><br></div><div>Kazem, if you apply that pull request locally then I think it should fix your issues. Please let us know if you need a hand compiling a version of yt from git with that fix applied.<br><div><br></div><div>One last note: the isocontour algorithm yt uses, marching cubes, needs data that bracket the chosen isocontour value. That means that if you create a sphere data object and then ask yt to create an isocontour from that sphere data object at the radius of the sphere, you will not get the results you expect. To generate a valid isosurface, yt needs field values both above and below the field value, but the sphere data object will only return values smaller than its own radius.</div><div><br></div><div>In practice, that means you shouldn't generate isocontours like this:</div><div><br></div><div> sp = ds.sphere(ds.domain_center, (0.5, 'code_length')</div></div><div> surface = ds.surface(sp, 'radius', (0.5, 'code_length'))</div><div><br></div><div>But instead should do something like this:</div><div><br></div><div> ad = ds.all_data()</div><div> surface = ds.surface(ad, 'radius', (0.5, 'code_length'))</div><div><br></div><div>Note that I didn't necessarily need to use `ds.all_data()` there, I could have used a sphere with a sufficiently large radius as to capture all of the bracketing radius values we need to generate the isosurface correctly. Using ds.all_data() just ensures that we will select all of the data we need no matter what.</div><div><br></div><div>Unfortunately I'm not sure there's a way yt can detect that you are in this situation and print a useful error or warning.</div><div><br></div><div>I hope that clears things up :)</div><div><br></div><div>-Nathan</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 19, 2017 at 9:53 AM, Britton Smith <span dir="ltr"><<a href="mailto:brittonsmith@gmail.com" target="_blank">brittonsmith@gmail.com</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">Hi Kazem, Nathan,<div><br></div><div>I've opened two separate issues for this, the one Nathan linked to above and this one:</div><div><a href="https://github.com/yt-project/yt/issues/1407" target="_blank">https://github.com/yt-project/<wbr>yt/issues/1407</a><br></div><div><br></div><div>Issue #1406 is just that the surface value and surface fields are never verified to be in the same units. That should be trivial to fix. Issue #1407 gets at the real problem, which is that surfaces made with the radius field are almost always empty. I say almost because as I was trying to recreate this issue this morning, I happened to make one that had some values. However, if I twiddle things just a little, the size is again zero, which seems incorrect to me.</div><div><br></div><div>Anyway, keep an eye on those issues and we will post updates there as we're able.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Britton</div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 19, 2017 at 7:38 AM, Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</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">Ah, it looks like Britton filed an issue while I was writing that e-mail:<div><br></div><div><a href="https://github.com/yt-project/yt/issues/1406" target="_blank">https://github.com/yt-project/<wbr>yt/issues/1406</a><br></div></div><div class="m_4887521608690275410HOEnZb"><div class="m_4887521608690275410h5"><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 19, 2017 at 9:36 AM, Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</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">Hi Kazem,<div><br></div><div>Sorry, I'm not quite sure what you were asking, I haven't actually done anything to help out here yet. As Britton said there is likely a bug that needs to be fixed. I asked him for a runnable test script I could use to debug the issue. If you feel like it, you could create a runnable test script as well and then file an issue on github:</div><div><br></div><div><a href="https://github.com/yt-project/yt/issues/new" target="_blank">https://github.com/yt-project/<wbr>yt/issues/new</a><br></div><div><br></div><div>There are a number of public test datasets you can use in your script at <a href="http://yt-project.org/data" target="_blank">yt-project.org/data</a>. Using a public dataset is important because it ensures anyone will be able to reproduce the issue you're seeing (sometimes problems depend on the exact type of data you are looking at).</div><div><br></div><div>Hopefully once we have an example to debug with one of us will be able to come up with a bugfix you can apply to your local yt installation and that others can use once the next yt release goes out.</div><span class="m_4887521608690275410m_3320620644062979149HOEnZb"><font color="#888888"><div><br></div><div>-Nathan</div></font></span></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_4887521608690275410m_3320620644062979149h5">On Fri, May 19, 2017 at 7:25 AM, Kazem Ardaneh <span dir="ltr"><<a href="mailto:kazem.arrdaneh@gmail.com" target="_blank">kazem.arrdaneh@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="m_4887521608690275410m_3320620644062979149h5"><div style="word-wrap:break-word"><i><b>Hi <span style="white-space:pre-wrap;background-color:rgb(255,255,255)">Nathan, </span></b></i><div><span style="white-space:pre-wrap;background-color:rgb(255,255,255)"><i><b><br></b></i></span></div><div><span style="white-space:pre-wrap;background-color:rgb(255,255,255)"><i><b>Could you kindly let me know how do you define:</b></i></span></div><div><span style="white-space:pre-wrap;background-color:rgb(255,255,255)"><i><b><br></b></i></span></div><div><span style="white-space:pre-wrap;background-color:rgb(255,255,255)"><i><b>1. sphere (sph object),</b></i></span></div><div><i><b><span style="white-space:pre-wrap;background-color:rgb(255,255,255)">2. </span><span style="white-space:pre-wrap;background-color:rgb(255,255,255)">surface (Surface structure)</span></b></i></div><div><span style="white-space:pre-wrap;background-color:rgb(255,255,255)"><i><b><br></b></i></span></div><div><span style="white-space:pre-wrap;background-color:rgb(255,255,255)"><i><b>It seems me you defined in code unit.</b></i></span></div><div><span style="white-space:pre-wrap;background-color:rgb(255,255,255)"><i><br></i></span></div><div><span style="white-space:pre-wrap;background-color:rgb(255,255,255)"><i><br></i></span><div>
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div>Best Regards,</div><div><br></div><div>kazem</div><div><br></div></div><br class="m_4887521608690275410m_3320620644062979149m_-8995136943301366197m_-4945421463629910545Apple-interchange-newline"><br class="m_4887521608690275410m_3320620644062979149m_-8995136943301366197m_-4945421463629910545Apple-interchange-newline">
</div>
<br></div></div><br></div></div><span>______________________________<wbr>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">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/lis<wbr>tinfo.cgi/yt-users-spacepope.o<wbr>rg</a><br>
<br></span></blockquote></div><br></div>
</blockquote></div><br></div>
</div></div><br>______________________________<wbr>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">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/lis<wbr>tinfo.cgi/yt-users-spacepope.<wbr>org</a><br>
<br></blockquote></div><br></div>
</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>