<div dir="ltr"><div>Traceback (most recent call last):</div><div> File "logP.py", line 11, in <module></div><div> a=dd["gas","logP_gradient_y"]</div><div> File "/Users/yuanli/source/yt3/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 253, in __getitem__</div><div> self.get_data(f)</div><div> File "/Users/yuanli/source/yt3/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 751, in get_data</div><div> self._generate_fields(fields_to_generate)</div><div> File "/Users/yuanli/source/yt3/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 771, in _generate_fields</div><div> fd = self._generate_field(field)</div><div> File "/Users/yuanli/source/yt3/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 290, in _generate_field</div><div> tr = self._generate_fluid_field(field)</div><div> File "/Users/yuanli/source/yt3/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 308, in _generate_fluid_field</div><div> rv = self._generate_spatial_fluid(field, ngt_exception.ghost_zones)</div><div> File "/Users/yuanli/source/yt3/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 318, in _generate_spatial_fluid</div><div> rv = self.ds.arr(np.empty(self.ires.size, dtype="float64"), units)</div><div> File "/Users/yuanli/source/yt3/yt-x86_64/src/yt-hg/yt/units/yt_array.py", line 332, in __new__</div><div> units = Unit(input_units, registry=registry)</div><div> File "/Users/yuanli/source/yt3/yt-x86_64/src/yt-hg/yt/units/unit_object.py", line 204, in __new__</div><div> unit_data = _get_unit_data_from_expr(unit_expr, registry.lut)</div><div> File "/Users/yuanli/source/yt3/yt-x86_64/src/yt-hg/yt/units/unit_object.py", line 496, in _get_unit_data_from_expr</div><div> unit_data = _get_unit_data_from_expr(expr, unit_symbol_lut)</div><div> File "/Users/yuanli/source/yt3/yt-x86_64/src/yt-hg/yt/units/unit_object.py", line 504, in _get_unit_data_from_expr</div><div> "objects." % str(unit_expr))</div><div>yt.units.unit_object.UnitParseError: Cannot parse for unit data from '()'. Please supply an expression of only Unit, Symbol, Pow, and Mulobjects.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 13, 2015 at 3:28 PM, Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</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"><br><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Thu, Aug 13, 2015 at 2:22 PM, Yuan Li <span dir="ltr"><<a href="mailto:yuan@astro.columbia.edu" target="_blank">yuan@astro.columbia.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">Hi all:<br><br>I am having trouble creating a gradient field in yt3. What I want is the gradient of logP, where logP is defined as:<br><br>@derived_field(name="logP",units="")<br>def _logP(field,data):<br> return numpy.log(numpy.array(data["gas","pressure"]))<br><br>Then I add the gradient field:<br><br>pf.add_gradient_fields(("gas","logP"))</div></blockquote><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br>which seems to work and I can find it in pf.derived_field_list.<br><br>But if I try to use it, I get this error message:<br>yt.units.unit_object.UnitParseError: Cannot parse for unit data from '()'. Please supply an expression of only Unit, Symbol, Pow, and Mulobjects.<br></div></blockquote><div><br></div></span><div>Can you pastebin the full traceback?</div><span class=""><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br>Is this because logP does not have a unit? <br><br>Also, I had another unrelated question about the gradient fields: <br>Can I add the gradient fields globally instead of to a specific dataset? <br><br></div></blockquote><div><br></div></span><div>Not right now, no.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br>Thanks!<span><font color="#888888"><br>Yuan</font></span></div>
<br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></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" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>