我需要一点帮助来创建一个派生场,它将返回引力势的负梯度,因此引力加速度。 我试过这样的事情:
@derived_field(name = "gradient_gravpot", units = " cm/s**2")
def _gradient (field, data):
G = add_gradient_fields(("gas"," gravitational_potental"))
return (- data["gravitational_potential_gradient_x"])
我收到了这个错误
NameError: global name 'add_gradient_fields' is not defined
请建议我定义该字段的正确方法。感谢
答案 0 :(得分:1)
看起来你得到了答案on the yt project mailing list。也就是说,我将在此重复add_gradient_fields
是Dataset
类的方法,因此您需要在加载的数据集上调用它。例如:
ds = yt.load('path/to/dataset')
ds.add_gradient_fields(('gas', 'gravitational_potential'))