<div dir="ltr">Dear all,<br><div class="gmail_quote"><div dir="ltr"><div class="gmail_quote"><div dir="ltr"><div><br></div><div>I am trying to create a muti plot of density's at different red-shift with the following code. </div>
<div><br></div><div><div>
<font face="courier new, monospace" size="4">from yt.mods import * </font></div>
<div><font face="courier new, monospace" size="4">import matplotlib.colorbar as cb</font></div><div><font face="courier new, monospace" size="4"><br></font></div><div><font face="courier new, monospace" size="4">fig, axes, colorbars = get_multi_plot( 3, 3, colorbar='vertical', bw = 3)</font></div>
<div><font face="courier new, monospace" size="4">file1="/app/run/csarkar/surajit/simulation_shock/simulation5/RD00"</font></div><div><font face="courier new, monospace" size="4">file2="/RedshiftOutput00"</font></div>
<div><font face="courier new, monospace" size="4"><br></font></div>
<div><font face="courier new, monospace" size="4">L = [11,15,18,20,22,25,30,35,40]</font></div><div><font face="courier new, monospace" size="4"><br></font></div><div><font face="courier new, monospace" size="4">#for i in range (0, 8):</font></div>
<div><font face="courier new, monospace" size="4">for kk in L: </font></div><div><font face="courier new, monospace" size="4"> fn_i=str(kk)</font></div><div><font face="courier new, monospace" size="4"> fn=file1+fn_i+file2+fn_i </font></div>
<div><font face="courier new, monospace" size="4"> print fn</font></div><div><font face="courier new, monospace" size="4"> pf = load(fn) </font></div><div><font face="courier new, monospace" size="4"> pc = PlotCollection(pf, center=[0.5, 0.5, 0.5])</font></div>
<div><font face="courier new, monospace" size="4"> for i in range (0,3):</font></div><div><font face="courier new, monospace" size="4"> for j in range (0, 3):</font></div><div><font face="courier new, monospace" size="4"> p = pc.add_slice("Density", 2, axes = axes[i][j],figure = fig, use_colorbar=False)</font></div>
<div><font face="courier new, monospace" size="4"> p.set_log_field(True)</font></div><div><font face="courier new, monospace" size="4"> p.set_zlim(1e-31,5e-26)</font></div><div><font face="courier new, monospace" size="4"> p.set_cmap("jet") </font></div>
<div><font face="courier new, monospace" size="4"> p.modify["contour"]("Density", ncont=10, plot_args = {'colors':'w'},clim=(2e-28,2e-26)) </font></div><div><font face="courier new, monospace" size="4"> p.modify["contour"]("Temperature", ncont=10, plot_args = {'colors':'k'},clim=(3e7,1e8))</font></div>
<div><font face="courier new, monospace" size="4"> #pc.set_width(3.0, 'mpc') </font></div><div><font face="courier new, monospace" size="4"> #break</font></div><div><font face="courier new, monospace" size="4"> # i = i+1</font></div>
<div><font face="courier new, monospace" size="4"><br></font></div><div><font face="courier new, monospace" size="4"> </font></div><div><font face="courier new, monospace" size="4">for p, cax in zip(pc.plots, colorbars):</font></div>
<div><font face="courier new, monospace" size="4"> cbar = cb.Colorbar(cax, p.image, orientation='vertical')</font></div><div><font face="courier new, monospace" size="4"> p.colorbar = cbar</font></div><div>
<font face="courier new, monospace" size="4"> p._autoset_label()</font></div>
<div><font face="courier new, monospace" size="4">fig.savefig("multi_plot_3_3" % pf)</font></div></div><div><font face="courier new, monospace" size="4"><br></font></div><div><font face="courier new, monospace" size="4"><br>
</font></div><div><font face="arial, helvetica, sans-serif">But I am getting a wrong result. It is printing same figure in all panels.</font></div><div><br></div><div><font face="arial, helvetica, sans-serif">I think my logic(algorithm) is wrong. We have to reverse first two for loop. But that also making problem. Please help me. </font><span style="font-family:arial,helvetica,sans-serif"> </span></div>
<div><span style="font-family:arial,helvetica,sans-serif"><br></span></div><div><font face="arial, helvetica, sans-serif">And also if I include the following line in the code. (in my code I commented it out)</font></div>
<div><font face="arial, helvetica, sans-serif"><br>
</font></div><div><span style="font-family:'courier new',monospace;font-size:large">pc.set_width(3.0, 'mpc')</span><font face="arial, helvetica, sans-serif"><br></font></div><div><br></div><div><font face="arial, helvetica, sans-serif">it is showing following error. </font></div>
<div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif"><div><i>Traceback (most recent call last):</i></div><div><i> File "multi_plot_loop.py", line 24, in <module></i></div>
<div><i> pc.set_width(3.0, 'mpc') </i></div><div><i> File "/data1/pdf/csurajit/yt-x86_64/src/yt-hg/yt/visualization/plot_collection.py", line 292, in set_width</i></div><div><i> plot.set_width(width, unit)</i></div>
<div><i> File "/data1/pdf/csurajit/yt-x86_64/src/yt-hg/yt/visualization/plot_types.py", line 458, in set_width</i></div><div><i> self._refresh_display_width()</i></div><div><i> File "/data1/pdf/csurajit/yt-x86_64/src/yt-hg/yt/visualization/plot_types.py", line 476, in _refresh_display_width</i></div>
<div><i> self._redraw_image()</i></div><div><i> File "/data1/pdf/csurajit/yt-x86_64/src/yt-hg/yt/visualization/plot_types.py", line 423, in _redraw_image</i></div><div><i> self._run_callbacks()</i></div>
<div>
<i> File "/data1/pdf/csurajit/yt-x86_64/src/yt-hg/yt/visualization/plot_types.py", line 263, in _run_callbacks</i></div><div><i> cb(self)</i></div><div><i> File "/data1/pdf/csurajit/yt-x86_64/src/yt-hg/yt/visualization/plot_modifications.py", line 312, in __call__</i></div>
<div><i> self.ncont = np.linspace(self.clim[0], self.clim[1], self.ncont)</i></div><div><i> File "/data1/pdf/csurajit/yt-x86_64/lib/python2.7/site-packages/numpy/core/function_base.py", line 74, in linspace</i></div>
<div><i> num = int(num)</i></div><div><i>TypeError: only length-1 arrays can be converted to Python scalars</i></div><div><br></div></font></div><div><br></div><font color="#000000" face="Arial, Liberation Sans, DejaVu Sans, sans-serif"><span style="font-size:14px;line-height:18px">Please suggest me the most efficient way to achieve my goal (</span></font>create a muti plot of density's at different red-shift<font color="#000000" face="Arial, Liberation Sans, DejaVu Sans, sans-serif"><span style="font-size:14px;line-height:18px">).</span></font><span class="HOEnZb"><font color="#888888"><span><font color="#888888"><br clear="all">
<div><br></div></font></span></font></span></div></div></div></div><br><br clear="all"><div><br></div>-- <br><div dir="ltr">Reju Sam John</div>
</div>