import yt
ds = yt.load('IsolatedGalaxy/galaxy0030/galaxy0030')

p = yt.ParticlePlot(ds, 'particle_position_x', 'particle_position_y')
p.show()

p.set_width(500, 'kpc')
p.show()

p = yt.ParticlePlot(ds, 'particle_position_x', 'particle_position_y', color='g', width=(0.1, 0.1))
p.show()

p = yt.ParticlePlot(ds, 'particle_position_x', 'particle_position_y', 'particle_mass', width=(0.5, 0.5))
p.set_unit('particle_mass', 'Msun')
p.show()

p.zoom(32)
p.annotate_title('Zoomed-in Particle Plot')

p = yt.ParticlePlot(ds, 'particle_position_x', 'particle_velocity_z')
p.show()

p.set_unit('particle_position_x', 'Mpc')
p.set_unit('particle_velocity_z', 'km/s')
p.show()

p = yt.ParticlePlot(ds, 'particle_position_x', 'particle_velocity_z', ['particle_mass'])
p.set_unit('particle_position_x', 'Mpc')
p.set_unit('particle_velocity_z', 'km/s')
p.set_unit('particle_mass', 'Msun')
p.show()

p = yt.ParticlePlot(ds, 'particle_velocity_x', 'particle_velocity_y', 'particle_mass')
p.set_unit('particle_velocity_x', 'km/s')
p.set_unit('particle_velocity_y', 'km/s')
p.set_unit('particle_mass', 'Msun')
p.set_ylim(-400, 400)
p.set_xlim(-400, 400)
p.show()

p = yt.ParticlePlot(ds, 'particle_velocity_x', 'particle_velocity_y', 'particle_mass', weight_field='particle_mass')
p.set_unit('particle_velocity_x', 'km/s')
p.set_unit('particle_velocity_y', 'km/s')
p.set_unit('particle_mass', 'Msun')
p.set_ylim(-400, 400)
p.set_xlim(-400, 400)
p.show()

p = yt.ParticlePlot(ds, 'particle_velocity_x', 'particle_velocity_y', 'particle_ones')
p.set_unit('particle_velocity_x', 'km/s')
p.set_unit('particle_velocity_y', 'km/s')
p.set_ylim(-400, 400)
p.set_xlim(-400, 400)
p.show()

p = yt.ParticlePlot(ds, 'particle_position_x', 'particle_position_y', 'particle_mass', 
                    weight_field='particle_mass', width=(0.01, 0.01))
p.show()