%matplotlib inline import yt ds = yt.load('IsolatedGalaxy/galaxy0030/galaxy0030') slc = yt.SlicePlot(ds, 2, 'x-velocity', width=(20, 'kpc')) ax = slc.plots['x-velocity'].axes ax.patch.set_facecolor('black') slc