Hi Andrew,<br><br>I'm not sure about the output from get_source, so I'll let someone else comment on that.  In general, the units of the fields are in cgs.  The exceptions to this are positions, both grid cell and particle positions, and any field that has units implied by the name, such as CellMassMsun.  There is a units conversion dict, pf.units, that you can use to convert positions into usable units.  For example, multiply position by pf.units['mpc'] to get proper Mpc.<br>
<br>Britton<br><br><div class="gmail_quote">On Wed, Aug 25, 2010 at 9:47 AM, Andrew Davis <span dir="ltr"><<a href="mailto:andrew.davis@yale.edu">andrew.davis@yale.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Hi yt users,<br>
<br>
I have some basic questions about the units.  I'm trying to follow through the<br>
cookbook and help site, but some of it isn't working.  Namely, the part about<br>
having it tell me what units some of the fields are in.<br>
<br>
When I follow the commands on the page:<br>
<a href="http://yt.enzotools.org/doc/howto/objects.html" target="_blank">http://yt.enzotools.org/doc/howto/objects.html</a>  This is what happens:<br>
<br>
In: pf=load("EZ2_0000")<br>
In: print(pf.h.field_list)<br>
Out: ['z-velocity', 'Temperature', 'Density', 'particle_velocity_z',<br>
'TotalEnergy', 'Dark_Matter_Density', 'particle_position_z', 'x-velocity',<br>
'particle_position_y', 'y-velocity', 'particle_index', 'particle_position_x',<br>
'particle_mass', 'particle_velocity_x', 'particle_velocity_y']<br>
<br>
So, I see DM particle fields, which is what I'm interested in.  However, when I<br>
try to find their units, I get:<br>
In: print(pf.field_info["particle_velocity_x"].units)<br>
Out: ---------------------------------------------------------------------------<br>
AttributeError                            Traceback (most recent call last)<br>
<br>
/home/adavis/Research/YT/src/yt-1.7-svn/scripts/iyt in <module>()<br>
----> 1<br>
      2<br>
      3<br>
      4<br>
      5<br>
<br>
AttributeError: 'DerivedField' object has no attribute 'units'<br>
<br>
I get the same error for the following fields:<br>
[xyz]-velocity, particle_velocity_[xyz], particle_position_[xyz], Density,<br>
TotalEnergy, Dark_Matter_Density, and particle_index (ok, this shouldn't have<br>
units)<br>
<br>
When I try the other method of getting at the units:<br>
pf.field_info["particle_velocity_x"].get_source() , there is nothing there that<br>
says anything about units.<br>
<br>
Out: def _Particles(field, data):<br>
        io = <a href="http://data.hierarchy.io" target="_blank">data.hierarchy.io</a><br>
        if not data.NumberOfParticles > 0:<br>
            return na.array([], dtype=dtype)<br>
        try:<br>
            return io._read_data_set(data, p_field).astype(dtype)<br>
        except io._read_exception:<br>
            pass<br>
        # This is bad.  But it's the best idea I have right now.<br>
        return data._read_data(p_field.replace("_"," ")).astype(dtype)<br>
<br>
And, while field_info["Density"].get_source() works, the same command for<br>
Temperature tells me that the units are g/cm^3, which is obviously wrong.<br>
<br>
So, my major question: particle positions seem to be in the [0-1] unit system,<br>
so conversion is simple for them.  Particle velocities, however, are unclear.<br>
They look like they're in cm/s, but I'd like to be sure.<br>
<br>
thanks!<br>
<br>
Andrew<br>
<br>
<br>
<br>
*************************************************************<br>
**  Andrew J. Davis             <a href="mailto:andrew.davis@yale.edu">andrew.davis@yale.edu</a>      **<br>
**  Dept. of Astronomy          203-432-5119               **<br>
**  Yale University             <a href="http://www.astro.yale.edu/adavis" target="_blank">www.astro.yale.edu/adavis</a>  **<br>
*************************************************************<br>
<br>
<br>
<br>
_______________________________________________<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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</blockquote></div><br>