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!

-- Sent from a computer.