Source

yt.geometry_tests / test_ramses_io.py

Full commit
from yt.mods import *
import time

pf = load("output_00045/info_00045.txt")
pf.h
t1 = time.time()
for dom in pf.h.domains:
    print dom.domain_id, len(dom.hydro_offset)
t2 = time.time()
print "%0.3e" % (t2-t1)

dd = pf.h.all_data()
t1 = time.time()
rho = dd["Density"]
t2 = time.time()
print "%0.3e" % (t2-t1)
print float((rho == 0).sum())/rho.size