import yt
ds = yt.load('IsolatedGalaxy/galaxy0030/galaxy0030')
slc = yt.SlicePlot(ds, 2, 'density')
slc
# Unfortunately the default labels aren't great...
slc.set_axes_unit('1000*kpc')
slc.set_xlabel(r'$\rm{x}\ \ (10^3\ \rm{kpc})$')
slc.set_ylabel(r'$\rm{y}\ \ (10^3\ \rm{kpc})$')