1import yt 2import yt.units as u 3 4ds = yt.load("HiresIsolatedGalaxy/DD0044/DD0044") 5 6center = [0.53, 0.53, 0.53] 7normal = [0, 0, 1] 8radius = 40 * u.kpc 9height = 5 * u.kpc 10 11disk = ds.disk(center, [0, 0, 1], radius, height) 12 13profile = yt.create_profile( 14 data_source=disk, 15 bin_fields=[("index", "radius")], 16 fields=[("gas", "velocity_cylindrical_theta")], 17 n_bins=256, 18 units=dict(radius="kpc", velocity_cylindrical_theta="km/s"), 19 logs=dict(radius=False), 20 weight_field=("gas", "mass"), 21 extrema=dict(radius=(0, 40)), 22) 23 24plot = yt.ProfilePlot.from_profiles(profile) 25 26plot.set_log(("gas", "velocity_cylindrical_theta"), False) 27plot.set_ylim(("gas", "velocity_cylindrical_theta"), 60, 160) 28 29plot.save() 30