<div dir="ltr">Hi Madison,<div><br></div><div>For whatever reason, yt doesn't convert units that have the same dimensions within maximum so you'll have to manually put your zero velocity in the same units as the data['radial_velocity']. You can either do this:</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>zero_vel</span> <span style="color:rgb(102,102,102)">=</span> <span>data</span><span style="color:rgb(102,102,102)">.</span><span>ds</span><span style="color:rgb(102,102,102)">.</span><span>quan</span><span>(</span><span style="color:rgb(64,160,112)">0.0</span><span>,</span> <span style="color:rgb(64,112,160)">'cm/s'</span><span>)</span>
<span>new_field</span> <span style="color:rgb(102,102,102)">=</span> <span>np</span><span style="color:rgb(102,102,102)">.</span><span>maximum</span><span>(</span><span>data</span><span>[</span><span style="color:rgb(64,112,160)">'radial_velocity'</span><span>],</span> <span>zero_vel.in_units('code_velocity)</span><span>)</span></pre>or since you're just comparing to a value of zero:</div><div><br></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"><span>zero_vel</span> <span style="color:rgb(102,102,102)">=</span> <span>data</span><span style="color:rgb(102,102,102)">.</span><span>ds</span><span style="color:rgb(102,102,102)">.</span><span>quan</span><span>(</span><span style="color:rgb(64,160,112)">0.0</span><span>,</span> <span style="color:rgb(64,112,160)">'code_velocity'</span><span>)</span></pre><div dir="ltr"><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"><span>new_field</span> <span style="color:rgb(102,102,102)">=</span> <span>np</span><span style="color:rgb(102,102,102)">.</span><span>maximum</span><span>(</span><span>data</span><span>[</span><span style="color:rgb(64,112,160)">'radial_velocity'</span><span>],</span> <span>zero_vel</span><span>)</span></pre></div><div dir="ltr"><br><br>Not sure what you mean by it reverses the order of the problem, but hopefully that helps!</div><div dir="ltr"><br></div><div dir="ltr">-Hilary</div><div dir="ltr"><br><div class="gmail_quote"><div dir="ltr">On Wed, Jul 22, 2015 at 10:05 AM Madison Fitzgerald <<a href="mailto:fitzg152@msu.edu" target="_blank">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>Hi Hilary and other yt users,</div><div><br></div>It solved my convert_to_units problem! However, it does raise the error you said it would, the unit comparison problem, and even using your suggestion, I’m not able to resolve the issue. I’ve also tried using the suggestions from the FAQ (<a href="http://yt-project.org/docs/3.1/faq/index.html?highlight=convert%20units#how-do-i-get-the-convert-between-code-units-and-physical-units-for-my-dataset" target="_blank">http://yt-project.org/docs/3.1/faq/index.html?highlight=convert%20units#how-do-i-get-the-convert-between-code-units-and-physical-units-for-my-dataset</a>), but when I do that, it then reverses the order of the problem (saying code_velocity and then cm/s in my error output).<div><br></div><div>Updated code (<a href="http://paste.yt-project.org/show/5736/" target="_blank">http://paste.yt-project.org/show/5736/</a>) and updated error message (<a href="http://paste.yt-project.org/show/5735/" target="_blank">http://paste.yt-project.org/show/5735/</a>).</div><div><br></div><div>Thank you,</div><div>Madison</div><div><br><div></div></div></div><div style="word-wrap:break-word"><div><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></div></div><div style="word-wrap:break-word"><div><div><div><blockquote type="cite"><div>On Jul 22, 2015, at 11:01 AM, Hilary Egan <<a href="mailto:hilaryye@gmail.com" target="_blank">hilaryye@gmail.com</a>> wrote:</div><br><div><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;line-height:normal"><span>new_field</span> <span style="color:rgb(102,102,102)">=</span> <span>np</span><span style="color:rgb(102,102,102)">.</span><span>maximum</span><span>(</span><span>data['radial_velocity']</span><span>,</span> <span style="color:rgb(64,160,112)">0.0</span><span>)</span></pre><pre style="margin-top:0px;margin-bottom:0px;padding:5px 0px;font-family:'Bitstream Vera Sans Mono',monospace;font-size:13px;line-height:normal"><span><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;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><span>new_field</span> <span style="color:rgb(102,102,102)">=</span> <span>np</span><span style="color:rgb(102,102,102)">.</span><span>maximum</span><span>(</span><span>data['radial_velocity']</span><span>, zero_vel</span></span><span>)</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" target="_blank">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="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="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>
_______________________________________________<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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br></div></blockquote></div><br></div></div></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>