import glob
import re
import numpy as np
import os

redshift0=15.0
os.chdir('../300eV')
filelist = glob.glob("Delta_Temperature_*.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 < 9: continue
 if z1 >= 9: continue
 oldtemperature = np.fromfile("../Temperature_512_%05.2f.dat"%redshift0)
 oldtemperature = oldtemperature*((1.0+z1)/(1.0+redshift0))**2
 oldtemperature0 = oldtemperature
 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+"/Final_Temperature_%05.2f_%05.2f.dat"%(z0,z1))-oldtemperature0
   oldtemperature = oldtemperature + deltatem*weight
 oldtemperature.tofile('Final_Temperature_%05.2f_%05.2f.dat'%(z0,z1))