gistfile1.txt

In [1]: import yt

In [2]: ds = yt.load('ArepoBullet/snapshot_150.hdf5', smoothing_factor=1.5)
yt : [INFO     ] 2020-02-27 15:45:18,078 Omega Lambda is 0.0, so we are turning off Cosmology.
yt : [INFO     ] 2020-02-27 15:45:18,138 Parameters: current_time
= 1.5380632172269546
yt : [INFO     ] 2020-02-27 15:45:18,138 Parameters: domain_dimensions
= [1 1 1]
yt : [INFO     ] 2020-02-27 15:45:18,139 Parameters: domain_left_edge
= [0. 0. 0.]
yt : [INFO     ] 2020-02-27 15:45:18,139 Parameters: domain_right_edge
= [40000. 40000. 40000.]
yt : [INFO     ] 2020-02-27 15:45:18,139 Parameters: cosmological_simulation
= 0

In [3]: ds.index
yt : [INFO     ] 2020-02-27 15:45:21,668 Allocating for 2.653e+07 particles
Initializing coarse index : 100%|██████████████| 51/51 [00:02<00:00, 18.57it/s]
Initializing refined index: 100%|██████████████| 51/51 [00:11<00:00,  4.53it/s]
Out[3]: <yt.frontends.sph.data_structures.SPHParticleIndex at 0x7f152f726390>

In [4]: ds.find_max(('gas', 'density'))
yt : [INFO     ] 2020-02-27 15:45:45,373 Max Value is 3.96176e-26 at 18600.6230468750000000 19997.5761718750000000 20007.3691406250000000
Out[4]:
(unyt_quantity(3.96175712e-26, 'g/cm**3'),
 unyt_array([18600.62304688, 19997.57617188, 20007.36914062], 'code_length'))