<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
Hi Matt.
<div><br>
</div>
<div>The Script you sent me has been of great help! I have now managed to get yt to plot projections for my data. </div>
<div><br>
</div>
<div>However I am still having small trouble. Firstly script only plots the 1 quarter of each projection(I can change which quarter), [Graph below] and secondly I would like to set the colorbar to be linear for my data. How can I go about that?</div>
<div><br>
</div>
<div>Thanks</div>
<div><br>
</div>
<div>Kearn</div>
<div><br>
</div>
<div><img height="199" width="320" apple-width="yes" apple-height="yes" apple-inline="yes" id="AC21915C-7105-4D0A-B14F-617309B8D833" src="cid:2ABCAB42-4A79-4A32-83DD-88FE8CEBB9E9@surrey.ac.uk"></div>
<div>
<div>
<div>On 11 Nov 2013, at 22:39, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div>Hi Kearn,<br>
<br>
Oops, I guess my reply was less than helpful! Sorry about that --<br>
been one of those days.<br>
<br>
On Mon, Nov 11, 2013 at 3:30 PM, <<a href="mailto:k.grisdale@surrey.ac.uk">k.grisdale@surrey.ac.uk</a>> wrote:<br>
<blockquote type="cite">Hi Matt,<br>
<br>
Thanks for the reply. The code that I used below is my attempt at adapting<br>
the script on the yt pages explaining how to use Multi-Plot Slice and<br>
Projections. As I couldn’t figure out how make multi-plots of yt Projections<br>
I thought I’d try the adapting route. Which means there bits there I don’t<br>
get but have included as I thought they where needed for the multi plot.<br>
<br>
</blockquote>
<br>
I've spent a bit of time looking at this and I've been able to make a<br>
multi plot that I think matches what you're looking for. I've put the<br>
script here, and I've also tried to note where there might be issues<br>
with it that could be important when you adapt it to your data.<br>
<br>
<a href="http://paste.yt-project.org/show/4047/">http://paste.yt-project.org/show/4047/</a><br>
<br>
Let me know if that helps out, or if I can help in any other way.<br>
<br>
-Matt<br>
<br>
<blockquote type="cite">Thanks again<br>
<br>
Kearn<br>
<br>
<br>
<br>
On 11 Nov 2013, at 20:13, Matthew Turk <matthewturk@gmail.com> wrote:<br>
<br>
Hi Kearn,<br>
<br>
<br>
On Mon, Nov 11, 2013 at 10:18 AM, <k.grisdale@surrey.ac.uk> wrote:<br>
<blockquote type="cite"><br>
Hi everyone<br>
<br>
I am trying to make a multi plot using data from Ramses. I have outputs<br>
from Ramses at 6 different times and I want to make projections at each time<br>
and place than in a multi plot to show evolution of the system. However when<br>
the code that I have written gives the results below.<br>
<br>
<br>
<Screen Shot 2013-11-11 at 15.12.13.png><br>
<br>
<br>
Can anyone offer any suggestions as to wear my code is going wrong?<br>
<br>
Thanks<br>
<br>
Kearn<br>
<br>
The code I’m using ( the for dax lines have been commented out because<br>
every time they are included I get the error at the bottom)<br>
<br>
from yt.mods import * # set up our namespace<br>
import matplotlib.colorbar as cb<br>
from matplotlib.colors import LogNorm<br>
<br>
p1 = load("output_00001/info_00001.txt")<br>
p2 = load("output_00002/info_00002.txt")<br>
p3 = load("output_00003/info_00003.txt")<br>
p4 = load("output_00004/info_00004.txt")<br>
p5 = load("output_00005/info_00005.txt")<br>
p6 = load("output_00006/info_00006.txt")<br>
<br>
cen = [0,0,0]<br>
orient = 'vertical'<br>
<br>
fig, axes, colorbars = get_multi_plot(3, 2, colorbar=orient, bw = 4)<br>
proj1 = p1.h.proj(2, "Density", center=cen)<br>
proj2 = p2.h.proj(2, "Density", center=cen)<br>
proj3 = p3.h.proj(2, "Density", center=cen)<br>
proj4 = p4.h.proj(2, "Density", center=cen)<br>
proj5 = p5.h.proj(2, "Density", center=cen)<br>
proj6 = p6.h.proj(2, "Density", center=cen)<br>
<br>
proj1_frb = proj1.to_frb((0.7,"cm"),512)<br>
proj2_frb = proj2.to_frb((0.7,"cm"),512)<br>
proj3_frb = proj3.to_frb((0.7,"cm"),512)<br>
proj4_frb = proj4.to_frb((0.7,"cm"),512)<br>
proj5_frb = proj5.to_frb((0.7,"cm"),512)<br>
proj6_frb = proj6.to_frb((0.7,"cm"),512)<br>
<br>
dens_axes = [axes[0][0], axes[1][0]]<br>
<br>
#for dax in zip(dens_axes) :<br>
</blockquote>
<br>
<br>
I'm not sure what you're doing here -- this will yield:<br>
<br>
[(axes[0][0],), (axes[1][0],)]<br>
<br>
So you'll be getting back a tuple of length one inside the loop. You should<br>
be able to do dax[0]? But I'm not sure you need to zip in the first place,<br>
either.<br>
<br>
-Matt<br>
<br>
<blockquote type="cite"><br>
<br>
# dax.xaxis.set_visible(False)<br>
# dax.yaxis.set_visible(False)<br>
<br>
plots = [dens_axes[0].imshow(proj1_frb["Density"],origin ='lower'),<br>
dens_axes[0].imshow(proj2_frb["Density"],origin ='lower'),<br>
dens_axes[0].imshow(proj3_frb["Density"],origin ='lower'),<br>
dens_axes[0].imshow(proj4_frb["Density"],origin ='lower'),<br>
dens_axes[0].imshow(proj5_frb["Density"],origin ='lower'),<br>
dens_axes[0].imshow(proj6_frb["Density"],origin ='lower')]<br>
<br>
plots[0].set_clim((0.0,0.34))<br>
plots[0].set_cmap("YlOrRd")<br>
plots[1].set_clim((0.0,0.345))<br>
plots[1].set_cmap("YlOrRd")<br>
plots[2].set_clim((0.0,0.34))<br>
plots[2].set_cmap("YlOrRd")<br>
plots[3].set_clim((0.0,0.34))<br>
plots[3].set_cmap("YlOrRd")<br>
plots[4].set_clim((0.0,0.34))<br>
plots[4].set_cmap("YlOrRd")<br>
plots[5].set_clim((0.0,0.34))<br>
plots[5].set_cmap("YlOrRd")<br>
<br>
fig.savefig("test”)<br>
<br>
<br>
error:<br>
<br>
Traceback (most recent call last):<br>
File "plotterm2.py", line 34, in <module><br>
dax.xaxis.set_visible(False)<br>
AttributeError: 'tuple' object has no attribute 'xaxis'<br>
<br>
_______________________________________________<br>
yt-users mailing list<br>
yt-users@lists.spacepope.org<br>
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org<br>
<br>
</blockquote>
<br>
_______________________________________________<br>
yt-users mailing list<br>
yt-users@lists.spacepope.org<br>
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org<br>
<br>
<br>
_______________________________________________<br>
yt-users mailing list<br>
yt-users@lists.spacepope.org<br>
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org<br>
<br>
</blockquote>
_______________________________________________<br>
yt-users mailing list<br>
yt-users@lists.spacepope.org<br>
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org<br>
<br>
</div>
</div>
<br>
</div>
</body>
</html>