<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 11, 2016 at 12:47 PM, Chang-Goo Kim <span dir="ltr"><<a href="mailto:cgkim@astro.princeton.edu" target="_blank">cgkim@astro.princeton.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div>Thanks, there's no hurry for this specific case. <br><br></div>I have another problem with magnetic field units. I set unit conversion factors with 'unit_override' keyword. These are the units I'm using<br><br><div style="margin-left:40px">>>> print ya.unit_base<br><br>{'length_unit': (1.0, 'pc'),<br> 'mass_unit': (2.38858753789e-24, 'g/cm**3*pc**3'),<br> 'time_unit': (1.0, 's*pc/km')}<br></div><br></div>This should set the magnetic field unit as follows:<br><div><br><div style="margin-left:40px">>>> mu=yt.YTQuantity(ya.unit_base[<wbr>'mass_unit'][0],ya.unit_base['<wbr>mass_unit'][1])<br>>>> lu=yt.YTQuantity(ya.unit_base[<wbr>'length_unit'][0],ya.unit_<wbr>base['length_unit'][1])<br>>>> tu=yt.YTQuantity(ya.unit_base[<wbr>'time_unit'][0],ya.unit_base['<wbr>time_unit'][1])<br>>>> mag_unit=(np.sqrt(4*np.pi*mu/<wbr>lu)/tu)<br>>>> print mag_unit.convert_to_units('<wbr>gauss')<br><br>5.4786746797e-07 gauss<br></div><br></div><div>But, what I found in the unit_registry is a bit strange number.<br></div><div><br><div style="margin-left:40px">>>> print ds.unit_registry['code_<wbr>magnetic']<br><br>(3.5449077018110318,<br> sqrt((mass))/(sqrt((length))*<wbr>(time)),<br> 0.0,<br> '\\rm{code\\ magnetic}')<br></div><br></div><div>It seems that "magnetic_unit" was set correctly in yt/frontends/athena/data_<wbr>structures.py:<br><br>        self.magnetic_unit = np.sqrt(4*np.pi * self.mass_unit /<br>                              <wbr>       (self.time_unit**2 * self.length_unit))<br>        self.magnetic_unit.convert_to_<wbr>units("gauss")<br><br>I don't know what's wrong here. Do you have any idea?<br></div></div></blockquote><div><br></div><div>Sounds like a bug. It might be that the overriden units aren't applied until after the magnetic_unit is calculated.</div><div><br></div><div>Would you mind opening an issue to track this?</div><div><br></div><div><a href="https://bitbucket.org/yt_analysis/yt/issues/new">https://bitbucket.org/yt_analysis/yt/issues/new</a><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><br></div><div>Thanks,<br></div><div>Chang-Goo<br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Thu, Aug 11, 2016 at 12:12 PM, John Zuhone <span dir="ltr"><<a href="mailto:jzuhone@gmail.com" target="_blank">jzuhone@gmail.com</a>></span> wrote:<br></span><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">



<div style="word-wrap:break-word"><span class=""><span>
We can add in a particle_file keyword argument like we do for FLASH data. 
<div><br>
</div>
<div>I can help out with this, but I might not be able to get to it until early next week. </div>
</span></span><div><br>
<div>
<blockquote type="cite"><span class=""><span>
<div>On Aug 11, 2016, at 12:11 PM, Nathan Goldbaum <<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a>> wrote:</div>
<br>
</span></span><div><div class="h5"><div><div><div>
<div dir="ltr" style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Thu, Aug 11, 2016 at 11:03 AM, Chang-Goo Kim<span> </span><span dir="ltr"><<a href="mailto:cgkim@astro.princeton.edu" target="_blank">cgkim@astro.princeton.edu</a><wbr>></span><span> </span>wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div dir="ltr">
<div>
<div>Hi Nathan,<br>
<br>
</div>
Thanks for the reply. Adding particles to the athena frontend will not be easy since there is no unified (and well settled) data dump for particles. Star particles are not in the public release. Let me find other ways at this point.<br>
<br>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>Ah, in that case we should probably add a way for you to supply your own particle data. You won't get chunked particle IO but if that's sufficient it should probably be straightforward to load the data in as a keyword argument for the load() command.</div>
<div><br>
</div>
<div>John might have some insight on what's best to do here, he's the maintainer of the athena frontend.</div>
<div><br>
</div>
<div>-Nathan</div>
<div> </div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div dir="ltr">
<div></div>
<div>Thanks,<br>
</div>
Chang-Goo<br>
<div><br>
</div>
</div>
<div>
<div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Thu, Aug 11, 2016 at 11:35 AM, Nathan Goldbaum<span> </span><span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.co<wbr>m</a>></span><span> </span>wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div dir="ltr">I think we'd need to add support for reading particle data to the athena frontend. I think the main reason this hasn't been done yet is lack of available test datasets. If you'd like to make one or more test datasets available to us
 (they will live on<span> </span><a href="http://yt-project.org/data" target="_blank">yt-project.org/data</a><span> </span>so we can use them for testing and debugging going forward), someone
 could try to add support for reading particle data to the Athena frontend.
