我正在尝试使用 yt-project 函数“yt.PhasePlot”制作相图（参见此处，在“2D Phase Plots”段落中）。我可以在 y 轴上创建压力图，在 x 轴上创建密度图，这样做：
ds = yt.load('file')
ad = ds.all_data()
plot = yt.PhasePlot(ad, ("gas", "density"), ("gas", "pressure"), [("gas", "cell_mass")], weight_field = None, fractional=(True))
``
类似于提供的链接中显示的示例。我需要做的是通过玻尔兹曼常数 kB 来缩放 y 轴上的压力，这是一个浮点数。由于 ("gas", "pressure") 是我所做的元组：
tuple(ad["gas", "pressure"]/kB)
但这需要很多时间，当我在 PhasePlot 函数中提供这段代码时，它不起作用。我还尝试将压力的 cgs 单位设置为 K/cm^3，以强制将压力除以 kB，但效果不佳。还有其他方法吗？