Hi, Matt-
 
It's neither a bug nor a feature; "density" is represented as a
"code_mass / code_length**3" primitive.

So what does code_density do?

Thanks a bunch!
d.
 

-Matt

On Fri, May 22, 2015 at 9:20 AM, David Collins <dcollins4096@gmail.com> wrote:
> Hi, Everyone--
>
> I'm trying to get code units from enzo for a cosmology run, but
> in_units('code_density')
> seems to do nothing.
> 'code_mass/code_length**3'
> seems to work.  Is this a bug, or a feature?
>
>
>>>> ad['density'].in_units('code_mass/code_length**3')
> YTArray([ 0.9970247 ,  0.99774731,  0.99854185, ...,  0.99993769,
>         0.99914176,  0.99841755]) code_mass/code_length**3
>>>> ad['density'].in_units('code_density')
> YTArray([  1.87321000e-29,   1.87456764e-29,   1.87606042e-29, ...,
>          1.87868293e-29,   1.87718755e-29,   1.87582689e-29]) code_density
>
> Thanks!
> d.
>
> --
> -- Sent from a computer.
>
> _______________________________________________
> yt-users mailing list
> yt-users@lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org



--
-- Sent from a computer.