#!/usr/bin/python

import glob

filelist=glob.glob('RD*_escape_fraction.txt')
filelist+=glob.glob('../../../../Normal_region1/Ionization/FOF/RD*_escape_fraction.txt')
filelist+=glob.glob('../../../../Rarepeak/Rarepeak_LWB/Ionization/FOF/RD*_escape_fraction.txt')

halo=[]

fo=open('largehalo1.txt','w')

for file in filelist:
    f=open(file,'r')
    lines = f.readlines()
    for line in lines[1:]:
        if float(line.split()[2])>= 8.0e8:
           halo.append(line)
           fo.write(line)

fo.close()