Hello,<div><br></div><div>I was trying to get the position and mass information for the dark matter (and stellar) particles in an Enzo simulation.  I have tried a couple of different ways and I was curious about the units on the following method. In this case I am presumably accessing the dark matter particles (particle_type = 1)</div>
<div><br></div><div><i>from yt.mods import *</i></div><div><i>pf = load("DD0252/DD0252")</i></div><div><i>array = pf.h.find_particles_by_type(ptype=1, max_num=3)</i></div><div><i>print "(%.5f, %.5f, %.5f) = %f\n" %(array['particle_position_x'][0], array['particle_position_y'][0], array['particle_position_z'][0], array['particle_mass'][0])</i></div>
<div><div><br></div>This returns <b>(0.01055, 0.01873, 0.00224) = 0.835448. </b>What are the mass units on the <b>0.835448</b>?  I don't think that they are cgs, and there doesn't appear to be a conversion factor available via pf[' * '].  If I continue and enter </div>
<div><br></div><div><i>print pf.field_info['Density'].get_units()</i></div><div><div><i>print pf.field_info['particle_mass'].get_units()</i></div></div><div><i><br></i></div><div>The 'Density' line returns <b>\rm{g}/\rm{cm}^3, </b>but the 'particle_mass' line returns nothing. So I can't get the units from here. I can, however, access the same particle (or at least a particle at the same location) by using</div>
<div><br></div><div><i>from yt.mods import *</i></div><div><i>pf = load("DD0252/DD0252")</i></div><div><i>dd = pf.h.all_data()</i></div><div><i>print "(%f %f %f) =  %f Msun, ptype= %i\n"%(dd["particle_position_x"][0], dd["particle_pos    ition_y"][0], dd["particle_position_z"][0], dd["ParticleMassMsun"][0],dd["particle_type"][0])</i></div>
<div><br></div><div>Which returns <b>(0.010554 0.018727 0.002239) =  82255596.670652 Msun, ptype= 1</b></div><div><br></div><div>This seems to give reasonable values. Both methods should be equivalent ways of to access the particle (both stellar and dark matter) data, correct?  So how do I get my units in solar masses when using <i>pf.h.find_particles_by_type()</i>?  Thanks for your time!</div>
<div><br></div><div>~ali</div>