import numpy as np

n=93

f1=open('RedshiftOutput%04i-halos.dat'%n,'r')
f2=open('escape_fraction_RD%02i.txt'%n,'r')

fo=open('RD%04i_escape_fraction.txt'%n,'w')
fo.write("%8s %12s %12s %12s %12s %12s %12s %16s\n"%\
   ('halo #','M_{tot}','M_{vir}','M_{gas}','M_{star}','M_{astar}','r_{vir}','escape fraction'))


lines1=f1.readlines()
lines2=f2.readlines()

f1.close()
f2.close()

i=1
j=0
for i in xrange(1,len(lines1)):
    if float(lines1[i].split()[0]) < 0.0: continue
    if not lines2[j].strip(): 
       j=j+1
       continue
    splited1=lines1[i].split()
    splited2=lines2[j].split()
    j=j+1
    fo.write("%8s %12s %12s %12s %12s %12s %12s %16s\n"%\
         (splited2[1].replace(':',''),splited1[5],splited1[4],splited1[7],splited1[6],\
          splited2[7].replace(',',''),splited1[8],splited2[10]))
    if j >= len(lines2): break

fo.close()