#!/usr/bin/env python from scipy import * a=zeros((132, 100), dtype=float) #a=loadtxt('spin.dat') #print a ll=0 f=open('spin.dat', 'r') print '#' a=[] #f.readline() for l in f: a=[float(val) for val in l.split()] b=int(a[4]) for i in range(b): print ll, a[i+b+15] ll+=1 f.close()