<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">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><div>Can you pastebin the full traceback?</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>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><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 class="HOEnZb"><font color="#888888"><br>Yuan</font></span></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></div>