<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 style class=""><div class=""><div class=""><pre><span class="">---------------------------------------------------------------------------</span>
<span class="">YTUnitConversionError</span> Traceback (most recent call last)
<span class=""><ipython-input-137-18016774c3d0></span> in <span class=""><module><span class="">()</span>
<span class="">----> 1<span class=""> </span>pdf</span><span class="">=</span></span>ProfilePlot<span class="">(</span>dd<span class="">,</span><span class="">"temp"</span><span class="">,</span><span class="">[</span><span class="">'cell_mass'</span><span class="">,</span><span class="">'cell_volume'</span><span class="">]</span><span class="">,</span>weight_field<span class="">=</span>None<span class="">,</span>fractional<span class="">=</span>True<span class="">)</span><span class=""></span>
<span class=""> 2</span> pdf<span class="">.</span>show<span class="">(</span><span class="">)</span><span class=""></span>
<span class="">/u/cgkim/Sources/yt-x86_64/src/yt-hg/yt/visualization/profile_plotter.pyc</span> in <span class="">__init__<span class="">(self, data_source, x_field, y_fields, weight_field, n_bins, accumulation, fractional, label, plot_spec, profiles)</span>
<span class=""> 216</span> </span>self<span class="">.</span>plot_spec <span class="">=</span> <span class="">[</span>self<span class="">.</span>plot_spec<span class="">.</span>copy<span class="">(</span><span class="">)</span> <span class="">for</span> p <span class="">in</span> self<span class="">.</span>profiles<span class="">]</span><span class=""></span>
<span class=""> 217</span> <span class=""></span>
<span class="">--> 218<span class=""> </span>self</span><span class="">.</span>_setup_plots<span class="">(</span><span class="">)</span><span class=""></span>
<span class=""> 219</span> <span class=""></span>
<span class=""> 220</span> <span class="">def</span> save<span class="">(</span>self<span class="">,</span> name<span class="">=</span>None<span class="">)</span><span class="">:</span><span class=""></span>
<span class="">/u/cgkim/Sources/yt-x86_64/src/yt-hg/yt/visualization/profile_plotter.pyc</span> in <span class="">_setup_plots<span class="">(self)</span>
<span class=""> 309</span> </span>self<span class="">.</span>axes <span class="">=</span> AxesContainer<span class="">(</span>self<span class="">.</span>figures<span class="">)</span><span class=""></span>
<span class=""> 310</span> <span class="">for</span> i<span class="">,</span> profile <span class="">in</span> enumerate<span class="">(</span>self<span class="">.</span>profiles<span class="">)</span><span class="">:</span><span class=""></span>
<span class="">--> 311<span class=""> <span class="">for</span> </span>field</span><span class="">,</span> field_data <span class="">in</span> profile<span class="">.</span>items<span class="">(</span><span class="">)</span><span class="">:</span><span class=""></span>
<span class=""> 312</span> self.axes[field].plot(np.array(profile.x),
<span class=""> 313</span> np<span class="">.</span>array<span class="">(</span>field_data<span class="">)</span><span class="">,</span><span class=""></span>
<span class="">/u/cgkim/Sources/yt-x86_64/src/yt-hg/yt/data_objects/profiles.pyc</span> in <span class="">items<span class="">(self)</span>
<span class=""> 854</span> <span class=""></span></span>
<span class=""> 855</span> <span class="">def</span> items<span class="">(</span>self<span class="">)</span><span class="">:</span><span class=""></span>
<span class="">--> 856<span class=""> <span class="">return</span> <span class="">[</span><span class="">(</span></span>k</span><span class="">,</span>self<span class="">[</span>k<span class="">]</span><span class="">)</span> <span class="">for</span> k <span class="">in</span> self<span class="">.</span>field_data<span class="">.</span>keys<span class="">(</span><span class="">)</span><span class="">]</span><span class=""></span>
<span class=""> 857</span> <span class=""></span>
<span class=""> 858</span> <span class="">def</span> __iter__<span class="">(</span>self<span class="">)</span><span class="">:</span><span class=""></span>
<span class="">/u/cgkim/Sources/yt-x86_64/src/yt-hg/yt/data_objects/profiles.pyc</span> in <span class="">__getitem__<span class="">(self, field)</span>
<span class=""> 851</span> <span class="">raise</span> </span>KeyError<span class="">(</span>field<span class="">)</span><span class=""></span>
<span class=""> 852</span> <span class="">else</span><span class="">:</span><span class=""></span>
<span class="">--> 853<span class=""> <span class="">return</span> </span>self</span><span class="">.</span>field_data<span class="">[</span>fname<span class="">]</span><span class="">.</span>in_units<span class="">(</span>self<span class="">.</span>field_units<span class="">[</span>fname<span class="">]</span><span class="">)</span><span class=""></span>
<span class=""> 854</span> <span class=""></span>
<span class=""> 855</span> <span class="">def</span> items<span class="">(</span>self<span class="">)</span><span class="">:</span><span class=""></span>
<span class="">/u/cgkim/Sources/yt-x86_64/src/yt-hg/yt/units/yt_array.pyc</span> in <span class="">in_units<span class="">(self, units)</span>
<span class=""> 358</span> <span class=""></span></span>
<span class=""> 359</span> """
<span class="">--> 360<span class=""> </span>new_units</span> <span class="">=</span> self<span class="">.</span>_unit_repr_check_same<span class="">(</span>units<span class="">)</span><span class=""></span>
<span class=""> 361</span> conversion_factor <span class="">=</span> self<span class="">.</span>units<span class="">.</span>get_conversion_factor<span class="">(</span>new_units<span class="">)</span><span class=""></span>
<span class=""> 362</span> <span class=""></span>
<span class="">/u/cgkim/Sources/yt-x86_64/src/yt-hg/yt/units/yt_array.pyc</span> in <span class="">_unit_repr_check_same<span class="">(self, units)</span>
<span class=""> 315</span> <span class="">if</span> <span class="">not</span> </span>self<span class="">.</span>units<span class="">.</span>same_dimensions_as<span class="">(</span>units<span class="">)</span><span class="">:</span><span class=""></span>
<span class=""> 316</span> raise YTUnitConversionError(
<span class="">--> 317<span class=""> self.units, self.units.dimensions, units, units.dimensions)
</span><span class=""> 318</span> <span class=""></span></span>
<span class=""> 319</span> <span class="">return</span> units<span class=""></span>
<span class="">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>