import glob
import os
import numpy as np
import re
import math

tc = 197.5/772
nx=512
Xray_r=np.zeros([nx**3])
filelist = glob.glob("Temperature_512_*.dat")
for i,name in enumerate(filelist):
  redshift = float(re.findall(r"[-+]?\d*\.\d+|\d+",name)[1])
  #if redshift >= 12.50: continue
  #if redshift <= 7.5: continue
  f = open(name,'r')
  #f.seek(12,os.SEEK_SET)
  tem=(np.fromfile(f,dtype='float64'))
  tem[tem<1000.0*((1+redshift)/100)**2] = tem[tem<1000.0*((1+redshift)/100)**2]*tc
  print "Reading %s"%name
  f.close()
  os.system('mv %s %s_old'%(name,name))
  tem.tofile(name)
  del tem
  print "Writing %s"%name