<div dir="ltr">Hi Chang-Goo,<div><br></div><div>Thanks for the bug report. I've issued a pull request to fix the issue here:</div><div><br></div><div><a href="https://bitbucket.org/yt_analysis/yt/pull-request/881/fixing-an-issue-with-the-profileplot/diff">https://bitbucket.org/yt_analysis/yt/pull-request/881/fixing-an-issue-with-the-profileplot/diff</a><br>
</div><div><br></div><div>I expect that to be merged soon. Once it's in you should be able to pull the change and update to get the fix.</div><div><br></div><div>-Nathan</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Tue, May 6, 2014 at 10:43 AM, Chang-Goo Kim <span dir="ltr"><<a href="mailto:cgkim@astro.princeton.edu" target="_blank">cgkim@astro.princeton.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 all,<br><br>I tried to draw PDFs using ProfilePlot with "fractional" keyword. But, it returned unit incompatibility. I have looked at the code, and realized that the fractional keyword makes field_data dimensionless, while it tried to return data with units of original data. I think it must be easy to fix. Error message is attached below.<br>
<br></div>Thanks,<br><br>Chang-Goo<br><div><div><div><div><pre><span>---------------------------------------------------------------------------</span>
<span>YTUnitConversionError</span> Traceback (most recent call last)
<span><ipython-input-137-18016774c3d0></span> in <span><module><span>()</span>
<span>----> 1<span> </span>pdf</span><span>=</span></span>ProfilePlot<span>(</span>dd<span>,</span><span>"temp"</span><span>,</span><span>[</span><span>'cell_mass'</span><span>,</span><span>'cell_volume'</span><span>]</span><span>,</span>weight_field<span>=</span>None<span>,</span>fractional<span>=</span>True<span>)</span><span></span>
<span> 2</span> pdf<span>.</span>show<span>(</span><span>)</span><span></span>
<span>/u/cgkim/Sources/yt-x86_64/src/yt-hg/yt/visualization/profile_plotter.pyc</span> in <span>__init__<span>(self, data_source, x_field, y_fields, weight_field, n_bins, accumulation, fractional, label, plot_spec, profiles)</span>
<span> 216</span> </span>self<span>.</span>plot_spec <span>=</span> <span>[</span>self<span>.</span>plot_spec<span>.</span>copy<span>(</span><span>)</span> <span>for</span> p <span>in</span> self<span>.</span>profiles<span>]</span><span></span>
<span> 217</span> <span></span>
<span>--> 218<span> </span>self</span><span>.</span>_setup_plots<span>(</span><span>)</span><span></span>
<span> 219</span> <span></span>
<span> 220</span> <span>def</span> save<span>(</span>self<span>,</span> name<span>=</span>None<span>)</span><span>:</span><span></span>
<span>/u/cgkim/Sources/yt-x86_64/src/yt-hg/yt/visualization/profile_plotter.pyc</span> in <span>_setup_plots<span>(self)</span>
<span> 309</span> </span>self<span>.</span>axes <span>=</span> AxesContainer<span>(</span>self<span>.</span>figures<span>)</span><span></span>
<span> 310</span> <span>for</span> i<span>,</span> profile <span>in</span> enumerate<span>(</span>self<span>.</span>profiles<span>)</span><span>:</span><span></span>
<span>--> 311<span> <span>for</span> </span>field</span><span>,</span> field_data <span>in</span> profile<span>.</span>items<span>(</span><span>)</span><span>:</span><span></span>
<span> 312</span> self.axes[field].plot(np.array(profile.x),
<span> 313</span> np<span>.</span>array<span>(</span>field_data<span>)</span><span>,</span><span></span>
<span>/u/cgkim/Sources/yt-x86_64/src/yt-hg/yt/data_objects/profiles.pyc</span> in <span>items<span>(self)</span>
<span> 854</span> <span></span></span>
<span> 855</span> <span>def</span> items<span>(</span>self<span>)</span><span>:</span><span></span>
<span>--> 856<span> <span>return</span> <span>[</span><span>(</span></span>k</span><span>,</span>self<span>[</span>k<span>]</span><span>)</span> <span>for</span> k <span>in</span> self<span>.</span>field_data<span>.</span>keys<span>(</span><span>)</span><span>]</span><span></span>
<span> 857</span> <span></span>
<span> 858</span> <span>def</span> __iter__<span>(</span>self<span>)</span><span>:</span><span></span>
<span>/u/cgkim/Sources/yt-x86_64/src/yt-hg/yt/data_objects/profiles.pyc</span> in <span>__getitem__<span>(self, field)</span>
<span> 851</span> <span>raise</span> </span>KeyError<span>(</span>field<span>)</span><span></span>
<span> 852</span> <span>else</span><span>:</span><span></span>
<span>--> 853<span> <span>return</span> </span>self</span><span>.</span>field_data<span>[</span>fname<span>]</span><span>.</span>in_units<span>(</span>self<span>.</span>field_units<span>[</span>fname<span>]</span><span>)</span><span></span>
<span> 854</span> <span></span>
<span> 855</span> <span>def</span> items<span>(</span>self<span>)</span><span>:</span><span></span>
<span>/u/cgkim/Sources/yt-x86_64/src/yt-hg/yt/units/yt_array.pyc</span> in <span>in_units<span>(self, units)</span>
<span> 358</span> <span></span></span>
<span> 359</span> """
<span>--> 360<span> </span>new_units</span> <span>=</span> self<span>.</span>_unit_repr_check_same<span>(</span>units<span>)</span><span></span>
<span> 361</span> conversion_factor <span>=</span> self<span>.</span>units<span>.</span>get_conversion_factor<span>(</span>new_units<span>)</span><span></span>
<span> 362</span> <span></span>
<span>/u/cgkim/Sources/yt-x86_64/src/yt-hg/yt/units/yt_array.pyc</span> in <span>_unit_repr_check_same<span>(self, units)</span>
<span> 315</span> <span>if</span> <span>not</span> </span>self<span>.</span>units<span>.</span>same_dimensions_as<span>(</span>units<span>)</span><span>:</span><span></span>
<span> 316</span> raise YTUnitConversionError(
<span>--> 317<span> self.units, self.units.dimensions, units, units.dimensions)
</span><span> 318</span> <span></span></span>
<span> 319</span> <span>return</span> units<span></span>
<span>YTUnitConversionError</span>: Unit dimensionalities do not match. Tried to convert between dimensionless (dim 1) and code_length**3 (dim (length)**3).
</pre></div></div></div><br></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>