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

filelist = glob.glob("XrayLuminosity_512*")
filelist.sort(reverse=True)

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

  redshift=float(re.findall(r"[-+]?\d*\.\d+|\d+",name)[1])
  Density=np.fromfile('Density_%04i_%4.1f_normal.dat'%(n,redshift)).reshape(512,512,512)
  filtered = gaussian_filter(Density,6)

  np.save("Density_filtered_6_%4.2f"%redshift,filtered)

  n+=1