import yt

ds = yt.load('Enzo_64/DD0043/data0043')

ds.add_deposited_particle_field(("io", "particle_position_cylindrical_radius"), 'cic')
ds.add_deposited_particle_field(("io", "particle_position_cylindrical_theta"), 'cic')
ds.add_deposited_particle_field(("io", "particle_position_cylindrical_z"), 'cic')
ds.add_deposited_particle_field(("io", "particle_position_spherical_radius"), 'cic')
ds.add_deposited_particle_field(("io", "particle_position_spherical_theta"), 'cic')
ds.add_deposited_particle_field(("io", "particle_position_spherical_phi"), 'cic')

test_fields = []
for field in ds.derived_field_list:
    if field[0] == 'deposit' and ('cylindrical' in field[1] or "spherical" in field[1]) and 'io_' in field[1]:
        test_fields.append(field)

from pprint import pprint
pprint(test_fields)

for field in test_fields:
    print field
    for d in 'xyz':
        plot = yt.SlicePlot(ds, d, field)
        plot.set_figure_size(5)
        plot.show()