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