import numpy as np
import glob
import re

filelist1 = glob.glob("XrayLuminosity_512_??.??.dat")
#filelist2 = glob.glob("XrayLuminosity_WV_??.??.dat")
filelist1.sort(reverse=True)
#filelist2.sort(reverse=True)
filelist3 = glob.glob("XrayLuminosity_512_?.??.dat")
#filelist4 = glob.glob("XrayLuminosity_WV_?.??.dat")
filelist3.sort(reverse=True)

filelist1+=filelist3

n=0
outarray=np.empty(shape=[0,3])
for name in filelist1:

  redshift=float(re.findall(r"[-+]?\d*\.\d+|\d+",name)[1])
  Xray=np.fromfile(name).reshape(512,512,512)
  if n<15:
      Density=np.fromfile('Density_%04i_%4.1f_normal.dat'%(n,redshift)).reshape(512,512,512)
  else: 
      Density=np.fromfile('Density_%04i_%4.2f_normal.dat'%(n,redshift)).reshape(512,512,512)

  #index = np.argwhere(Xray > 0)
  #xmin = index[:][:,0].min()
  #xmax = index[:][:,0].max()
  #ymin = index[:][:,1].min()
  #ymax = index[:][:,1].max()
  #zmin = index[:][:,2].min()
  #zmax = index[:][:,2].max()

  Density = Density[Xray>0]
  X=Xray[Xray>0]
  z=np.ones_like(X)*redshift
  temp = np.array([z,Density,X])
  outarray=np.append(outarray,temp.T,axis=0)
  print temp.T
  #print outarray
  #Xray=np.zeros([512,512,512])
  #Xray[Density>Density_mean] = Xmean

  #np.savetxt('Density_Xray_%4.2f.dat'%redshift,outarray.T)
  n+=1

np.savetxt('Density_Xray.dat',outarray)