Source

yt.geometry_tests / spatial_fields.py

Full commit
1
2
3
4
5
from yt.mods import *

pf = load("DD0039/output_0039")
sp = pf.h.sphere("c", (10.0, 'pc'))
print "STATS", sp["DivV"].max(), sp["DivV"].size, sp["DivV"].sum(dtype="float64")