<div dir="ltr"><br><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Jun 24, 2014 at 8:57 AM, Desika Narayanan <span dir="ltr"><<a href="mailto:dnarayan@haverford.edu" target="_blank">dnarayan@haverford.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">Hi All,<div><br></div><div>I'm having some issues with really basic halo finding on gadget snapshots in the bleeding edge version. What I want to do is ID halos in a volume, save a list of the halos, and plot a projection plot with circles over the halos (all of which I think should be possible).</div>
<div><br></div><div>The principle script is: <a href="http://paste.yt-project.org/show/4824/" target="_blank">http://paste.yt-project.org/show/4824/</a></div><div><br></div><div>The error I get when trying to save the halos is:</div>
<div>
<br></div><div><a href="http://paste.yt-project.org/show/4822/" target="_blank">http://paste.yt-project.org/show/4822/</a></div></div></blockquote><div><br></div><div>I don't think I can speak to this beyond that this error indicates that something is going wrong with field detection. It might be helpful to rerun the script with yt's log level threshold decreased to 10, which will show debug messages. See this page: <a href="http://yt-project.org/docs/dev-3.0/reference/configuration.html">http://yt-project.org/docs/dev-3.0/reference/configuration.html</a></div>
<div><br></div><div>Running with debug logs will allow us to see if any fields are raising exceptions during field detection.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div dir="ltr"><div><br></div><div><br></div><div>and in trying to annotate the halos on the projection plot is:</div><div><br></div><div>
<a href="http://paste.yt-project.org/show/4823/" target="_blank">http://paste.yt-project.org/show/4823/</a><br></div><div><br></div><div>Am I either doing something wrong in my commands? (hopefully as simple as that...) or are these not supported yet in the bleeding edge yt?</div>
</div></blockquote><div><br></div><div>This looks like a straight up typo in the HaloCatalogCallback. You could try commenting out the line it's failing on.</div><div><br></div><div>It looks like Hillary Egan was the last one to touch this in May.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">
<div><br></div><div>thanks</div><span class=""><font color="#888888"><div>desika</div></font></span><div><br></div><div>ps. On another (related) note, is it possible to manually draw circles on a projection plot? Or arrows, or tick marks, or anything that I could use to visually show off individual halos? </div>
</div></blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div> </div></div></blockquote>
<div><br></div><div>You can access the underlying matplotlib primitives like so:</div><div><br></div><div>slc = SlicePlot(pf, 0, 'density')</div><div>plot = slc.plots['density']</div><div>fig = plot.figure # matplotlib figure object</div>
<div>ax = plot.axes # matplotlib axes object</div><div>cax = plot.cax # matplotlib colorbar axes object</div><div>cb = plot.cb # matplotlib colorbar object</div><div>image = plot.image # matplotlib AxesImage object</div>
<div><br></div><div>Using these objects you should be able to customize the plot however you'd like.</div><div><br></div><div>-Nathan</div><div><br></div><div><br></div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div dir="ltr">
<div><br></div><div><br></div><div><br></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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div></div>