<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>I have just applied the pull request using the command:</div>
<div><br>
</div>
<div>cd $YT_DEST/src/yt-hg</div>
<div>hg pull && hg update yt</div>
<div><br>
</div>
<div>But now when i tried to load RAMSES data into the YT I get the following error:</div>
<div><br>
</div>
<div>
<div>yt : [ERROR ] 2013-11-13 10:22:49,486 Couldn't figure out output type for output_00001/info_00001.txt</div>
<div>Traceback (most recent call last):</div>
<div> File "<stdin>", line 1, in <module></div>
<div> File "/Applications/Code/yt/yt-x86_64/src/yt-hg/yt/convenience.py", line 90, in load</div>
<div> raise YTOutputNotIdentified(args, kwargs)</div>
<div>yt.utilities.exceptions.YTOutputNotIdentified: Supplied ('output_00001/info_00001.txt',) {}, but could not load!</div>
</div>
<div><br>
</div>
<div>Any idea what has happened?</div>
<div><br>
</div>
<div>Thanks</div>
<div><br>
</div>
<div>Kearn</div>
<div><br>
</div>
<div><br>
</div>
<div>
<div>Traceback (most recent call last):</div>
<div> File "plotterm3.py", line 7, in <module></div>
<div> p1 = load("output_00001/info_00001.txt")</div>
<div> File "/Applications/Code/yt/yt-x86_64/src/yt-hg/yt/convenience.py", line 90, in load</div>
<div> raise YTOutputNotIdentified(args, kwargs)</div>
<div>yt.utilities.exceptions.YTOutputNotIdentified: Supplied ('output_00001/info_00001.txt',) {}, but could not load!</div>
</div>
<div><br>
</div>
<div><br>
<div>
<div>On 12 Nov 2013, at 14:57, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div>
<div dir="ltr">Hi Kearn,
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Tue, Nov 12, 2013 at 9:47 AM, <span dir="ltr"><<a href="mailto:k.grisdale@surrey.ac.uk" target="_blank">k.grisdale@surrey.ac.uk</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 style="word-wrap:break-word">Dear Matt,
<div><br>
</div>
<div>Thanks again. The normalize command was just what I needed. </div>
<div><br>
</div>
<div>Does the projection plot used in the multi-plot work differently to the normal ProjectionPlot?</div>
</div>
</blockquote>
<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 style="word-wrap:break-word">
<div><br>
</div>
<div>I ask because with the code below I get the image below, while using the same centre and width with the multi plot I get the image I sent in my last email. I have tried playing around with different centres, and projection radii but this only moves the
image around or makes the square that it does plot different sizes(this is still the case even after changing to a linear scale). Both plots are using the same data set so i am explicating the same plot.</div>
</div>
</blockquote>
<div><br>
</div>
<div>They sort of work the same. What you're doing with the FRBs is much lower-level, though, even though it's what the ProjectionPlot uses. I think the issue you're running into is that the FixedResolutionBuffer is *not* periodic by default, and to_frb()
doesn't turn it on. I've issued a pull request that adds this functionality to the to_frb() command, as I think it is a problem not to have it!</div>
<div><br>
</div>
<div><a href="https://bitbucket.org/yt_analysis/yt/pull-request/639/adding-periodic-default-false-to-to_frb">https://bitbucket.org/yt_analysis/yt/pull-request/639/adding-periodic-default-false-to-to_frb</a></div>
<div><br>
</div>
<div>-Matt</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 style="word-wrap:break-word">
<div><br>
</div>
<div>Thanks Kearn</div>
<div><br>
</div>
<div><br>
</div>
<div>
<div>from yt.mods import *</div>
<div class="im">
<div>p6 = load("output_00006/info_00006.txt")</div>
<div>cen = [0,0,0]</div>
</div>
<div>pr6 = ProjectionPlot(p6,'z','Density',center=cen)</div>
<div>pr6.set_cmap('Density','YlOrRd')</div>
<div>pr6.set_zlim('Density', 0.0,0.35)</div>
<div>pr6.save("p6")</div>
</div>
<div><br>
</div>
<div><br>
</div>
<div><span><p6_Projection_z_Density.png></span></div>
<div><br>
</div>
<div>
<div>
<div>
<div class="im">
<div>On 12 Nov 2013, at 11:33, Matthew Turk <<a href="mailto:matthewturk@gmail.com" target="_blank">matthewturk@gmail.com</a>> wrote:</div>
<br>
</div>
<div>
<div dir="ltr">Hi Kearn,
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">
<div class="im">On Tue, Nov 12, 2013 at 3:46 AM, <span dir="ltr"><<a href="mailto:k.grisdale@surrey.ac.uk" target="_blank">k.grisdale@surrey.ac.uk</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 style="word-wrap:break-word">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>
</blockquote>
<div><br>
</div>
<div>If you change LogNorm to Normalize (and change the import) it should change to linear. You can probably then just remove the LogLocator ticker. To change the FOV, the 0.25 "unitary" is what I set, and you should be able to change it to whatever you like.
Just be sure to also change the center value to be the center of the domain! :)</div>
<div><br>
</div>
<div>-Matt</div>
<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 style="word-wrap:break-word">
<div><br>
</div>
<div>Thanks</div>
<div><br>
</div>
<div>Kearn</div>
<div><br>
</div>
<div><span><Screen Shot 2013-11-12 at 08.46.03.png></span></div>
<div>
<div class="h5">
<div>
<div>
<div>
<div>
<div>On 11 Nov 2013, at 22:39, Matthew Turk <<a href="mailto:matthewturk@gmail.com" target="_blank">matthewturk@gmail.com</a>> wrote:</div>
<br>
<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" target="_blank">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/" target="_blank">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 <<a href="mailto:matthewturk@gmail.com" target="_blank">matthewturk@gmail.com</a>> wrote:<br>
<br>
Hi Kearn,<br>
<br>
<br>
On Mon, Nov 11, 2013 at 10:18 AM, <<a href="mailto:k.grisdale@surrey.ac.uk" target="_blank">k.grisdale@surrey.ac.uk</a>> 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>
<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>
<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>
<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>
_______________________________________________<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>
</div>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div class="h5"><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>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
<div>
<div class="h5">_______________________________________________<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>
</div>
</div>
</div>
</div>
<br>
</div>
</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>
</div>
</div>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org<br>
</div>
</div>
<br>
</div>
</body>
</html>