<div dir="ltr">ok, one issue for Maestro is that there are 2 temperatures.  For low Mach flows, you don't need an energy equation formally, but we carry one in Maestro anyway.  So you can get the temperature with or without the energy (enthalpy actually).  So it would be nice if we could specify which of the temperature fields we want to use with yt (they are named "tfromp" and "tfromh").</div>
<div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Mar 28, 2014 at 8:00 PM, Matthew Turk <span dir="ltr"><<a href="mailto:matthewturk@gmail.com" target="_blank">matthewturk@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 class="">On Fri, Mar 28, 2014 at 7:59 PM, Michael Zingale<br>
<<a href="mailto:michael.zingale@stonybrook.edu">michael.zingale@stonybrook.edu</a>> wrote:<br>
> Thanks.  I see in there<br>
><br>
> def _temperature(field,data):<br>
>     mu = data.get_field_parameter("mu")<br>
>     return ( (data.pf.gamma-1.0) * mu * mh *<br>
>              data["thermal_energy"] / (kboltz * data["density"]) )<br>
><br>
> is it the case that yt will assume a gamma-law gas?<br>
<br>
</div>Only if temperature isn't defined -- the field system is set up now to<br>
be able to dynamically swap these in and out, so we definitely should.<br>
<br>
-Matt<br>
<div><div class="h5"><br>
><br>
><br>
> On Fri, Mar 28, 2014 at 7:53 PM, Nathan Goldbaum <<a href="mailto:nathan12343@gmail.com">nathan12343@gmail.com</a>><br>
> wrote:<br>
>><br>
>><br>
>><br>
>> On Fri, Mar 28, 2014 at 4:52 PM, Michael Zingale<br>
>> <<a href="mailto:michael.zingale@stonybrook.edu">michael.zingale@stonybrook.edu</a>> wrote:<br>
>>><br>
>>> ok, final update -- I am able to do a basic volume rendering, but I need<br>
>>> to set<br>
>>><br>
>>> no_ghost=True<br>
>>><br>
>>> if this is False, then there is a units error in sanitize_units_add.<br>
>>> Perhaps when the ghost-cells are added, that is done as a normal ndarray and<br>
>>> not a YTarray?<br>
>>><br>
>>> Similar issues arise if I try to draw_coordinate_vectors()<br>
>>><br>
>>><br>
>><br>
>> Can you report these issues on our bug tracker?  A test script that fails<br>
>> for you would be especially helpful.<br>
>><br>
>>><br>
>>><br>
>>> On Fri, Mar 28, 2014 at 7:26 PM, Michael Zingale<br>
>>> <<a href="mailto:michael.zingale@stonybrook.edu">michael.zingale@stonybrook.edu</a>> wrote:<br>
>>>><br>
>>>> I made a bit more progress here, but I get an error in<br>
>>>><br>
>>>> amr_kdtree/amr_kdtree.py<br>
>>>><br>
>>>> in line 257, it does<br>
>>>><br>
>>>>         dds = grid.dds.ndarray_view()<br>
>>>><br>
>>>> but grid.dds is already an ndarray, so this should be<br>
>>>><br>
>>>>         dds = grid.dds<br>
>>>><br>
>>>> still not through the cam.snapshot() part of my script, but I got<br>
>>>> everything up to there working.<br>
>>>><br>
>>>><br>
>>>><br>
>>>> On Fri, Mar 28, 2014 at 4:23 PM, Michael Zingale<br>
>>>> <<a href="mailto:michael.zingale@stonybrook.edu">michael.zingale@stonybrook.edu</a>> wrote:<br>
>>>>><br>
>>>>> On Fri, Mar 28, 2014 at 3:37 PM, Nathan Goldbaum<br>
>>>>> <<a href="mailto:nathan12343@gmail.com">nathan12343@gmail.com</a>> wrote:<br>
>>>>>><br>
>>>>>> Hi Mike,<br>
>>>>>><br>
>>>>>> You've hit a bunch of big changes we've been making recently.  If you<br>
>>>>>> want to get back to the version of yt before we merged these big changes in,<br>
>>>>>> you'll need to update to the 'development' bookmark.  Right now you're on<br>
>>>>>> the 'experimental' bookmark.<br>
>>>>>><br>
>>>>><br>
>>>>> Oh... lucky me.<br>
>>>>><br>
>>>>>><br>
>>>>>><br>
>>>>>> On Fri, Mar 28, 2014 at 12:10 PM, Michael Zingale<br>
>>>>>> <<a href="mailto:michael.zingale@stonybrook.edu">michael.zingale@stonybrook.edu</a>> wrote:<br>
>>>>>>><br>
>>>>>>> Is there any updated examples for using the yt 3.0 volume rendering?<br>
>>>>>>><br>
>>>>>>> For instance, looking at this:<br>
>>>>>>><br>
>>>>>>><br>
>>>>>>> <a href="http://yt-project.org/docs/dev/cookbook/simple_plots.html#simple-volume-rendering" target="_blank">http://yt-project.org/docs/dev/cookbook/simple_plots.html#simple-volume-rendering</a><br>

