#import h5py
import re
import glob
import numpy as np

nx=512 
ev_erg = 1.60217657e-12

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

position=np.zeros([3])

filelist = glob.glob("PopIII_haloes_xray-sources_*")

#f=open('PopIII_haloes_xray-sources_lifetime_1KeV.txt','r')

for name in filelist:
  redshift=float(re.findall(r"[-+]?\d*\.\d+|\d+",name)[0])
  f=open(name,'r')
  Xray = np.zeros([nx,nx,nx])

  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_Fortran.tofile("XrayLuminosity_%3.1f.dat"%redshift)
  Xray.tofile("XrayLuminosity_512_%3.1f.dat"%redshift)

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