Colourdiscon.py

From CUC3
Revision as of 19:55, 8 May 2009 by import>Mp466
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search
  1. !/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')