<div dir="ltr">Hi Madison,<div><br></div><div>When using the data object within a derived field, you have to specify what field you want to access. So line 13 should look like </div><div><br></div><div><pre style="margin-top:0px;margin-bottom:0px;padding:5px 0px;font-family:'Bitstream Vera Sans Mono',monospace;font-size:13px;color:rgb(0,0,0);line-height:normal"><span class="n">new_field</span> <span class="o" style="color:rgb(102,102,102)">=</span> <span class="n">np</span><span class="o" style="color:rgb(102,102,102)">.</span><span class="n">maximum</span><span class="p">(</span><span class="n">data['radial_velocity']</span><span class="p">,</span> <span class="mf" style="color:rgb(64,160,112)">0.0</span><span class="p">)</span></pre><pre style="margin-top:0px;margin-bottom:0px;padding:5px 0px;font-family:'Bitstream Vera Sans Mono',monospace;font-size:13px;color:rgb(0,0,0);line-height:normal"><span class="p"><br></span></pre>Right now the field is trying to find the maximum of the YTSphere object rather than the radial velocity data contained within the YTSphere, which isn't well defined. </div><div><br></div><div>One other quick comment, once that's fixed your script is going to try and compare data['radial_velocity'] which has units of velocity, and 0.0 which is unitless. In order for that to be a valid comparison, you should give 0.0 in the appropriate units by doing something like this:</div><div><br></div><div><pre style="font-size:13px;margin-top:0px;margin-bottom:0px;padding:5px 0px;font-family:'Bitstream Vera Sans Mono',monospace;color:rgb(0,0,0);line-height:normal">zero_vel = data.ds.quan(0.0, 'km/s')</pre><pre style="font-size:13px;margin-top:0px;margin-bottom:0px;padding:5px 0px;font-family:'Bitstream Vera Sans Mono',monospace;line-height:normal"><span style="color:rgb(0,0,0)"><span class="n">new_field</span> <span class="o" style="color:rgb(102,102,102)">=</span> <span class="n">np</span><span class="o" style="color:rgb(102,102,102)">.</span><span class="n">maximum</span><span class="p">(</span><span class="n">data['radial_velocity']</span><span class="p">, zero_vel</span></span><span class="p" style="color:rgb(0,0,0)">)</span></pre></div><div><br></div>Hope that helps!<div>-Hilary</div><div><br><div class="gmail_quote"><div dir="ltr">On Wed, Jul 22, 2015 at 8:49 AM Madison Fitzgerald <<a href="mailto:fitzg152@msu.edu">fitzg152@msu.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>Hello,</div><div><br></div><div>I am trying to create two new derived fields and then make a simple phase plot using one of these fields, but I am running into a pretty consistent error involving “convert_to_units". I get basically the same error when using the ds.sphere object and a ds.all_data object, so if there’s something I’ve messed up on, it’s not obvious to me and I would appreciate any help figuring out what’s wrong.</div><div><br></div><div>I am running yt on a MacBook Pro with OS 10.10.4. I have the newest version of stable yt (3.1 and changeset 665e2b6cdf78), and you can find my script at this link (<a href="http://paste.yt-project.org/show/5733/" target="_blank">http://paste.yt-project.org/show/5733/</a>) and my error output at this link (<a href="http://paste.yt-project.org/show/5734/" target="_blank">http://paste.yt-project.org/show/5734/</a>).</div><div><br></div><div>Thank you very much,</div><div>Madison Fitzgerald</div><div><br></div><div><br></div><br><div>
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div>Madison Fitzgerald</div><div><br></div><div>Astrophysics | Lyman Briggs College</div><div>Women’s & Gender Studies | College of Arts & Letters</div><div>LGBTQ and Sexuality Studies Specialization</div><div><div>Women in Science, President</div><div>Undergraduate Teaching Assistant</div><div>Undergraduate Research Assistant</div><div>Department of Physics and Astronomy</div></div><div>Honors College</div><div>Michigan State University</div></div></div></div></div></div></div></div>
</div>
<br></div>_______________________________________________<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>
</blockquote></div></div></div>