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)
   a[206:306,206:306,206:306] = 0.0
   print a.sum()/(512**3-100**3)