gointersects.py

   1 #
   2 #  gointersects.py
   3 #  ENRAGE
   4 #
   5
   6 import sys
   7
   8 print "gointersects: version 1.1"
   9
  10 def main(argv=None):
  11     if not argv:
  12         argv = sys.argv
  13
  14     if len(argv) < 4:
  15         print "usage: python %s gogidfile gidfile outfile" % argv[0]
  16         sys.exit(1)
  17
  18     gogidfilename = argv[1]
  19     gidfilename = argv[2]
  20     outfilename = argv[3]
  21
  22     gointersects(gogidfilename, gidfilename, outfilename)
  23
  24
  25 def gointersects(gogidfilename, gidfilename, outfilename):
  26     gidList = []
  27     gogidfile = open(gogidfilename)
  28     for line in gogidfile:
  29         fields = line.split()
  30         gidList.append(fields[0])
  31
  32     gogidfile.close()
  33
  34     gidfile = open(gidfilename)
  35     outfile = open(outfilename, "w")
  36     for line in gidfile:
  37         fields = line.split()
  38         if fields[0] in gidList:
  39             outfile.write(line)
  40
  41     gidfile.close()
  42     outfile.close()
  43
  44
  45 if __name__ == "__main__":
  46     main(sys.argv)