import glob
import re
import numpy as np

filelist = glob.glob("Final_Temperature_*.dat")
filelist.sort()
filelist.reverse()
for name in filelist:
 z = float(re.findall(r"[-+]?\d*\.\d+|\d+",name)[1])
 if z < 25.0 :
   a=np.fromfile(name).reshape(512,512,512)
   b=np.fromfile('../Density_512_%05.2f.dat'%z).reshape(512,512,512)
   #a[206:306,206:306,206:306] = 0.0
   b[206:306,206:306,206:306] = 0.0
   print (a*b).sum()/b.sum()