#!/usr/bin/env python
# coding: utf-8

# In[1]:


import yt
ds = yt.load('enzo_cosmology_plus/DD0046/DD0046')
p = yt.ProjectionPlot(ds, 2, 'density')


# In[2]:


density_plot = p.plots['density']


# In[3]:


density_plot.image.get_clim()


# In[4]:


density_plot


# In[ ]: