compute_Yields.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/usr/bin/env python


	import libcris
	import sys,os
	from numpy import *


	Ms = libcris.Mcode.keys()
	Zs = libcris.Zcode.keys()

	Ms.sort()
	Zs.sort()


	# first, get elts

	filename = "Net/yields_net_%s_20130626_134213.txt"%(libcris.Mcode[Ms[0]] + libcris.Zcode[Zs[0]])
	Elts,Masses = libcris.ReadTable(filename)




	for Z in Zs:

	outputname = "cristallo-Yields-%6.4f.txt"%Z

	# open file and write header
	f = open(outputname,'w')
	f.write("# M")
	for elt in Elts:
	f.write(" %14s"%elt)
	f.write("\n")



	for M in Ms:

	if M==3.0 and Z<=0.0003: # avoid missing files
	Z=0.0010

	filename = "Net/yields_net_%s_20130626_134213.txt"%(libcris.Mcode[M] + libcris.Zcode[Z])
	Elts,Masses = libcris.ReadTable(filename)

	# write yields
	f.write("%5.2f"%M)
	for i,elt in enumerate(Elts):
	f.write(" %14.7e"%Masses[i])
	f.write("\n")

	f.close()