Just in case you haven't updated in a while, I would do so since there was a fix included in yt 3.2.1 that might solve your problem.<div><br></div><div><a href="https://bitbucket.org/yt_analysis/yt/pull-requests/1698/ensure-the-frb-is-valid-before-executing/diff">https://bitbucket.org/yt_analysis/yt/pull-requests/1698/ensure-the-frb-is-valid-before-executing/diff</a><br><div><br></div><div>If you're already on yt 3.2.1 <span></span>I'll take a closer look at this soon.<br><div><br></div><div>Nathan<br><br>On Saturday, September 12, 2015, Guido granda muñoz <<a href="mailto:guidogranda@gmail.com">guidogranda@gmail.com</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Hello there,<br><br></div>I'm trying to do a multipanel plot of projected density with titles in each axis and annotations but I couldn't figure out how to do it, could you please help me.<br><br></div><div>I'm using AxesGrid as in your example(<a href="http://yt-project.org/doc/cookbook/complex_plots.html#multipanel-with-axes-labels" target="_blank">http://yt-project.org/doc/cookbook/complex_plots.html#multipanel-with-axes-labels</a>)  however<br></div><div>something is wrong. I'm sending my output image as a attached file.<br><br></div>My code is the following:<br><br><br><br><br>import yt<br>import matplotlib.pyplot as plt<br>import sys<br>import numpy as np<br>from mpl_toolkits.axes_grid1 import AxesGrid<br>import matplotlib.colorbar as cb<br>from matplotlib.colors import LogNorm<br><br>data=[]<br>data.append(sys.argv[1])<br>data.append(sys.argv[2])<br>data.append(sys.argv[3])<br>data.append(sys.argv[4])<br><br>pf1=yt.load(data[0])<br>pf2=yt.load(data[1])<br>pf3=yt.load(data[2])<br>pf4=yt.load(data[3])<br><br>pfs=[]<br>pfs.append(pf1)<br>pfs.append(pf2)<br>pfs.append(pf3)<br>pfs.append(pf4)<br><br><br>fig=plt.figure()<br>grid = AxesGrid(fig, (0.075,0.075,0.85,0.85),nrows_ncols = (2, 2),axes_pad = 1.0,label_mode = "all",share_all = True,cbar_location="right",cbar_mode="single",cbar_size="5%",cbar_pad="5%",cbar_set_cax=True)<br><br><br>axe=['y','y','y','y']<br>field=['density','density','density','density']<br>colors=[]<br>title=['100 Mpc/h box size', '75 Mpc/h box size','100 Mpc/h box size','50 Mpc/h box size']<br><br>for i in range(len(data)):<br>    p   =yt.ProjectionPlot(pfs[i],axe[i],field[i],center='c',width=(25,'Mpc/h'))<br>        p.set_zlim('density',4.4601764879407653e-5,5.0e-2)<br>        p.annotate_timestamp(corner='upper_left', redshift=True,draw_inset_box=True)<br>        p.annotate_title(title[i])<br>        plot=p.plots[field]<br>        plot.figure=fig<br>        plot.axes=grid[i].axes<br>        plot.cax=grid.cbar_axes[i]<br>        p._setup_plots()<br><br><br>fig.savefig('plot1_f.png',pdi=300)<br clear="all"><div><div><div><br><br><br><br></div><div>Kind regards,<br></div><div>-- <br><div>Guido</div>
</div></div></div></div>
</blockquote></div></div></div>