Issue #1326: Inconsistency between code_length / code_time and code_velocity in enzo_tiny_cosmology (yt_analysis/yt)
New issue 1326: Inconsistency between code_length / code_time and code_velocity in enzo_tiny_cosmology https://bitbucket.org/yt_analysis/yt/issues/1326/inconsistency-between-code_... John ZuHone: `code_velocity` and `code_length/code_time` should be the same units, but in `enzo_tiny_cosmology` this does not appear to be the case. Someone with more knowledge of Enzo should have a look at this--I tried to look but for cosmo sims the units are set in a way that is rather enzo-specific. ```python import yt ds = yt.load("enzo_tiny_cosmology/DD0046/DD0046") dd = ds.all_data() print(dd["velocity_divergence"].to("1/code_time")) print(dd["velocity_divergence"].to("code_velocity/code_length")) print(dd["velocity_magnitude"].to("code_velocity")) print(dd["velocity_magnitude"].to("code_length/code_time")) ``` ```pycon # velocity_divergence check YTArray([ 0.00186079, 0.00200255, 0.00191645, ..., -0.03508909, -0.04609586, -0.02762048]) 1/code_time YTArray([ 0.09485985, 0.10208633, 0.09769683, ..., -1.78877676, -2.34988201, -1.40804124]) code_velocity/code_length # velocity_magnitude check YTArray([ 0.01237452, 0.01096058, 0.00971304, ..., 0.02039839, 0.0223968 , 0.02642769]) code_velocity YTArray([ 0.00024274, 0.00021501, 0.00019053, ..., 0.00040014, 0.00043934, 0.00051841]) code_length/code_time ```
participants (1)
-
John ZuHone