<div dir="ltr">Dear Britton<div>Thank you very much.  </div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Feb 22, 2017 at 9:02 PM, Britton Smith <span dir="ltr"><<a href="mailto:brittonsmith@gmail.com" target="_blank">brittonsmith@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 dir="ltr">Hi Reju,<div><br></div><div>When you create a halo catalog in the above manner, the field will be called "particle_position".  This dataset will also be saved to disk after the call to hc.create(), so you can use yt.load to reload the new catalog without having to remake it.</div><div><br></div><div>Britton</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Tue, Feb 21, 2017 at 11:25 PM, Reju Sam John <span dir="ltr"><<a href="mailto:rejusamjohn@gmail.com" target="_blank">rejusamjohn@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Dear All,<div><br></div><div>I am very new in YT 3.3's HaloCatalog functionality. With the following piece of code, I am able to get paerical mass and virial radius. </div><div><font face="monospace, monospace" size="4"><br></font></div><div><font face="monospace, monospace" size="4" color="#ff00ff">import yt</font><div><font face="monospace, monospace" size="4" color="#ff00ff">from yt.analysis_modules.halo_analy<wbr>sis.api import HaloCatalog</font></div><div><font face="monospace, monospace" size="4" color="#ff00ff">data_ds = yt.load('/run/media/john/Seaga<wbr>te_Expansion_Drive/cosmo-sim_<wbr>20/RD0057/RedshiftOutput0057')</font></div><div><font face="monospace, monospace" size="4" color="#ff00ff">hc = HaloCatalog(data_ds=data_ds, finder_method='hop')</font></div><div><font face="monospace, monospace" size="4" color="#ff00ff">hc.create()</font></div><div><font face="monospace, monospace" size="4" color="#ff00ff">ad = hc.halos_ds.all_data()</font></div><div><font face="monospace, monospace" size="4" color="#ff00ff">pm = ad['particle_mass'][:]</font></div><div><div><font face="monospace, monospace" size="4" color="#ff00ff">vr = ad['virial_radius'][:]</font></div><div><br></div></div><div>I would like to calculate the center of mass of the halo. I tried with the following line <br><font face="monospace, monospace" size="4" color="#ff00ff">com = ad['center_of_mass'][:]</font></div><div>since like <font face="monospace, monospace" size="4"><font color="#ff00ff">virial_radius()</font>, <font color="#ff00ff">center_of_mas<wbr>s()</font> </font>was available in YT2.6.</div><div>But this is not working. </div><div><br></div><div>Could you please suggest me how to calculate the center of mass of a halo in yt 3.3.3 ?</div><span class="m_966475192198582347HOEnZb"><font color="#888888"><div><br></div><div><br></div>-- <br><div class="m_966475192198582347m_-3807444873011230699gmail_signature"><div dir="ltr">Reju Sam John</div></div>
</font></span></div></div>
<br></div></div>______________________________<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.<wbr>org</a><br>
<br></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><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Reju Sam John</div></div>
</div>