import yt

ds = yt.load('maestro_subCh_plt00248')

for normal in [2, [.01, .01, 1]]:
    yt.SlicePlot(ds, normal, 'density', width=(3e4, 'km')).show()