在項目中創建場的梯度作爲派生字段yt

Question

我需要一些幫助來創建派生字段，該字段將返回引力勢的負梯度，即重力加速度。 我想是這樣的：

@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 

請建議我去定義字段正確的方法。謝謝

Answer 1

它看起來像你有一個答案on the yt project mailing list。這就是說，我在這裏重複說明add_gradient_fields是Dataset類的一種方法，因此您需要在加載的數據集上調用它。例如：

ds = yt.load('path/to/dataset') 
ds.add_gradient_fields(('gas', 'gravitational_potential'))