<div dir="ltr">Hi Molly,<div><br></div><div>I think this is a bug in yt - the issue is that the field_info associated with the original profile object is not being saved to the profile dataset.</div><div><br></div><div>I've filed this issue here:</div><div><br></div><div><a href="https://bitbucket.org/yt_analysis/yt/issues/1267">https://bitbucket.org/yt_analysis/yt/issues/1267</a><br></div><div><br></div><div>For now, you should be able to create a PhasePlot using the profile object you created from the original dataset:</div><div><br></div><div><div> nhT_prof_new = yt.create_profile(sph, [('gas', 'H_number_density'),('gas','temperature')], fields='cell_mass',</div><div> extrema=extrema,fractional=True, n_bins=6)</div><div> plot = yt.PhasePlot.from_profile(nht_prof_new)</div></div><div><br></div><div>-Nathan</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Aug 29, 2016 at 12:45 PM, Molly Peeples <span dir="ltr"><<a href="mailto:molly@stsci.edu" target="_blank">molly@stsci.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">
Hi folks,
<div><br>
</div>
<div>I’m trying to construct and save a 2d profile container to make a phase plot. Here’s what I have:</div>
<div>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
nhT_prof_new = yt.create_profile(sph, [('gas', 'H_number_density'),('gas','<wbr>temperature')], fields='cell_mass', extrema=extrema,fractional=<wbr>True, n_bins=6)</p>
</div>
<div>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
nhT_prof_new.save_as_dataset()</p>
</div>
<div>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
prof_ds = yt.load("snap_047_z000p000_<wbr>Profile2D.h5")</p>
</div>
<div><br>
</div>
<div>(just using a small number of bins for now to play with…)</div>
<div><br>
</div>
<div>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#35a327">In [</span><span style="color:#2de621">115</span><span style="color:#35a327">]:
</span>prof_ds.data[<span style="color:#ce7924">"temperature"</span>]</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(206,121,36);background-color:rgb(0,0,0)">
<span style="color:#35a327"> ...: </span><span style="color:#f5f5f5">prof_ds.data[</span>"H_number_<wbr>density"<span style="color:#f5f5f5">]</span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#35a327"> ...: </span>prof_ds.data[<span style="color:#ce7924">"cell_mass"</span>]</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(53,163,39);background-color:rgb(0,0,0)">
...: </p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(178,54,34);background-color:rgb(0,0,0)">
Out[<span style="color:#ff3b1d">115</span>]: </p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
YTArray([[ 0.01685665, 0.0405085 , 0.03575519, 0.03091732, 0.02552936,</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
0.0230379 ],</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
[ 0.03770747, 0.05441198, 0.04859202, 0.03285325, 0.02385346,</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
0.04674817],</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
[ 0.03022001, 0.05150457, 0.03591368, 0.01974481, 0.01052408,</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
0.00331175],</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
[ 0.04278126, 0.01889035, 0.0177686 , 0.01644028, 0.01598129,</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
0.00578262],</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
[ 0.05010083, 0.06573628, 0.04247398, 0.03914208, 0. ,</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
0. ],</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
[ 0. , 0.11691224, 0. , 0. , 0. ,</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
0. ]]) dimensionless</p>
</div>
<div><br>
</div>
<div>but when I try to create a PhasePlot:</div>
<div>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#35a327">In [</span><span style="color:#2de621">116</span><span style="color:#35a327">]:
</span>p = yt.PhasePlot(prof_ds.data, <span style="color:#ce7924">"H_number_density"</span>,
<span style="color:#ce7924">"temperature"</span>,<span style="color:#ce7924">"cell_mass"</span>,<wbr>weight_field=<span style="color:#35a327">None</span>)</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(53,163,39);background-color:rgb(0,0,0)">
...: </p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(195,55,32);background-color:rgb(0,0,0)">
------------------------------<wbr>------------------------------<wbr>---------------</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#c33720">KeyError</span> Traceback (most recent call last)</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(52,189,38);background-color:rgb(0,0,0)">
<ipython-input-116-<wbr>37b81d75b563><span style="color:#f5f5f5"> in </span><span style="color:#34bbc7"><module></span><span style="color:#5330e1">()</span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(83,48,225);background-color:rgb(0,0,0)">
<span style="color:#34bd26">----> 1</span><span style="color:#c33720"> </span><span style="color:#f5f5f5">p
</span>=<span style="color:#f5f5f5"> yt</span>.<span style="color:#f5f5f5">PhasePlot</span>(<span style="color:#f5f5f5">prof_ds</span>.<span style="color:#f5f5f5">data</span>,<span style="color:#f5f5f5">
</span>"H_number_density",<span style="color:#f5f5f5"> </span>"temperature","cell_mass",<span style="color:#f5f5f5">weig<wbr>ht_field</span>=<span style="color:#f5f5f5">None</span>)</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0);min-height:14px">
<br>
</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(83,48,225);background-color:rgb(0,0,0)">
<span style="color:#34bd26">/Users/molly/anaconda2/envs/<wbr>astroconda/lib/python2.7/site-<wbr>packages/yt/visualization/<wbr>profile_plotter.pyc</span><span style="color:#f5f5f5"> in
</span><span style="color:#34bbc7">__init__</span>(self, data_source, x_field, y_field, z_fields, weight_field, x_bins, y_bins, accumulation, fractional, fontsize, figure_size)</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(45,230,33);background-color:rgb(0,0,0)">
756<span style="color:#f5f5f5"> </span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#2de621"> 757</span> type(self)._initialize_<wbr>instance(self, data_source, profile, fontsize,</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(195,55,32);background-color:rgb(0,0,0)">
<span style="color:#34bd26">--> 758</span> figure_size)</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(45,230,33);background-color:rgb(0,0,0)">
759<span style="color:#f5f5f5"> </span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#2de621"> 760</span> <span style="color:#5330e1">@</span>classmethod</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0);min-height:14px">
<br>
</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(52,189,38);background-color:rgb(0,0,0)">
/Users/molly/anaconda2/envs/<wbr>astroconda/lib/python2.7/site-<wbr>packages/yt/visualization/<wbr>profile_plotter.pyc<span style="color:#f5f5f5"> in
</span><span style="color:#34bbc7">_initialize_instance</span><span style="color:#5330e1">(cls, obj, data_source, profile, fontsize, figure_size)</span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#2de621"> 773</span> obj<span style="color:#5330e1">.</span>profile
<span style="color:#5330e1">=</span> profile</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#2de621"> 774</span> super<span style="color:#5330e1">(</span>PhasePlot<span style="color:#5330e1">,</span> obj<span style="color:#5330e1">).</span>__init__<span style="color:#5330e1">(</span>data_source<span style="color:#5330e1">,</span>
figure_size<span style="color:#5330e1">,</span> fontsize<span style="color:#5330e1">)</span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#34bd26">--> 775</span><span style="color:#c33720"> </span>
obj<span style="color:#5330e1">.</span>_setup_plots<span style="color:#5330e1">()</span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#2de621"> 776</span> obj<span style="color:#5330e1">.</span>_initfinished
<span style="color:#5330e1">=</span> True</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#2de621"> 777</span> <span style="color:#34bd26">
return</span> obj</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0);min-height:14px">
<br>
</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(52,189,38);background-color:rgb(0,0,0)">
/Users/molly/anaconda2/envs/<wbr>astroconda/lib/python2.7/site-<wbr>packages/yt/visualization/<wbr>profile_plotter.pyc<span style="color:#f5f5f5"> in
</span><span style="color:#34bbc7">_setup_plots</span><span style="color:#5330e1">(self)</span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#2de621"> 854</span> cax <span style="color:#5330e1">
=</span> self<span style="color:#5330e1">.</span>plots<span style="color:#5330e1">[</span>f<span style="color:#5330e1">].</span>cax</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(45,230,33);background-color:rgb(0,0,0)">
855<span style="color:#f5f5f5"> </span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#34bd26">--> 856</span><span style="color:#c33720">
</span>x_scale<span style="color:#5330e1">,</span> y_scale<span style="color:#5330e1">,</span> z_scale
<span style="color:#5330e1">=</span> self<span style="color:#5330e1">.</span>_get_field_log<span style="color:#5330e1">(</span>f<span style="color:#5330e1">,</span> self<span style="color:#5330e1">.</span>profile<span style="color:#5330e1">)</span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#2de621"> 857</span> x_title<span style="color:#5330e1">,</span> y_title<span style="color:#5330e1">,</span> z_title
<span style="color:#5330e1">=</span> self<span style="color:#5330e1">.</span>_get_field_title<span style="color:#5330e1">(</span>f<span style="color:#5330e1">,</span> self<span style="color:#5330e1">.</span>profile<span style="color:#5330e1">)</span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(45,230,33);background-color:rgb(0,0,0)">
858<span style="color:#f5f5f5"> </span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0);min-height:14px">
<br>
</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(52,189,38);background-color:rgb(0,0,0)">
/Users/molly/anaconda2/envs/<wbr>astroconda/lib/python2.7/site-<wbr>packages/yt/visualization/<wbr>profile_plotter.pyc<span style="color:#f5f5f5"> in
</span><span style="color:#34bbc7">_get_field_log</span><span style="color:#5330e1">(self, field_z, profile)</span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(45,230,33);background-color:rgb(0,0,0)">
813<span style="color:#f5f5f5"> </span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#2de621"> 814</span> <span style="color:#34bd26">def</span> _get_field_log<span style="color:#5330e1">(</span>self<span style="color:#5330e1">,</span> field_z<span style="color:#5330e1">,</span> profile<span style="color:#5330e1">):</span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#34bd26">--> 815</span><span style="color:#c33720"> </span>
zfi <span style="color:#5330e1">=</span> profile<span style="color:#5330e1">.</span>field_info<span style="color:#5330e1">[</span>field_z<span style="color:#5330e1">]</span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#2de621"> 816</span> <span style="color:#34bd26">
if</span> self<span style="color:#5330e1">.</span>x_log <span style="color:#34bd26">
is</span> None<span style="color:#5330e1">:</span></p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#2de621"> 817</span> x_log <span style="color:#5330e1">
=</span> profile<span style="color:#5330e1">.</span>x_log</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0);min-height:14px">
<br>
</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)">
<span style="color:#c33720">KeyError</span>: ('data', 'cell_mass')</p>
<p style="margin:0px;font-size:10px;font-family:Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0);min-height:14px">
<br>
</p>
</div>
<div><br>
</div>
<div><br>
</div>
<div>Suggestions?</div><span class="HOEnZb"><font color="#888888">
<div>—Molly</div>
</font></span></div>
<br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
<br></blockquote></div><br></div>