<div>Hi Matt and Stephen,</div><div><br></div><div>I was able to fork YT (<span class="Apple-style-span" style="color: rgb(82, 82, 82); font-family: 'Bitstream Vera Sans Mono', 'DejaVu Sans Mono', Monaco, monospace; font-size: 12px; line-height: 18px; white-space: pre; ">hg clone <a href="https://gsiisg@bitbucket.org/gsiisg/el-yt">https://gsiisg@bitbucket.org/gsiisg/el-yt</a></span>) and implemented the</div>
<div><br></div><div>class AMREllipsoidBase(AMR3DData):</div><div><br></div><div>in data_containers.py, and right now I'm trying to put in the different methods of finding the ellipsoid parameters from the halo parameters.  I'm going to put in the geometry method first, and that requires the halo's</div>
<div><br></div><div>-center of mass x,y,z</div><div>-particle positions x,y,z</div><div>-grid edges left and right (for checking periodicity)</div><div><br></div><div>From Matt's previous email:</div><div>"I think how this would work the best would be to make two separate</div>
<div>methods on the Halo base class."</div><div><br></div><div>I was wondering where to put the addition to Halo base class, should it be in the halo_objects.py, or make a new file?</div><div><br></div><div>I'm still need to figure out (and this is where help would save me a ton of time), how to access the two files inside a class in halo_objects.py.  Both of which are created by</div>
<div><br></div><div>haloes.dump()</div><div><br></div><div>One of the file (.out) has the center of mass, and the other (.h5) has the particle positions.  In my script they were either in memory already when I ran the halo finder parallelHF, or loaded into memory by the single function LoadHaloes().  I did a h5ls on the .h5 files and don't see the center of mass information, so I assume it opened both the .out and .h5 and loaded them into memory.</div>
<div><br></div><div>And last but not least, I need to make another "dump" mechanism, or add to the existing one so it output the ellipsoid parameters.  Yes I did take Britton's advice and cut out the e2xyz, e3xyz and use a single tilt angle to replace them, which should cut out 5 columns of extra information.  The ellipsoid object will take *center_of_mass, A,B,C,*e1vector,tilt_angle as input, the ones with * means the input is a numpy array.</div>
<div><br></div><div>From</div><div>G.S.</div>