<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body dir="auto"><div><br></div><div><br></div><div><meta http-equiv="Content-Type" content="text/html charset=utf-8">Hello all,<div class=""><br class=""></div><div class="">I’m having trouble getting yt to plot 2 gas-phase diagrams side-by-side in a grid. I can do this fine for a projection plot, phase plots seems to behave differently. Here’s the code — any idea how to get these two plots side by side sharing the same scale/color-bar?</div><div class=""><br class=""></div><blockquote style="margin: 0 0 0 40px; border: none; padding: 0px;" class=""><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""># fractional=True keyword -- divide the profile data by the sum. </span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class="">fns = [ "/Users/earnric/RAMSES-Data/coolingFixedTau/output_00010/info_00010.txt",</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> "/Users/earnric/RAMSES-Data/noCool/output_00010/info_00010.txt"]</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""><br class=""></span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class="">fig = plt.figure()</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class="">grid = AxesGrid(fig, (0.075,0.075,0.85,0.85),</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> nrows_ncols = (1, 2),</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> axes_pad = 0.05,</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> label_mode = "L",</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> share_all = True,</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> cbar_location="right",</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> cbar_mode="single",</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> cbar_size="3%",</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> cbar_pad="0%")</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class="">fig.suptitle("Gas Phase")</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""><br class=""></span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class="">for i, fn in enumerate(fns):</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> ds = yt.load(fn)</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> ad = ds.all_data()</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> p = yt.PhasePlot(ad, "density", "temperature", "cell_mass", x_bins=256, y_bins=256, weight_field=None, fractional=True )</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> plot = p.plots['cell_mass']</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> if i % 2 == 0:</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> ptype = "Gas phase ${H_2}$ cooling "</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> else:</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> ptype = "Gas phase Standard cooling "</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> zStr = ptype + "z = %.1f" % ds.current_redshift </span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> plot.figure = fig</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> plot.axes = grid[i].axes</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> plot.cax = grid.cbar_axes[i]</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> #plot.annotate_title(zStr) # YOU HAVE TO DO THIS AFTER THE AXES calls</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""><br class=""></span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> p['cell_mass'].axes.set_ylim (1e-1, 1e6)</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> p['cell_mass'].axes.set_xlim (1e-29, 1e-24)</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class=""> p['cell_mass'].axes.set_title(zStr, fontsize=20)</span></font></div></div></div><div class=""><div class=""><div class=""><font face="Courier New" class=""><span style="font-size: 12px;" class="">p.show()</span></font></div></div></div></blockquote><div class=""><div class=""><br class=""></div><div class=""><br class=""></div><div apple-content-edited="true" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="font-size: 15px;" class=""><i class="">Cheers,</i></div><div style="font-size: 15px;" class=""><i class=""><br class=""></i></div><div style="font-size: 15px;" class=""><i class="">Rick Sarmento</i></div><div class=""><br></div><div class=""><br class=""></div></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"><br class="Apple-interchange-newline">
</div>
<br class=""></div></div></body></html>