Erm, *cough*, please ignore me completely. (why do I always figure things out after sending emails).<br><br>Sorry to bother you all.<br><br>Libby<br><br>ps - if youa re interested p.set_cmap("cold") was the offending part, doh!!<br>
<br><div class="gmail_quote">On 4 November 2010 12:04, Elizabeth Harper-Clark <span dir="ltr"><<a href="mailto:h-clark@astro.utoronto.ca">h-clark@astro.utoronto.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Ah, upon further testing the error seems to be related to plotting HII_fraction rather than 3 plots. However, I am able to plot HII fraction slices fine with the below script:<br><br>from yt.mods import *<br>import pylab<br>
import sys<br>import cmath<br>import matplotlib.ticker<br><br>rho_min = 1e-24<br>rho_max = 2e-24<br><br>ss = 1<div class="im"><br>centre = (0.511901, 0.452331, 0.434998)<br></div>frame_template = "aaSN/%s/ionisationframezoom_%01i_%04i"<br>
<br>for n in range(132,401,1): # make density projection with particles on<br> print "Dump:", n<div class="im"><br> pf = load("SciNet/%s/DD%04i/data%04i" % (sys.argv[-1],n, n))<br></div><div class="im">
pc = PlotCollection(pf, center=centre)<br></div>
Timemyr = pf["InitialTime"]*pf["years"]/1.e6<br> print "Time(code units):", pf["InitialTime"]<br> for ax in range(3):<br> p = pc.add_slice("HII_Fraction", ax)<div class="im">
<br>
p.set_width(ss, 'pc')<br></div> sp = pf.h.all_data()<br> if na.any(pf.h.grid_particle_count):<br> colours = sp["ParticleMassMsun"]<br> p.modify["particles"](ss, p_size=6.0, col='k', minimum_mass= 6e-26)<br>
print " XXXXX Particles present XXXXX"<br> p.modify["text"]((0.101,0.899),"Time = %8.4f Myr" % (Timemyr), text_args = {'color':'k'})<br> p.modify["text"]((0.1,0.9),"Time = %8.4f Myr" % (Timemyr), text_args = {'color':'w'})<br>
pc.save(frame_template % (sys.argv[-1],ax,n), override=True)<br><br>????<div><div></div><div class="h5"><br><br><div class="gmail_quote">On 4 November 2010 11:56, Elizabeth Harper-Clark <span dir="ltr"><<a href="mailto:h-clark@astro.utoronto.ca" target="_blank">h-clark@astro.utoronto.ca</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi all,<br><br>I have adapted the multi plot script from the cookbook an dit all works as long as I only want 2 plot. When I try to increase to three I get:<br>
<br>[h-clark@tpb4 my_scripts]$ python2.6 MultiPlot.py GMCsetup3<br>
1.0 0.714285714286<br>yt INFO 2010-11-04 11:53:30,115 Created plot collection with default plot-center = [0.51190100000000005, 0.45233099999999998, 0.434998]<br>yt INFO 2010-11-04 11:53:30,127 Getting the binary hierarchy<br>
yt INFO 2010-11-04 11:53:30,160 Finished with binary hierarchy reading<br>yt INFO 2010-11-04 11:53:30,253 Adding PotentialField to list of fields<br>yt INFO 2010-11-04 11:53:30,254 Adding AccelerationField1 to list of fields<br>
yt INFO 2010-11-04 11:53:30,254 Adding HeII_gamma to list of fields<br>yt INFO 2010-11-04 11:53:30,254 Adding RadAccel1 to list of fields<br>yt INFO 2010-11-04 11:53:30,254 Adding RadAccel2 to list of fields<br>
yt INFO 2010-11-04 11:53:30,255 Adding RadAccel3 to list of fields<br>yt INFO 2010-11-04 11:53:30,255 Adding HI_kph to list of fields<br>yt INFO 2010-11-04 11:53:30,255 Adding HeII_kph to list of fields<br>
yt INFO 2010-11-04 11:53:30,255 Adding HI_gamma to list of fields<br>[0 0 0]<br>yt INFO 2010-11-04 11:53:31,962 Added slice of Density at x = 0.511901 with 'center' = [0.51190100000000005, 0.45233099999999998, 0.434998]<br>
yt INFO 2010-11-04 11:53:33,265 Added slice of Temperature at x = 0.511901 with 'center' = [0.51190100000000005, 0.45233099999999998, 0.434998]<br>yt INFO 2010-11-04 11:53:34,933 Added slice of HII_Fraction at x = 0.511901 with 'center' = [0.51190100000000005, 0.45233099999999998, 0.434998]<br>
Traceback (most recent call last):<br> File "MultiPlot.py", line 34, in <module><br> fig.savefig(frame_template % (sys.argv[-1],ax,n))<br> File "/home/h-clark/YT/lib/python2.6/site-packages/matplotlib/figure.py", line 1032, in savefig<br>
self.canvas.print_figure(*args, **kwargs)<br> File "/home/h-clark/YT/lib/python2.6/site-packages/matplotlib/backend_bases.py", line 1476, in print_figure<br> **kwargs)<br> File "/home/h-clark/YT/lib/python2.6/site-packages/matplotlib/backends/backend_agg.py", line 358, in print_png<br>
FigureCanvasAgg.draw(self)<br> File "/home/h-clark/YT/lib/python2.6/site-packages/matplotlib/backends/backend_agg.py", line 314, in draw<br> self.figure.draw(self.renderer)<br> File "/home/h-clark/YT/lib/python2.6/site-packages/matplotlib/artist.py", line 46, in draw_wrapper<br>
draw(artist, renderer, *args, **kwargs)<br> File "/home/h-clark/YT/lib/python2.6/site-packages/matplotlib/figure.py", line 773, in draw<br> for a in self.axes: a.draw(renderer)<br> File "/home/h-clark/YT/lib/python2.6/site-packages/matplotlib/artist.py", line 46, in draw_wrapper<br>
draw(artist, renderer, *args, **kwargs)<br> File "/home/h-clark/YT/lib/python2.6/site-packages/matplotlib/axes.py", line 1701, in draw<br> im.draw(renderer)<br> File "/home/h-clark/YT/lib/python2.6/site-packages/matplotlib/artist.py", line 46, in draw_wrapper<br>
draw(artist, renderer, *args, **kwargs)<br> File "/home/h-clark/YT/lib/python2.6/site-packages/matplotlib/image.py", line 237, in draw<br> im = self.make_image(renderer.get_image_magnification())<br> File "/home/h-clark/YT/lib/python2.6/site-packages/matplotlib/image.py", line 182, in make_image<br>
x = self.to_rgba(self._A, self._alpha)<br> File "/home/h-clark/YT/lib/python2.6/site-packages/matplotlib/cm.py", line 167, in to_rgba<br> x = self.cmap(x, alpha=alpha, bytes=bytes)<br>TypeError: 'NoneType' object is not callable<br>
<br clear="all"><br><b>My Script is:</b><br><br>from yt.mods import *<br>import matplotlib.colorbar as cb<br><br>orient = 'horizontal'<br>frame_template = "aaSN/%s/multiplot_%01i_%04i"<br>centre = (0.511901, 0.452331, 0.434998)<br>
ax = 0<br>ss = 1<br>n = 134<br><br>pf = load("SciNet/%s/DD%04i/data%04i" % (sys.argv[-1],n, n))<br>fig, axes, colorbars = get_multi_plot( 3, 1, colorbar=orient, bw = 4)<br>pc = PlotCollection(pf, center=centre)<br>
<br>p = pc.add_slice("Density", ax, figure = fig, axes = axes[0][0], use_colorbar=False)<br>p.set_cmap("bds_highcontrast")<br>p = pc.add_slice("Temperature", ax, figure=fig, axes=axes[0][1], use_colorbar=False)<br>
p.set_cmap("hot")<br>p = pc.add_slice("HII_Fraction", ax, figure=fig, axes=axes[0][2], use_colorbar=False)<br>p.set_cmap("cold")<br><br>p.set_width(ss, 'pc')<br>for p, cax in zip(pc.plots, colorbars):<br>
cbar = cb.Colorbar(cax, p.image, orientation=orient)<br> p.colorbar = cbar<br> p._autoset_label()<br>fig.savefig(frame_template % (sys.argv[-1],ax,n))<br><br><br><br>What am I missing?<br><br>many thanks,<br><br>
Libby<br><font color="#888888"><br>-- <br>Elizabeth Harper-Clark MA MSci<br>PhD Candidate, Canadian Institute for Theoretical Astrophysics, UofT<br>Sciences and Engineering Coordinator, Teaching Assistants' Training Program, UofT<br>
<br><a href="http://www.astro.utoronto.ca/%7Eh-clark" target="_blank">www.astro.utoronto.ca/~h-clark</a><br>
<a href="mailto:h-clark@cita.utoronto.ca" target="_blank">h-clark@cita.utoronto.ca</a><br>Astronomy office phone: +1-416-978-5759<br>
</font></blockquote></div><br><br clear="all"><br>-- <br>Elizabeth Harper-Clark MA MSci<br>PhD Candidate, Canadian Institute for Theoretical Astrophysics, UofT<br>Sciences and Engineering Coordinator, Teaching Assistants' Training Program, UofT<br>
<br><a href="http://www.astro.utoronto.ca/%7Eh-clark" target="_blank">www.astro.utoronto.ca/~h-clark</a><br><a href="mailto:h-clark@cita.utoronto.ca" target="_blank">h-clark@cita.utoronto.ca</a><br>Astronomy office phone: +1-416-978-5759<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Elizabeth Harper-Clark MA MSci<br>PhD Candidate, Canadian Institute for Theoretical Astrophysics, UofT<br>Sciences and Engineering Coordinator, Teaching Assistants' Training Program, UofT<br>
<br><a href="http://www.astro.utoronto.ca/%7Eh-clark" target="_blank">www.astro.utoronto.ca/~h-clark</a><br><a href="mailto:h-clark@cita.utoronto.ca" target="_blank">h-clark@cita.utoronto.ca</a><br>Astronomy office phone: +1-416-978-5759<br>