<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_analysis.api import HaloCatalog</font></div><div><font face="monospace, monospace" size="4" color="#ff00ff">data_ds = yt.load('/run/media/john/Seagate_Expansion_Drive/cosmo-sim_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_mass()</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><div><br></div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr">Reju Sam John</div></div>
</div></div>