<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="ltr">Hey Joseph,<div><br></div><div>I don't know if this will work, but maybe try (sorry about the weird indenting...trying to type while also watching the kids in bath)</div><div><br></div><div><ol style="margin:0px;padding:0px 0px 0px 48px;background-color:rgb(248,248,248)"><li class="" style="color:rgb(172,172,172);font-family:Consolas,Menlo,Monaco,'Lucida Console','Liberation Mono','DejaVu Sans Mono','Bitstream Vera Sans Mono',monospace,serif;font-size:12px;line-height:21px"><div class="" style="padding:0px 5px;vertical-align:top;color:rgb(0,0,0);border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);background:rgb(255,255,255)">ds <span class="" style="color:rgb(102,204,102)">=</span> load<span class="">(</span><span class="" style="color:rgb(72,61,139)">'myfile.hdf5'</span><span class="" style="color:rgb(102,204,102)">,</span>over_refine_factor<span class="" style="color:rgb(102,204,102)">=</span><span class="" style="color:rgb(255,69,0)">1</span><span class="">)</span></div></li><li class="" style=""><div class="" style="padding:0px 5px;vertical-align:top;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);background-image:initial;background-color:rgb(255,255,255);background-repeat:initial"><span class="" style="font-size:12px;line-height:21px"><font color="#000000" face="Consolas, Menlo, Monaco, Lucida Console, Liberation Mono, DejaVu Sans Mono, Bitstream Vera Sans Mono, monospace, serif"><div class="" style="padding:0px 5px;vertical-align:top;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);background-image:initial;background-repeat:initial"> bbox_lim = 50</div><div class="" style="padding:0px 5px;vertical-align:top;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);background-image:initial;background-repeat:initial"> center = whatever_your_desired_center</div><div class="" style="padding:0px 5px;vertical-align:top;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);background-image:initial;background-repeat:initial"> </div><div class="" style="padding:0px 5px;vertical-align:top;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);background-image:initial;background-repeat:initial"> bbox1 = [[center[0]-bbox_lim,center[0]+bbox_lim],</div><div class="" style="padding:0px 5px;vertical-align:top;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);background-image:initial;background-repeat:initial"> [center[1]-bbox_lim,center[1]+bbox_lim],</div><div class="" style="padding:0px 5px;vertical-align:top;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);background-image:initial;background-repeat:initial"> [center[2]-bbox_lim,center[2]+bbox_lim]]</div></font></span></div></li><li class="" style=""><div class="" style="padding:0px 5px;vertical-align:top;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);background-image:initial;background-repeat:initial"><div class="" style="padding:0px 5px;vertical-align:top;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);background-image:initial;background-repeat:initial"> ds1 = load(fname,bounding_box=bbox1,n_ref = 64 ,over_refine_factor=1)<br></div></div></li></ol></div><div class="gmail_extra"><br></div><div class="gmail_extra">>ad1 = ds1.all_data()</div><div class="gmail_extra">and search within ad1?</div><div class="gmail_extra"><br></div><div class="gmail_extra">you might need to force periodicity off</div><div class="gmail_extra">> ds1.periodicity = (False,False,False)</div><div><br></div><div>but maybe not. </div><div><br></div><div>Anyways, I hope this works!</div><div>-d</div><div class="gmail_extra"><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, May 2, 2015 at 6:27 PM, Joseph Smidt <span dir="ltr"><<a href="mailto:josephsmidt@gmail.com" target="_blank">josephsmidt@gmail.com</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">Hey everyone,<br>
<br>
I am trying to find the densest point in Gadget data in the<br>
subvolume with left edge [700]*3 and right edge [800]*3 where 700 and<br>
800 are in code units.<br>
<br>
This is what I do [1], which was inspired by this page [2],<br>
except I try and take the all_data confined to a region. The output<br>
[3] says however where you will note the center is not inside the<br>
region. It found the densest point in the whole data set, not the<br>
region.<br>
<br>
So how would I do this instead for a subvolume of Gadget data? Thanks.<br>
<br>
<br>
[1] <a href="http://pastebin.com/kyUgEi3a" target="_blank">http://pastebin.com/kyUgEi3a</a><br>
<br>
[2] <a href="http://yt-project.org/doc/cookbook/gadget_notebook.html" target="_blank">http://yt-project.org/doc/cookbook/gadget_notebook.html</a><br>
<br>
[3] <a href="http://pastebin.com/E4aLb5VC" target="_blank">http://pastebin.com/E4aLb5VC</a><br>
<br>
<br>
<br>
--<br>
------------------------------------------------------------------------<br>
Joseph Smidt <<a href="mailto:josephsmidt@gmail.com">josephsmidt@gmail.com</a>><br>
<br>
Theoretical Division<br>
P.O. Box 1663, Mail Stop B283<br>
Los Alamos, NM 87545<br>
Office: <a href="tel:505-665-9752" value="+15056659752">505-665-9752</a><br>
Fax: <a href="tel:505-667-1931" value="+15056671931">505-667-1931</a><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>
</blockquote></div><br></div></div>