<div dir="ltr">Hi Kathy,<div><br></div><div>You don't need the load call. The ts has already loaded it for you. fn is actually a pf!</div><div><br></div><div>Just do</div><div><br></div><div>for i, pf in enumerate(ts):</div>
<div> sl = pf.h.slice(...)</div><div><br></div><div>and it'll work.</div><div><br></div><div>See here:</div><div><br></div><div><a href="http://yt-project.org/docs/dev/bootcamp/data_objects_and_time_series.html">http://yt-project.org/docs/dev/bootcamp/data_objects_and_time_series.html</a><br>
</div><div><br></div><div>for more examples.</div><div><br></div><div>Let me know if this doesn't help.</div><div><br></div><div>j</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Jan 27, 2014 at 7:20 PM, Kathy Eastwood <span dir="ltr"><<a href="mailto:kathy.eastwood@nau.edu" target="_blank">kathy.eastwood@nau.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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>TypeError</span> Traceback (most recent call last)
<span><ipython-input-6-cd09eb0aaf88></span> in <span><module><span>()</span>
<span> 1</span> <span>for</span> </span>i<span>,</span> fn <span>in</span> enumerate<span>(</span>ts<span>)</span><span>:</span><span></span>
<span>----> 2<span> </span>pf</span> <span>=</span> load<span>(</span>fn<span>)</span> <span># load data</span><span></span>
<span> 3</span> sl <span>=</span> pf<span>.</span>h<span>.</span>slice<span>(</span><span>0</span><span>,</span> pf<span>.</span>domain_center<span>[</span><span>0</span><span>]</span><span>,</span>fields<span>=</span><span>[</span><span>"Density"</span><span>,</span><span>"Temperature"</span><span>]</span> <span>)</span><span></span>
<span> 4</span> frb <span>=</span> sl<span>.</span>to_frb<span>(</span> <span>(</span><span>1.0</span><span>,</span> <span>'kpc'</span><span>)</span><span>,</span> <span>(</span><span>4096</span><span>,</span><span>1024</span><span>)</span><span>,</span>height<span>=</span><span>(</span><span>4.0</span><span>,</span> <span>'kpc'</span><span>)</span><span>)</span><span></span>
<span> 5</span> dens_axes <span>=</span> <span>[</span>axes<span>[</span><span>0</span><span>]</span><span>[</span>i<span>]</span><span>]</span><span></span>
<span>/Users/kde/yt-x86_64/src/yt-hg/yt/convenience.pyc</span> in <span>load<span>(*args, **kwargs)</span>
<span> 65</span> <span>try</span><span>:</span><span></span></span>
<span> 66</span> <span>from</span> yt<span>.</span>data_objects<span>.</span>time_series <span>import</span> TimeSeriesData<span></span>
<span>---> 67<span> </span>ts</span> <span>=</span> TimeSeriesData<span>.</span>from_filenames<span>(</span><span>*</span>args<span>,</span> <span>**</span>kwargs<span>)</span><span></span>
<span> 68</span> <span>return</span> ts<span></span>
<span> 69</span> <span>except</span> YTOutputNotIdentified<span>:</span><span></span>
<span>/Users/kde/yt-x86_64/src/yt-hg/yt/data_objects/time_series.pyc</span> in <span>from_filenames<span>(cls, filenames, parallel, **kwargs)</span>
<span> 263</span> </span>filenames <span>=</span> glob<span>.</span>glob<span>(</span>filenames<span>)</span><span></span>
<span> 264</span> filenames<span>.</span>sort<span>(</span><span>)</span><span></span>
<span>--> 265<span> </span>obj</span> <span>=</span> cls<span>(</span>filenames<span>[</span><span>:</span><span>]</span><span>,</span> parallel <span>=</span> parallel<span>,</span> <span>**</span>kwargs<span>)</span><span></span>
<span> 266</span> <span>return</span> obj<span></span>
<span> 267</span> <span></span>
<span>/Users/kde/yt-x86_64/src/yt-hg/yt/data_objects/static_output.pyc</span> in <span>__getitem__<span>(self, key)</span>
<span> 146</span> for d in [self.units, self.time_units, self.parameters, \
<span> 147</span> self.conversion_factors]:
<span>--> 148<span> <span>if</span> </span>key</span> <span>in</span> </span>d<span>:</span> <span>return</span> d<span>[</span>key<span>]</span><span></span>
<span> 149</span> <span>raise</span> KeyError<span>(</span>key<span>)</span><span></span>
<span> 150</span> <span></span>
<span>TypeError</span>: unhashable type
</pre><br></div><div><div>thanks<span class="HOEnZb"><font color="#888888"><br>kathy<br><br></font></span></div><span class="HOEnZb"><font color="#888888"><div>-- <br>Kathy DeGioia Eastwood, Ph.D.<br>Professor of Physics and Astronomy<br>
Northern Arizona University<br>Flagstaff, AZ 86011-6010<br>Ph: <a href="tel:928-523-7159" value="+19285237159" target="_blank">928-523-7159</a> FX: <a href="tel:928-523-1371" value="+19285231371" target="_blank">928-523-1371</a><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></font></span></div></div>
<br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>