我试图用屏蔽场屏蔽系统的边界密度
@derived_field(name = "bound_density", units = "g/cm**3")
def _get_ejected_density(field, data):
E = 0.5* data["cell_mass"]* (data["velx"]**2+data["vely"]**2+data["velz"]**2)+ data["gpot"]*data["cell_mass"]
return ((np.array(E)<0)*1)*data["density"]
它还返回带有白色条的未绑定密度以及绑定密度,但我不想在这里显示未绑定的内容。
我想知道是否有一种方法可以只显示该图中的绑定材料。另一种解决方案是将未绑定密度颜色条设置为与地板条匹配,因此为 1.0e0,这样虽然绘图仍显示未绑定密度,但它与绘图中的最低密度颜色条完全匹配,因此无法区分。