import sys
from yt.mods import *
import numpy as np
#from yt.analysis_modules.halo_finding.api import *

output = 38

#fn = 'output_%04i'%(output)
fn = 'RedshiftOutput%04i'%(output)

pf = load(fn)
vl=pf["RefineRegionLeftEdge"]
vr=pf["RefineRegionRightEdge"]
sv = pf.h.region((vl+vr)/2, vl, vr)
average_value = sv.quantities["WeightedAverageQuantity"]('H2I_kdiss', 'CellVolume')
print average_value