>>>>>><br>
>>>>>><br>
>>>>>> We haven't finished updating the docs yet.  In-development docs<br>
>>>>>> (including some out-of-date information that we need to update) are<br>
>>>>>> available here: <a href="http://yt-project.org/docs/dev-3.0/" target="_blank">http://yt-project.org/docs/dev-3.0/</a><br>
>>>>>><br>
>>>>>>><br>
>>>>>>><br>
>>>>>>> The field name is no longer simply "Density", but now, for my data, I<br>
>>>>>>> need to use a tuple like:<br>
>>>>>>><br>
>>>>>>> field =  ('boxlib', 'magvel')<br>
>>>>>><br>
>>>>>><br>
>>>>>> We no longer have a 'Density' field, since this introduces ambiguities<br>
>>>>>> for Enzo data and to create a consistent naming scheme for fields.  The new<br>
>>>>>> name is 'density'. See:<br>
>>>>>> <a href="http://ytep.readthedocs.org/en/latest/YTEPs/YTEP-0003.html" target="_blank">http://ytep.readthedocs.org/en/latest/YTEPs/YTEP-0003.html</a><br>
>>>>><br>
>>>>><br>
>>>>> thanks for the list.  Maestro has it's own fun with variables (two<br>
>>>>> definitions for temperature, several pressures, etc., magvel above is not<br>
>>>>> simply the vector magnitude of x_vel, y_vel, z_vel, etc), so I'll pay<br>
>>>>> attention to these definitions.<br>
>>>>><br>
>>>>><br>
>>>>><br>
>>>>>>><br>
>>>>>>><br>
>>>>>>> And then when I look at the extrema, I no longer take the [0] index,<br>
>>>>>>> as that just gives the min, but instead if I look at:<br>
>>>>>>><br>
>>>>>><br>
>>>>>> That true, we've simplified the derived quantity API somewhat.  This<br>
>>>>>> fixes a long-standing issue at the cost of breaking some old user scripts.<br>
>>>>>> As you can see we are planning to break backward compatibility somewhat for<br>
>>>>>> the final yt 3.0 release.<br>
>>>>>><br>
>>>>>>><br>
>>>>>>> dd.quantities["Extrema"](field)<br>
>>>>>>><br>
>>>>>>> I get:<br>
>>>>>>><br>
>>>>>>> (124.858613718 dimensionless, 18061918.0546 dimensionless)<br>
>>>>>><br>
>>>>>><br>
>>>>>>><br>
>>>>>>> (not sure why they are dimensionless, but that's a different<br>
>>>>>>> problem).<br>
>>>>>>><br>
>>>>>><br>
>>>>>> If you do<br>
>>>>>><br>
>>>>>> type(dd.quantities["Extrema"](field)[0])<br>
>>>>>><br>
>>>>>> you will see that it is a new type, YTQuantity, which is a new class<br>
>>>>>> we've defined to handle units for scalar data.  YTArray does the same thing<br>
>>>>>> for array data.  More detail here:<br>
>>>>>> <a href="http://yt-project.org/docs/dev-3.0/analyzing/units/index.html" target="_blank">http://yt-project.org/docs/dev-3.0/analyzing/units/index.html</a><br>
>>>>>><br>
>>>>>> These are 'dimensionless' because the boxlib frontend does not know<br>
>>>>>> what the appropriate units are for the 'magvel' field.  It may be the case<br>
>>>>>> that we could modify the frontend to detect the units<br>
>>>>>><br>
>>>>>>><br>
>>>>>>> Finally, if I try to make a transfer function, like:<br>
>>>>>>><br>
>>>>>>> tf = ColorTransferFunction((mi, ma))<br>
>>>>>>><br>
>>>>>>> I get:<br>
>>>>>>><br>
>>>>>>> NameError: name 'ColorTransferFunction' is not defined<br>
>>>>>><br>
>>>>>><br>
>>>>>> I think this was an oversight on Matt's part.  Thanks for the report!<br>
>>>>>> For now you should be able to import ColorTransferFunction from<br>
>>>>>> yt.visualization.volume_rendering.api, but this should be fixed soon.<br>
>>>>>><br>
>>>>>> The reason this happened is that we've changed the way imports work in<br>
>>>>>> the top-level yt namespace in the interest of yt behaving more like a<br>
>>>>>> 'normal' scientific python package.  More detail is available here:<br>
>>>>>><br>
>>>>>> <a href="http://ytep.readthedocs.org/en/latest/YTEPs/YTEP-0019.html" target="_blank">http://ytep.readthedocs.org/en/latest/YTEPs/YTEP-0019.html</a><br>
>>>>>><br>
>>>>>> as well as in the following pull request:<br>
>>>>>><br>
>>>>>> <a href="https://bitbucket.org/yt_analysis/yt/pull-request/729/ytep-0019/diff" target="_blank">https://bitbucket.org/yt_analysis/yt/pull-request/729/ytep-0019/diff</a><br>
>>>>>><br>
>>>>>>><br>
>>>>>>><br>
>>>>>>> I just synced up with the development version today.<br>
>>>>>>><br>
>>>>>>> Mike<br>
>>>>>>><br>
>>>>>>><br>
>>>>>><br>
>>>>>> We're working hard to update the docs and get the code into shape to<br>
>>>>>> release.  If you're ok with dealing with not having fully up-to-date docs,<br>
>>>>>> we'd really appreciate you trying out the changes to the code, letting us<br>
>>>>>> know what you think, and reporting any bugs you find.<br>
>>>>>><br>
>>>>><br>
>>>>> yeah, I'm happy to keep playing with this version if it helps.  The<br>
>>>>> reason I updated today is because Matt told me that he had incorporated the<br>
>>>>> Maestro field info that I sent him.<br>
>>>>><br>
>>>>><br>
>>>>>>><br>
>>>>>>><br>
>>>>>>><br>
>>>>>>> --<br>
>>>>>>> Michael Zingale<br>
>>>>>>> Associate Professor<br>
>>>>>>><br>
</div></div>>>>>>>> Dept. of Physics & Astronomy * Stony Brook University * Stony Brook,<br>
<div class="">>>>>>>> NY 11794-3800<br>
>>>>>>> phone:  <a href="tel:631-632-8225" value="+16316328225">631-632-8225</a><br>
>>>>>>> e-mail: <a href="mailto:Michael.Zingale@stonybrook.edu">Michael.Zingale@stonybrook.edu</a><br>
>>>>>>> web: <a href="http://www.astro.sunysb.edu/mzingale" target="_blank">http://www.astro.sunysb.edu/mzingale</a><br>
>>>>>>><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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
>>>>>>><br>
>>>>>><br>
>>>>>><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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
>>>>>><br>
>>>>><br>
>>>>><br>
>>>>><br>
>>>>> --<br>
>>>>> Michael Zingale<br>
>>>>> Associate Professor<br>
>>>>><br>
</div>>>>>> Dept. of Physics & Astronomy * Stony Brook University * Stony Brook, NY<br>
<div class="">>>>>> 11794-3800<br>
>>>>> phone:  <a href="tel:631-632-8225" value="+16316328225">631-632-8225</a><br>
>>>>> e-mail: <a href="mailto:Michael.Zingale@stonybrook.edu">Michael.Zingale@stonybrook.edu</a><br>
>>>>> web: <a href="http://www.astro.sunysb.edu/mzingale" target="_blank">http://www.astro.sunysb.edu/mzingale</a><br>
>>>><br>
>>>><br>
>>>><br>
>>>><br>
>>>> --<br>
>>>> Michael Zingale<br>
>>>> Associate Professor<br>
>>>><br>
</div>>>>> Dept. of Physics & Astronomy * Stony Brook University * Stony Brook, NY<br>
<div class="">>>>> 11794-3800<br>
>>>> phone:  <a href="tel:631-632-8225" value="+16316328225">631-632-8225</a><br>
>>>> e-mail: <a href="mailto:Michael.Zingale@stonybrook.edu">Michael.Zingale@stonybrook.edu</a><br>
>>>> web: <a href="http://www.astro.sunysb.edu/mzingale" target="_blank">http://www.astro.sunysb.edu/mzingale</a><br>
>>><br>
>>><br>
>>><br>
>>><br>
>>> --<br>
>>> Michael Zingale<br>
>>> Associate Professor<br>
>>><br>
</div>>>> Dept. of Physics & Astronomy * Stony Brook University * Stony Brook, NY<br>
<div class="">>>> 11794-3800<br>
>>> phone:  <a href="tel:631-632-8225" value="+16316328225">631-632-8225</a><br>
>>> e-mail: <a href="mailto:Michael.Zingale@stonybrook.edu">Michael.Zingale@stonybrook.edu</a><br>
>>> web: <a href="http://www.astro.sunysb.edu/mzingale" target="_blank">http://www.astro.sunysb.edu/mzingale</a><br>
>>><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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
>>><br>
>><br>
>><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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
>><br>
><br>
><br>
><br>
> --<br>
> Michael Zingale<br>
> Associate Professor<br>
><br>
</div>> Dept. of Physics & Astronomy * Stony Brook University * Stony Brook, NY<br>
<div class="HOEnZb"><div class="h5">> 11794-3800<br>
> phone:  <a href="tel:631-632-8225" value="+16316328225">631-632-8225</a><br>
> e-mail: <a href="mailto:Michael.Zingale@stonybrook.edu">Michael.Zingale@stonybrook.edu</a><br>
> web: <a href="http://www.astro.sunysb.edu/mzingale" target="_blank">http://www.astro.sunysb.edu/mzingale</a><br>
><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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
><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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div>Michael Zingale</div><div>Associate Professor</div><div><br></div><div>Dept. of Physics & Astronomy • Stony Brook University • Stony Brook, NY 11794-3800</div>
<div><i>phone</i>:  631-632-8225</div><div><i>e-mail</i>: <a href="mailto:Michael.Zingale@stonybrook.edu" target="_blank">Michael.Zingale@stonybrook.edu</a></div><div><i>web</i>: <a href="http://www.astro.sunysb.edu/mzingale" target="_blank">http://www.astro.sunysb.edu/mzingale</a></div>

</div>