import glob
import re
import numpy as np
import os

os.chdir('../300eV')
filelist = glob.glob("Delta_Electron_fraction_*.dat")
os.chdir('../synthetic')
filelist.sort()
for name in filelist:
 z0 = float(re.findall(r"[-+]?\d*\.\d+|\d+",name)[0])
 z1 = float(re.findall(r"[-+]?\d*\.\d+|\d+",name)[1])
 if z1 < 6: continue
 if z1 > 6: continue
 oldtemperature = np.fromfile("../770eV/HII_Fraction_512_%5.2f.dat"%z0)
 for folder,weight in zip(['../300eV','../500eV','../770eV','../1keV','../3kev'],[0.134,0.125,0.089,0.219,0.433]): #[ 0.288,0.307,0.200,0.154,0.051]):
   deltatem = np.fromfile(folder+"/Delta_Electron_fraction_%5.2f_%05.2f.dat"%(z0,z1))
   oldtemperature += deltatem*weight
 oldtemperature.tofile('Final_Electron_fraction_%5.2f_%05.2f.dat'%(z0,z1))