Dear all, I get the following error when trying to plot a radial profile of the density. Does anyone have an idea how to solve this ? Thanks a lot, Jean-Claude ##################################################################################################################### In [42]: ls CommonEnvelope0001 CommonEnvelope0001.configure CommonEnvelope0001.harrays CommonEnvelope0001.taskmap0000 CommonEnvelope0001.boundary CommonEnvelope0001.cpu0000 CommonEnvelope0001.hierarchy CommonEnvelope0001.taskmap0001 CommonEnvelope0001.boundary.hdf CommonEnvelope0001.cpu0001 CommonEnvelope0001.memorymap CommonEnvelope0001.yt In [43]: import yt.raven as R In [44]: import yt.lagos as L In [45]: filen = 'CommonEnvelope0001' In [46]: data = L.EnzoStaticOutput(filen) In [47]: plot_data = R.PlotCollection(data) yt INFO 2010-04-06 15:39:02,963 Created plot collection with default plot-center = [0.4921875, 0.4921875, 0.4921875] In [48]: plot_data.add_profile_sphere(0.5,'1',["Radius", "Density"]) yt INFO 2010-04-06 15:39:56,695 Getting field Radius from 1 yt INFO 2010-04-06 15:39:56,696 Getting field x from 1 yt INFO 2010-04-06 15:39:56,857 Getting field y from 1 yt INFO 2010-04-06 15:39:56,919 Getting field z from 1 yt INFO 2010-04-06 15:39:56,987 Getting field Density from 1 Out[48]: <yt.raven.PlotTypes.Profile1DPlot object at 0x3724d90> In [49]: plot_data.save('DensityProfile', format = 'png') --------------------------------------------------------------------------- OverflowError Traceback (most recent call last) /rpod2/jcpassy/Enzo/enzo/src/yt/scripts/iyt in <module>() ----> 1 2 3 4 5 /rpod2/jcpassy/Enzo/enzo/src/yt/yt/raven/PlotCollection.pyc in save(self, basename, format, override, force_save) 79 fn.append(plot.save_image(basename, \ 80 format=format, submit=self._run_id, ---> 81 override=override, force_save=force_save)) 82 if self.submit: 83 im = plot.im.copy() /rpod2/jcpassy/Enzo/enzo/src/yt/yt/raven/PlotTypes.pyc in save_image(self, prefix, format, submit, override, force_save) 110 *override* will force no filename generation beyond the prefix. 111 """ --> 112 self._redraw_image() 113 if not override: 114 self._generate_prefix(prefix) /rpod2/jcpassy/Enzo/enzo/src/yt/yt/raven/PlotTypes.pyc in _redraw_image(self) 711 func(self.data[self.fields[0]][indices], 712 self.data[self.fields[1]][indices], --> 713 **self.plot_options) 714 self._autoset_label(self.fields[0], self.set_x_label, 'x') 715 self._autoset_label(self.fields[1], self.set_y_label, 'y') /rpod2/jcpassy/yt-x86_64/lib/python2.6/site-packages/matplotlib/axes.pyc in loglog(self, *args, **kwargs) 3566 b = self._hold 3567 self._hold = True # we've already processed the hold -> 3568 l = self.plot(*args, **kwargs) 3569 self._hold = b # restore the hold 3570 /rpod2/jcpassy/yt-x86_64/lib/python2.6/site-packages/matplotlib/axes.pyc in plot(self, *args, **kwargs) 3435 3436 -> 3437 self.autoscale_view(scalex=scalex, scaley=scaley) 3438 return lines 3439 /rpod2/jcpassy/yt-x86_64/lib/python2.6/site-packages/matplotlib/axes.pyc in autoscale_view(self, tight, scalex, scaley) 1625 self.set_xbound(XL) 1626 if scaley and self._autoscaleYon: -> 1627 YL = self.yaxis.get_major_locator().view_limits(y0, y1) 1628 self.set_ybound(YL) 1629 /rpod2/jcpassy/yt-x86_64/lib/python2.6/site-packages/matplotlib/ticker.pyc in view_limits(self, vmin, vmax) 1117 vmin = minpos 1118 -> 1119 if not is_decade(vmin,self._base): vmin = decade_down(vmin,self._base) 1120 if not is_decade(vmax,self._base): vmax = decade_up(vmax,self._base) 1121 /rpod2/jcpassy/yt-x86_64/lib/python2.6/site-packages/matplotlib/ticker.pyc in is_decade(x, base) 1028 def is_decade(x,base=10): 1029 lx = math.log(x)/math.log(base) -> 1030 return lx==int(lx) 1031 1032 class LogLocator(Locator): OverflowError: cannot convert float infinity to integer #####################################################################################################################