import numpy as np
import commands
import re

tunit=6.88027e+14/3.156e7/1e6

for n, rs in zip([10,12,14,21,31,41],[20,19,18,17,16,15]):

 fi=open('Redshift%02i/rarepeak_galaxies_properties_luminosity_%04i'%(rs,n),'r')
 lines1 = fi.readlines()
 fi.close()

 fo=open('Redshift%02i/rarepeak_galaxies_properties_luminosity_noheader_%04i'%(rs,n),'w')

 region=lines1[0].split()[1].replace(',','')
 result=[float(x) for x in re.findall("[0-9.]+",lines1[0])]

 for line1 in lines1[2:]:
    line1=line1.rstrip('\r\n') 
    fo.write(line1+" %s %4.1f %12.8f\n"%(region,result[0],result[1]) )
    
 fo.close()