[yt-users] Making a Multiplot using Ramses data

k.grisdale at surrey.ac.uk k.grisdale at surrey.ac.uk
Mon Nov 11 07:18:20 PST 2013


Hi everyone

I am trying to make a multi plot using data from Ramses. I have outputs from Ramses at 6 different times and I want to make projections at each time and place than in a multi plot to show evolution of the system. However when the code that I have written gives the results below.


[cid:F88B06FF-D341-4589-9FD2-C8BEF5BE7CE8 at surrey.ac.uk]


Can anyone offer any suggestions as to wear my code is going wrong?

Thanks

Kearn

The code I’m using ( the for dax lines have been commented out because every time they are included I get the error at the bottom)

from yt.mods import * # set up our namespace
import matplotlib.colorbar as cb
from matplotlib.colors import LogNorm

p1 = load("output_00001/info_00001.txt")
p2 = load("output_00002/info_00002.txt")
p3 = load("output_00003/info_00003.txt")
p4 = load("output_00004/info_00004.txt")
p5 = load("output_00005/info_00005.txt")
p6 = load("output_00006/info_00006.txt")

cen = [0,0,0]
orient = 'vertical'

fig, axes, colorbars = get_multi_plot(3, 2, colorbar=orient, bw = 4)
proj1 = p1.h.proj(2, "Density", center=cen)
proj2 = p2.h.proj(2, "Density", center=cen)
proj3 = p3.h.proj(2, "Density", center=cen)
proj4 = p4.h.proj(2, "Density", center=cen)
proj5 = p5.h.proj(2, "Density", center=cen)
proj6 = p6.h.proj(2, "Density", center=cen)

proj1_frb = proj1.to_frb((0.7,"cm"),512)
proj2_frb = proj2.to_frb((0.7,"cm"),512)
proj3_frb = proj3.to_frb((0.7,"cm"),512)
proj4_frb = proj4.to_frb((0.7,"cm"),512)
proj5_frb = proj5.to_frb((0.7,"cm"),512)
proj6_frb = proj6.to_frb((0.7,"cm"),512)

dens_axes = [axes[0][0], axes[1][0]]

#for dax in zip(dens_axes) :

#    dax.xaxis.set_visible(False)
#    dax.yaxis.set_visible(False)

plots = [dens_axes[0].imshow(proj1_frb["Density"],origin ='lower'),
         dens_axes[0].imshow(proj2_frb["Density"],origin ='lower'),
         dens_axes[0].imshow(proj3_frb["Density"],origin ='lower'),
         dens_axes[0].imshow(proj4_frb["Density"],origin ='lower'),
         dens_axes[0].imshow(proj5_frb["Density"],origin ='lower'),
         dens_axes[0].imshow(proj6_frb["Density"],origin ='lower')]

plots[0].set_clim((0.0,0.34))
plots[0].set_cmap("YlOrRd")
plots[1].set_clim((0.0,0.345))
plots[1].set_cmap("YlOrRd")
plots[2].set_clim((0.0,0.34))
plots[2].set_cmap("YlOrRd")
plots[3].set_clim((0.0,0.34))
plots[3].set_cmap("YlOrRd")
plots[4].set_clim((0.0,0.34))
plots[4].set_cmap("YlOrRd")
plots[5].set_clim((0.0,0.34))
plots[5].set_cmap("YlOrRd")

fig.savefig("test”)


error:

Traceback (most recent call last):
  File "plotterm2.py", line 34, in <module>
    dax.xaxis.set_visible(False)
AttributeError: 'tuple' object has no attribute 'xaxis'
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20131111/c3aa52e1/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2013-11-11 at 15.12.13.png
Type: image/png
Size: 20138 bytes
Desc: Screen Shot 2013-11-11 at 15.12.13.png
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20131111/c3aa52e1/attachment.png>


More information about the yt-users mailing list