from yt.mods import *
from matplotlib import use; use('Agg')
import matplotlib.pyplot as plt
import h5py

pf = load("output_0171")
#pf = load("RedshiftOutput0002")

field_list = ['Density','Dark_Matter_Density','Temperature']
#field_list = ['Dark_Matter_Density']
nx = 512
c = [0.5, 0.5, 0.5]
cl = [0.0,0.0,0.0]
cr = [1.0,1.0,1.0] 

#region = pf.h.smoothed_covering_grid(12,left_edge=cl,right_edge=cr,dims=[nx,nx,nx])
region = pf.h.covering_grid(0,left_edge=cl,right_edge=cr,dims=[nx,nx,nx])

for field in field_list:
    region.get_data('%s'%field)
    #data_Fortran = na.asfortranarray(region['%s'%field])
    region['%s'%field].tofile("%s_512.dat"%field)


#f=h5py.File('XrayFlux.h5','w')
#dataset=f.create_dataset("XRayFlux",data=region['XRayFlux'])
#dataset=f.create_dataset("deltaTb",data=proj1['deltaTb']/l)
#f.close()

#region.write_out('XrayFlux.h5',fields=['XrayFlux'])