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

# In[1]:


import yt

ds = yt.load('IsolatedGalaxy/galaxy0030/galaxy0030')


# In[2]:


slc = yt.SlicePlot(ds, 2, 'density')

slc


# In[3]:


# Unfortunately the default labels aren't great...
slc.set_axes_unit('1000*kpc')


# In[4]:


slc.set_xlabel(r'$\rm{x}\ \ (10^3\ \rm{kpc})$')
slc.set_ylabel(r'$\rm{y}\ \ (10^3\ \rm{kpc})$')


# In[ ]: