# Read a set star lists and sort the by exposre time
# 

from FITS import *
import re
import os
from string import *
import fileinput

imagenames = '^H.*fts'
imagere = re.compile(imagenames)
outfiles = ['gsc15s-1.dat', 'gsc150s-1.dat', 'gsc15s-2.dat', 'gsc150s-2.dat', 'gsc15s-3.dat', 'gsc150s-3.dat', 'gsc15s-4.dat', 'gsc150s-4.dat', 'gsc15s-5.dat', 'gsc150s-5.dat']
times = [15, 150]
fps = []
for a in outfiles:
    fps.append(open(a, 'w'))
files = os.listdir('.')
for a in files:
    if imagere.search(a) != None:
    	if os.path.isfile(a):
	    #  Read keywords
	    keywords = getFITSkeywords(a)
	    exposure = float(keywords['EXPOSURE'][0])
	    print a, exposure
	    for i in range(len(times)):
		if int(exposure) == times[i]:
		    for j in range(1, 6):
		    	#  Create file name
		        datname = split(a, '.')[0] + ('-%d.dat' % (j))
		    	for line in fileinput.input(datname):
			    #  Skip the first line
			    if fileinput.lineno() == 1: continue
			    #  Copy line to output, if it has match to ACT star
			    try:
			        if split(line)[6] == 'g': fps[i + (j - 1) * 2].write(line)
			    except:
			        pass