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)
[('deposit', 'io_cic_position_cylindrical_radius'), ('deposit', 'io_cic_position_cylindrical_theta'), ('deposit', 'io_cic_position_cylindrical_z'), ('deposit', 'io_cic_position_spherical_phi'), ('deposit', 'io_cic_position_spherical_radius'), ('deposit', 'io_cic_position_spherical_theta')]
for field in test_fields:
print field
for d in 'xyz':
plot = yt.SlicePlot(ds, d, field)
plot.set_figure_size(5)
plot.show()
('deposit', 'io_cic_position_cylindrical_radius')
/usr/local/lib/python2.7/site-packages/matplotlib/axes/_base.py:1039: UnicodeWarning: Unicode equal comparison failed to convert both arguments to Unicode - interpreting them as being unequal if aspect == 'normal': /usr/local/lib/python2.7/site-packages/matplotlib/axes/_base.py:1044: UnicodeWarning: Unicode equal comparison failed to convert both arguments to Unicode - interpreting them as being unequal elif aspect in ('equal', 'auto'):
('deposit', 'io_cic_position_cylindrical_theta')
('deposit', 'io_cic_position_cylindrical_z')
('deposit', 'io_cic_position_spherical_phi')
('deposit', 'io_cic_position_spherical_radius')
('deposit', 'io_cic_position_spherical_theta')