#import h5py
import re
import glob
import numpy as np

nx=512
ev_erg = 1.60217657e-12
redshift = 17.5

Xray = np.zeros([nx,nx,nx])
dx = 1.0/nx

position=np.zeros([3])

f=open('PopIII_haloes_xray-sources_17.45.txt','r')

line=f.readline()
while line !='':
   line = f.readline()
   position[0]=float(line.split()[2])
   position[1]=float(line.split()[3])
   position[2]=float(line.split()[4])
   line = f.readline()
   luminosity = float(line.split()[2])
   for i in xrange(4):
     line = f.readline()
   photonenergy = float(line.split()[2])
   luminosity = luminosity*photonenergy*ev_erg
   Xray[int(position[0]/dx),int(position[1]/dx),int(position[2]/dx)] = \
     Xray[int(position[0]/dx),int(position[1]/dx),int(position[2]/dx)]+luminosity
   line=f.readline();line=f.readline()

f.close()
#Xray_Fortran = np.asfortranarray(Xray)
Xray.tofile("XrayLuminosity_%i_%3.1f.dat"%(nx,redshift))

#f=h5py.File('XrayLuminosity.h5','w')
#dataset=f.create_dataset("XRayLuminosity",data=Xray)
#f.close()