Colourdiscon.py

From CUC3
Revision as of 19:02, 8 May 2009 by import>Mp466
Jump to navigation Jump to search

Insert non-formatted text here #!/usr/bin/python import numpy import sys import getopt """ print files (test.out.int.*) comprised of a set of minima sorted into a set of bins by an order parameter Used for colouring discon graphs Order Parameter Minima # file format AMHQ 0.1653832 6 """

x=[] y=[] bins = [0.0, 0.4, 0.6, 0.8, 1.0]

filename = sys.argv[1]

f=open(filename,"r")

for line in f:

       tokens=line.split()
       x.append(float(tokens[1]))
       y.append(int(tokens[2]))

total=[]

for a, b in enumerate(bins):

   total.append([])
   for i in range(len(x)):
      if x[i] >= b and x[i] < bins[a+1] :
              total[a].append(y[i])

for a, b in enumerate(bins):

     filefile="test.out.int"+str(a)
     numpy.savetxt(filefile, total[a], fmt='%d')