<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="ltr">Thanks! <br></div><div class="gmail_extra"><br><div class="gmail_quote">2015-01-17 15:47 GMT-05:00 Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Sat, Jan 17, 2015 at 12:16 PM, Yang Haifeng <span dir="ltr"><<a href="mailto:hfyangpku@gmail.com" target="_blank">hfyangpku@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"><div dir="ltr"><div>Hi Nathan,<br><br></div>Here is the output of ds.print_stats(). I think it hasn't run off the edge yet. The refine ratio is 2 and top grid is 64*3. ad_c["Density"] will not produce any error.<br></div><div class="gmail_extra"><br><div class="gmail_quote">2015-01-17 12:48 GMT-05:00 Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a>></span>:<div><div><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 Haifeng,<br><div class="gmail_extra"><br><div class="gmail_quote"><span>On Sat, Jan 17, 2015 at 8:55 AM, Yang Haifeng <span dir="ltr"><<a href="mailto:hy4px@virginia.edu" target="_blank">hy4px@virginia.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"><div><div>Dear yt users,<br><br></div>I am trying to make covering_grid from an enzo simulation. The code is very simple: (somewhat sloppy since it's running with ipython. )<br><br>ds=yt.load("data0201")<br></div></div></blockquote><div><br></div></span><div>What is the output of "ds.print_stats()" for this dataset? I'm wondering if your covering grid is running off the edge of the simulation domain.</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><div><div dir="ltr"><div>ad_c=ds.covering_grid(level=3, left_edge=0.4375, dims=[64,64,64])<br></div></div></div></div></blockquote></div></div></div></blockquote></div></div></div></div></blockquote><div><br></div></span><div>Ah, I see what's wrong now. This should be left_edge=[0.4375, 0.4375, 0.4375] (or some other 3-tuple).</div><div><br></div><div>I've gone ahead and added some basic input handling code for the covering_grid constructor, so future users will get a less confusing error message:</div><div><br></div><div><a href="https://bitbucket.org/yt_analysis/yt/pull-request/1414/add-some-input-sanitizing-and-error/diff" target="_blank">https://bitbucket.org/yt_analysis/yt/pull-request/1414/add-some-input-sanitizing-and-error/diff</a><br></div><div><div class="h5"><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 class="gmail_extra"><div class="gmail_quote"><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 class="gmail_extra"><div class="gmail_quote"><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><div><div dir="ltr"><div></div>ad_c["x"]<br clear="all"><div><div><div><div><br></div><div>where the ad_c["x"] will give an Index Error on my desktop. However, the same piece of code runs great in a different machine. The version of yt are the same in both machines, which is 3.0.2. Both of them installed yt through anaconda. The error message I got was following. Thanks for your help!<br><br>---------------------------------------------------------------------------<br>IndexError Traceback (most recent call last)<br><ipython-input-6-ff78a8bccc54> in <module>()<br>----> 1 ad_c["x"]<br><br>/home/fxyang/software/anaconda/lib/python2.7/site-packages/yt/data_objects/data_containers.pyc in __getitem__(self, key)<br> 243 if f in self._container_fields:<br> 244 self.field_data[f] = \<br>--> 245 self.ds.arr(self._generate_container_field(f))<br> 246 return self.field_data[f]<br> 247 else:<br><br>/home/fxyang/software/anaconda/lib/python2.7/site-packages/yt/data_objects/construction_data_containers.pyc in _generate_container_field(self, field)<br> 587 np.multiply(rv, self.dds[2], rv)<br> 588 elif field == ("index", "x"):<br>--> 589 x = np.mgrid[self.left_edge[0] + 0.5*self.dds[0]:<br> 590 self.right_edge[0] - 0.5*self.dds[0]:<br> 591 self.ActiveDimensions[0] * 1j]<br><br>/home/fxyang/software/anaconda/lib/python2.7/site-packages/yt/units/yt_array.pyc in __getitem__(self, item)<br> 981 <br> 982 def __getitem__(self, item):<br>--> 983 ret = super(YTArray, self).__getitem__(item)<br> 984 if ret.shape == ():<br> 985 return YTQuantity(ret, self.units)<br><br>IndexError: too many indices for array<span><font color="#888888"><br><br></font></span></div><span><font color="#888888"><div><br>-- <br><div><div dir="ltr">Best,<br>Haifeng<br></div></div>
<br>
</div></font></span></div></div></div></div>
<br></div></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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div></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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div></div></div><span><font color="#888888"><br><br clear="all"><br>-- <br><div><div dir="ltr">Best,<br>Haifeng<br></div></div>
</font></span></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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div></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><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr">Best,<br>Haifeng<br></div></div>
</div>