import numpy as np
import glob
import re
from scipy.ndimage.filters import gaussian_filter

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)

  filtered = gaussian_filter(Density,6)
  np.save("Density_filtered_6_%4.2f"%redshift,filtered)

  n+=1