<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="ltr">Sorry, forgot the output:<br><br>level # grids # cells # cells^3<br>----------------------------------------------<br> 0 1 262144 64<br> 1 1 24576 30<br> 2 1 38720 34<br> 3 1 98304 47<br> 4 1 147456 53<br> 5 1 327680 69<br> 6 1 269120 65<br> 7 1 403200 74<br> 8 1 645504 87<br>----------------------------------------------<br> 9 2216704<br><br><br>t = 1.06267789e-01 = 1.14972625e+12 s = 3.64326264e+04 years<br><br>Smallest Cell:<br> Width: 9.890e-12 Mpc<br> Width: 9.890e-06 pc<br> Width: 2.040e+00 AU<br> Width: 3.052e+13 cm<br><br></div><div class="gmail_extra"><br><div class="gmail_quote">2015-01-17 15:16 GMT-05:00 Yang Haifeng <span dir="ltr"><<a href="mailto:hfyangpku@gmail.com" target="_blank">hfyangpku@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"><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 class="h5"><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc 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:0 0 0 .8ex;border-left:1px #ccc 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:0 0 0 .8ex;border-left:1px #ccc 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>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 class="HOEnZb"><font color="#888888"><br><br clear="all"><br>-- <br><div><div dir="ltr">Best,<br>Haifeng<br></div></div>
</font></span></div>
</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr">Best,<br>Haifeng<br></div></div>
</div>