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

output = 29

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)
svm=sv['CellMass'].sum()
wv = pf.h.region([0.5,0.5,0.5],(.0,.0,.0),(1.,1.,1.))
wvm = wv['CellMass'].sum()

f=open('massratio.txt','w')
f.write(str(svm)+" "+str(wvm))
#halo_list=parallelHF(pf,dm_only=True,subvolume = sv)
#center=np.array(halo_list[0].center_of_mass())
#halo_list.write_out("sv.out")
#halo_list.dump("RD%04i"%output)