<div dir="ltr">Hi John,<div><br></div><div>Yes, you should be able to do this. I think the clumps have an attribute that describes their value, and then you can supply that to the surface operator. I've gotten your bug report and will be looking into supplying clumps to surfaces.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jul 25, 2016 at 3:08 PM, John Regan <span dir="ltr"><<a href="mailto:johnanthonyregan@gmail.com" target="_blank">johnanthonyregan@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><p dir="ltr">Hi Matt,</p>
<p dir="ltr">Thanks for looking at this. Much appreciated! I'll file a bug report no problem. </p>
<p dir="ltr">For calling the surface operator with a clump as a data source is it possible to extract the isodensity value of the clump from the Clump object and pass that to the surface operator? </p>
<p dir="ltr">Cheers,<br>
John </p><div class="HOEnZb"><div class="h5">
<div class="gmail_extra"><br><div class="gmail_quote">On 25 Jul 2016 19:58, "Matthew Turk" <<a href="mailto:matthewturk@gmail.com" target="_blank">matthewturk@gmail.com</a>> wrote:<br type="attribution"><blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi John,<div><br></div><div>Sadly, we don't store information about which cells are on the boundary. This is potentially possible to compute (although the naive solution, which is to check the # of faces of a zone that have neighbors in the clump, will not work for clumps that have holes in them) but it's not stored just now. I don't have an alternate solution to present right now to that problem.</div><div><br></div><div>For the second part, I'm not sure what's going on. Could you file a bug with a repro script?</div><div><br></div><div>-Matt</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div>On Mon, Jul 25, 2016 at 8:33 AM, John Regan <span dir="ltr"><<a href="mailto:johnanthonyregan@gmail.com" target="_blank">johnanthonyregan@gmail.com</a>></span> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div dir="ltr"><div><div><div><div><div><div><div><div>Hi All, <br><br></div>I'm interested in calculating some surface properties of clumps found with the clump finder. Finding the clumps is straight forward via the clump finder API. Once I find a clump then I apply the surface operator<br><br>e.g.<br>surf = clump.data.ds.surface(clump, "density", 1e-24)<br><br></div>This works but isn't ideal. For one I don't want to specify a value for the density contour as by definition this is part of the clump - ideally I'd just be able to extract the surface cells from the clump with something like <br><br></div>surf2 = clump.surface() <br><br></div>Does YT store any information on which cells reside on the surface of a clump? <br><br></div>Even pushing ahead with the YTSurface object doesn't go very far anyway as even trying to access the triangles variable gives the error<br><br>triangles = surf.triangles<br><br>File "/cosma/home/dp004/dc-rega4/YT/yt-x86_64/src/yt-hg/yt/analysis_modules/level_sets/virial_validator.py", line 74, in _virial_validation<br> triangles = surf.triangles<br> File "/cosma/home/dp004/dc-rega4/YT/yt-x86_64/src/yt-hg/yt/data_objects/construction_data_containers.py", line 1223, in triangles<br> self.get_data()<br> File "/cosma/home/dp004/dc-rega4/YT/yt-x86_64/src/yt-hg/yt/data_objects/construction_data_containers.py", line 1108, in get_data<br> for io_chunk in parallel_objects(self.data_source.chunks([], "io")):<br>AttributeError: 'Clump' object has no attribute 'chunks'<br><br><br></div>Any tips/ideas greatly appreciated!<br><br></div>Cheers,<br></div>John<br><br><div><div><div><div><div><div><div><div><div><br></div></div></div></div></div></div></div></div></div></div>
<br></div>_______________________________________________<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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>
<br>_______________________________________________<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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>
</div></div><br>_______________________________________________<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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>