seems to do nothing.
>>> 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.
--