<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi Christine, <br>
<br>
thanks a lot for your script, it does exactly what I need. <br>
<br>
The last thing I would need is being able to overplot the particles
positions and the velocity field. I tried to add the following after
creating the slice:<br>
<br>
######################################################<br>
<address>p = pc0.plots[-1].modify["velocity"]()</address>
<address>p = pc0.plots[-1].modify["particles"](width=10,
p_size=10.0, col='g')</address>
######################################################<br>
<br>
Although it does not make the code break, the slices remain
unchanged. <br>
<br>
Do you have any idea how to make it work ? If not, I will follow
Britton's advice and dig a little bit more into the more complex
yt-advanced-plotting rep.<br>
<br>
Thanks a lot,<br>
<br>
JC<br>
<br>
<br>
On 22/03/11 06:24, Christine Simpson wrote:
<blockquote
cite="mid:1300735485.3504.108.camel@cressida.astro.columbia.edu"
type="cite">
<pre wrap="">Hi JC,
Matt says the yt developers are working on some more sophisticated
imaging capabilities, but in the meantime, this script does something
close to what you want. It's a bit cumbersome, but it may be helpful to
you. It plots density slices in the x,y,z directions for two different
datasets in side-by-side columns. In general, I think you just need to
set up separate plot collections for each dataset you want to plot. The
tricky part is getting the colorbars to work; the way i've done it
below, you can in can have different colorbars for the two different
columns.:
from yt.mods import * # set up our namespace
import matplotlib.colorbar as cb
fn0 = "DD0007/output_0007"
fn1 = "DD0010/output_0010"
pf0 = load(fn0) # load data
pf1 = load(fn1)
orient = 'horizontal'
fig, axes, colorbars = get_multi_plot( 2, 3, colorbar=orient, bw = 4)
pc0 = PlotCollection(pf0)
pc1 = PlotCollection(pf1)
for ax in range(3):
p = pc0.add_slice("Density", ax, figure = fig, axes = axes[ax][0],
use_colorbar=False)
p.set_cmap("bds_highcontrast")
p.set_zlim(10**(-27),10**(-25))
p = pc1.add_slice("Density", ax, figure=fig, axes=axes[ax][1],
use_colorbar=False)
p.set_cmap("bds_highcontrast")
p.set_zlim(10**(-30),10**(-25))
pc0.set_width(25.0, 'kpc')
pc1.set_width(25.0, 'kpc')
zip_plot = zip(pc0.plots,pc1.plots)
all_plots = [zip_plot[i][j] for i in range(3) for j in range(2)]
for p, cax in zip(all_plots, colorbars):
cbar = cb.Colorbar(cax, p.image, orientation=orient)
p.colorbar = cbar
p._autoset_label()
fig.savefig("Density_before_after_reion")
On Sun, 2011-03-20 at 17:06 -0700, Matthew Turk wrote:
</pre>
<blockquote type="cite">
<pre wrap="">Hi Chris and JC,
You're right; they are not up to date in the build. It was a mistake
in how I built the docs last. I updated all the *API* content, but
not the content of the *recipes*. When I get back to my desk tonight
(I'm out enjoying some Sunday afternoon recreation time right now) I
will rebuild the docs, but I'll build from the development branch so
they will reflect the upcoming 2.1 release in a few minor ways, like
minor shifts in the light cone API. Sam's release manager for 2.1,
and last I heard from him he's planning on late this week or early
next.
Also, JC, Britton, Sam and Jeff have been working on making it easier
to do more complicated things with montages and compositing of
arbitrary images. They might be able to chime in with more; multiplot
is kind of a sledge hammer and they're putting the finishing touches
on their nice set of drill bits and screw drivers.
Best,
Matt
On Sun, Mar 20, 2011 at 4:59 PM, <chris.m.malone@gmail.com> wrote:
</pre>
<blockquote type="cite">
<pre wrap="">I should say, the recipe is up to date (in the hg repository), but the
Cookbook website is not. Same goes for the regular Multiplot recipe.
Matt, could you update these?
Chris
On Mar 20, 2011 7:47pm, chris.m.malone@gmail.com wrote:
</pre>
<blockquote type="cite">
<pre wrap="">Hi Jean-Claude,
Some of the recipes in the Cookbook are a bit outdated - for example,
module names like "raven" have been replaced with more meaningful names.
Actually, now the get_multi_plot() method lives in the plot_collection
module, which is automatically loaded into the namespace when you `from
yt.mods import *`. In other words, change the `raven.get_mult_plot(...)` to
just `get_multi_plot(...)` and it should work ok, although I haven't tested
this.
Chris
On Mar 20, 2011 7:27pm, Jean-Claude Passy jcpassy@gmail.com> wrote:
</pre>
<blockquote type="cite">
<pre wrap="">
Hi all,
I would like to create a multi plot equivalent to the procedure
described in the Cookbook (Multi plot 3x2) BUT with
different datasets (1 slice for 1 dataset).
The first issue I am encountering is that the module raven is not
detected:
#############################################################
Traceback (most recent call last):
File "/rpod2/jcpassy/Enzo/python/MultiPlot.py", line 20,
in
fig, axes, colorbars = raven.get_multi_plot( 2, 4,
colorbar=orient, bw = 4)
NameError: name 'raven' is not defined
#############################################################
Then, is there a specific description somewhere of how to do a multi
plot with different datasets ?
Thanks a lot for your help,
JC
</pre>
</blockquote>
</blockquote>
<pre wrap="">_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
</pre>
</blockquote>
<pre wrap="">_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
</pre>
</blockquote>
<pre wrap="">
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
</pre>
</blockquote>
</body>
</html>