import yt from yt.units import kpc ds = yt.load('IsolatedGalaxy/galaxy0030/galaxy0030') my_galaxy = ds.disk(ds.domain_center, [0.0, 0.0, 1.0], 10*kpc, 3*kpc) p = yt.ParticlePlot(my_galaxy, 'particle_position_x', 'particle_position_y') p.set_ylog(False) p.set_xlog(False) p.set_line_property('markersize', 2) p.set_unit('particle_position_x', 'kpc') p.set_xlim(485, 515) p.set_unit('particle_position_y', 'kpc') p.set_ylim(485, 515) p.show()