import yt
ds = yt.load("GaussianCloud/data.0077.3d.hdf5")
ad = ds.all_data()
plot = yt.PhasePlot(ad, "density", "magnetic_energy", ["ones"],
                    fractional=True, weight_field=None)
plot.save()