convert_yields.py.svn-base
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/usr/bin/env python

	import sys,string
	from numpy import *

	file1 = sys.argv[1]
	file2 = sys.argv[2]


	f = open(file1,'r')
	lines = f.readlines()
	f.close()


	data = array([])

	for line in lines:
	line = string.split(line)
	data = concatenate((data,array(line)))

	data.shape = (len(lines),len(data)/len(lines))

	data[0,0] = '# M'

	f = open(file2,'w')

	# header
	f.write("%s\t"%data[0,0])

	for j in xrange(1,data.shape[0]):
	f.write("%s\t\t"%data[j,0])

	f.write("\n")


	# 10 M
	f.write("%s\t"%"10")
	for j in xrange(1,data.shape[0]):
	f.write("%s\t"%"0.00E-00")

	f.write("\n")



	for i in xrange(1,data.shape[1]):
	f.write("%s\t"%data[0,i])

	for j in xrange(1,data.shape[0]):
	f.write("%s\t"%data[j,i])

	f.write("\n")


	f.close()