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