convert_yields.py
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] = '# %s'%data[0,0]

	f = open(file2,'w')

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

	for j in xrange(1,data.shape[0]):
	f.write("%s\t\t"%data[j,0])
	f.write("%s\t\t"%'Metals') # metals
	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])
	s = sum(data[:,i].astype(float))
	f.write("%8.2E\t"%s) # metals
	f.write("\n")


	f.close()