<div dir="ltr"><div>hi, helpful folks<br><br><br></div><div>I am once again very confused, and would appreciate any help. I am still trying to make a plot which contains frb's created from datafiles at different time steps. This is what I have that I think should work, and doesn't:<br>
<br>from yt.mods import *<br>%matplotlib inline<br>import matplotlib.pyplot as plt<br>import matplotlib.colorbar as cb<br>from matplotlib.colors import LogNorm<br>import numpy as np<br><br>ts = TimeSeriesData.from_filenames("Strat_Box_hdf5_plt_cnt_01*")<br>
<br>print len(ts) #just checking<br>num = len(ts)<br><br>for index, val in enumerate(ts):<br> print "%s = %s" % (index, val) #again, just checking that this does what I think<br><br>orient = 'vertical'<br>
fig, axes, colorbars = get_multi_plot( num, 1, colorbar=orient, bw = 6)<br><br>plots = []<br><br>for i, fn in enumerate(ts):<br> pf = load(fn) # load data<br> sl = pf.h.slice(0, pf.domain_center[0],fields=["Density","Temperature"] )<br>
frb = sl.to_frb( (1.0, 'kpc'), (4096,1024),height=(4.0, 'kpc')) <br> dens_axes = [axes[0][i]]<br> for ax in dens_axes:<br> ax.xaxis.set_visible(False)<br> ax.yaxis.set_visible(False)<br>
plots = [dens_axes[0].imshow(frb["Density"], origin='lower', norm=LogNorm())]<br><br></div><div>These are the errors that I get;<br><pre><span class="">TypeError</span> Traceback (most recent call last)
<span class=""><ipython-input-6-cd09eb0aaf88></span> in <span class=""><module><span class="">()</span>
<span class=""> 1</span> <span class="">for</span> </span>i<span class="">,</span> fn <span class="">in</span> enumerate<span class="">(</span>ts<span class="">)</span><span class="">:</span><span class=""></span>
<span class="">----> 2<span class=""> </span>pf</span> <span class="">=</span> load<span class="">(</span>fn<span class="">)</span> <span class=""># load data</span><span class=""></span>
<span class=""> 3</span> sl <span class="">=</span> pf<span class="">.</span>h<span class="">.</span>slice<span class="">(</span><span class="">0</span><span class="">,</span> pf<span class="">.</span>domain_center<span class="">[</span><span class="">0</span><span class="">]</span><span class="">,</span>fields<span class="">=</span><span class="">[</span><span class="">"Density"</span><span class="">,</span><span class="">"Temperature"</span><span class="">]</span> <span class="">)</span><span class=""></span>
<span class=""> 4</span> frb <span class="">=</span> sl<span class="">.</span>to_frb<span class="">(</span> <span class="">(</span><span class="">1.0</span><span class="">,</span> <span class="">'kpc'</span><span class="">)</span><span class="">,</span> <span class="">(</span><span class="">4096</span><span class="">,</span><span class="">1024</span><span class="">)</span><span class="">,</span>height<span class="">=</span><span class="">(</span><span class="">4.0</span><span class="">,</span> <span class="">'kpc'</span><span class="">)</span><span class="">)</span><span class=""></span>
<span class=""> 5</span> dens_axes <span class="">=</span> <span class="">[</span>axes<span class="">[</span><span class="">0</span><span class="">]</span><span class="">[</span>i<span class="">]</span><span class="">]</span><span class=""></span>
<span class="">/Users/kde/yt-x86_64/src/yt-hg/yt/convenience.pyc</span> in <span class="">load<span class="">(*args, **kwargs)</span>
<span class=""> 65</span> <span class="">try</span><span class="">:</span><span class=""></span></span>
<span class=""> 66</span> <span class="">from</span> yt<span class="">.</span>data_objects<span class="">.</span>time_series <span class="">import</span> TimeSeriesData<span class=""></span>
<span class="">---> 67<span class=""> </span>ts</span> <span class="">=</span> TimeSeriesData<span class="">.</span>from_filenames<span class="">(</span><span class="">*</span>args<span class="">,</span> <span class="">**</span>kwargs<span class="">)</span><span class=""></span>
<span class=""> 68</span> <span class="">return</span> ts<span class=""></span>
<span class=""> 69</span> <span class="">except</span> YTOutputNotIdentified<span class="">:</span><span class=""></span>
<span class="">/Users/kde/yt-x86_64/src/yt-hg/yt/data_objects/time_series.pyc</span> in <span class="">from_filenames<span class="">(cls, filenames, parallel, **kwargs)</span>
<span class=""> 263</span> </span>filenames <span class="">=</span> glob<span class="">.</span>glob<span class="">(</span>filenames<span class="">)</span><span class=""></span>
<span class=""> 264</span> filenames<span class="">.</span>sort<span class="">(</span><span class="">)</span><span class=""></span>
<span class="">--> 265<span class=""> </span>obj</span> <span class="">=</span> cls<span class="">(</span>filenames<span class="">[</span><span class="">:</span><span class="">]</span><span class="">,</span> parallel <span class="">=</span> parallel<span class="">,</span> <span class="">**</span>kwargs<span class="">)</span><span class=""></span>
<span class=""> 266</span> <span class="">return</span> obj<span class=""></span>
<span class=""> 267</span> <span class=""></span>
<span class="">/Users/kde/yt-x86_64/src/yt-hg/yt/data_objects/static_output.pyc</span> in <span class="">__getitem__<span class="">(self, key)</span>
<span class=""> 146</span> for d in [self.units, self.time_units, self.parameters, \
<span class=""> 147</span> self.conversion_factors]:
<span class="">--> 148<span class=""> <span class="">if</span> </span>key</span> <span class="">in</span> </span>d<span class="">:</span> <span class="">return</span> d<span class="">[</span>key<span class="">]</span><span class=""></span>
<span class=""> 149</span> <span class="">raise</span> KeyError<span class="">(</span>key<span class="">)</span><span class=""></span>
<span class=""> 150</span> <span class=""></span>
<span class="">TypeError</span>: unhashable type
</pre><br></div><div><div>thanks<br>kathy<br><br></div><div>-- <br>Kathy DeGioia Eastwood, Ph.D.<br>Professor of Physics and Astronomy<br>Northern Arizona University<br>Flagstaff, AZ 86011-6010<br>Ph: 928-523-7159 FX: 928-523-1371<br>
<a href="mailto:Kathy.Eastwood@nau.edu" target="_blank">Kathy.Eastwood@nau.edu</a><br>deliveries: 602 S. Humphreys St., Bldg 19 Rm 209
</div></div></div>