<div><br>
</div>
<div>To answer your question, no, there isn't a way to add particle fields to an athena dataset without modifying the athena frontend.</div>
<div><br>
</div>
<div>That said, you might be able to cheat and load in the particle fields using the stream frontend to reload your data as an in-memory dataset - in particular load_uniform_grid or load_amr_grids. Of course this is sort of a hack, ideally we'd just
 add particle support to the athena frontend.</div>
<div><br>
</div>
<div>-Nathan<span><font color="#888888"><br>
<div><br>
</div>
<div>-Nathan</div>
</font></span></div>
</div>
<div>
<div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Thu, Aug 11, 2016 at 10:32 AM, Chang-Goo Kim<span> </span><span dir="ltr"><<a href="mailto:cgkim@astro.princeton.edu" target="_blank">cgkim@astro.princeton.edu</a><wbr>></span><span> </span>wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div dir="ltr">
<div>
<div>Sorry for confusion.<br>
<br>
</div>
I have a separate vtk file for star particles. Just wonder whether I can add this to original dataset as fields to handle them together.<br>
<br>
</div>
<div>Thanks,<br>
</div>
Chang-Goo<br>
</div>
<div>
<div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Wed, Aug 10, 2016 at 3:19 PM, John Zuhone<span> </span><span dir="ltr"><<a href="mailto:jzuhone@gmail.com" target="_blank">jzuhone@gmail.com</a>></span><span> </span>wro<wbr>te:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
Hi Chang-Goo,<br>
<span><br>
> When I load the dataset from Athena and I have star particle mass, position, velocity, age information separately,<br>
<br>
</span>Sorry, I’m a bit confused as to what you mean by this, because as far as I know we don’t have particle support for Athena data yet. Do you actually see the particle data when you load the dataset? Also, is this the old Athena or Athena++?<br>
<span><br>
> how can I add particle information to yt dataset.<br>
<br>
</span>I think by this you mean getting the particle fields into yt, correct? Just making sure.<br>
<br>
Is the data part of the Athena file itself or is it in a separate file?<br>
<br>
Best,<br>
<br>
John<br>
<div>
<div><br>
> On Aug 10, 2016, at 2:24 PM, Chang-Goo Kim <<a href="mailto:changgoo@princeton.edu" target="_blank">changgoo@princeton.edu</a>> wrote:<br>
><br>
> Hi all,<br>
><br>
> I'm analysing data from Athena simulations with star particles. When I load the dataset from Athena and I have star particle mass, position, velocity, age information separately, how can I add particle information to yt dataset. I tried to look up the hep,
 but I cannot find any good reference for this issue.<br>
><br>
> Thanks,<br>
> Chang-Goo<br>
</div>
</div>
> ______________________________<wbr>_________________<br>
> yt-users mailing list<br>
><span> </span><a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
><span> </span><a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/l<wbr>istinfo.cgi/yt-users-spacepope<wbr>.org</a><br>
<br>
______________________________<wbr>_________________<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/lis<wbr>tinfo.cgi/yt-users-spacepope.o<wbr>rg</a><br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
<br>
______________________________<wbr>_________________<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/lis<wbr>tinfo.cgi/yt-users-spacepope.o<wbr>rg</a><br>
<br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
<br>
______________________________<wbr>_________________<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/lis<wbr>tinfo.cgi/yt-users-spacepope.o<wbr>rg</a><br>
<br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
<br>
______________________________<wbr>_________________<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/lis<wbr>tinfo.cgi/yt-users-spacepope.o<wbr>rg</a><br>
<br>
</blockquote>
</div>
<br>
</div>
</div>
<span style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">______________________________<wbr>_________________</span><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
<span style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">yt-users
 mailing list</span><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
<a href="mailto:yt-users@lists.spacepope.org" style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px" target="_blank">yt-users@lists.spacepope.org</a><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px" target="_blank">http://lists.spacepope.org/lis<wbr>tinfo.cgi/yt-users-spacepope.<wbr>org</a></div>
</div></div></div></div></blockquote>
</div>
<br>
</div>
</div>

</blockquote></div><br></div>
<br>______________________________<wbr>_________________<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/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
<br></blockquote></div><br></div></div>