def dpdy(field, data): output = xo.grad(data, 'pressure', 1) return output
def dpdz(field,data): output = xo.grad(data,'gas_pressure',2) return output
def dpdy(field, data): output = xo.grad(data, YT_gas_pressure, 1) return output