I am using arbitrary_grid to access my data without flattening:
> dd=ds.arbitrary_grid(left_edge=[(com_x_1.value-30), (com_y_1.value-30),
> (com_z_1.value-30)]*kpc, right_edge=[(com_x_1.value+30), (com_y_1.value+30),
> (com_z_1.value+30)]*kpc, dims=[128, 128, 128])
and as I attempt to get the gas fields for this object,
> print(dd['gas', 'velocity_x'])
but end up getting the error: (which isn't the case for deposit fields)
yt.units.unit_registry.UnitParseError: Code units have not been defined.
> Try creating the array or quantity using ds.arr or ds.quan instead.
If I'm thinking right, is it because the arbitrary_grid is only applied to
the particles. In that case, is there any equivalent for gas fields where I
can access my 3D data for the gas as it is , without flattening? Thanks