Files F93056894

compute_Yields.py
No OneTemporary
Actions

Subscribers

None

File Metadata

Created: Mon, Nov 25, 21:12

compute_Yields.py
View Options

	#!/usr/bin/env python

	from numpy import *

	costSr = 3.16 * 88./138.
	costLa = 0.136* 1.
	costZr = 2.53 * 90./138.
	costEu = 0.117*151./138./2.
	costY = 1.625* 89./138./1.6

	Ms = array([7,8,10,12,15,30],float)





	#Ba = array([0.0,6.4e-06,3.0e-08,3.6e-08,4.5e-08,9.0e-08],float)

	v1=.8e-6/1.5
	v2=.3e-8/1.5

	X = array([0,v1,v2,v2,v2,v2],float)
	Ba = X*Ms

	Sr = Ba* costSr
	Y = Ba* costY
	Eu = Ba* costEu

	Elts = ["Ba","Sr","Y","Eu"]

	Ars = [Ba,Sr,Y,Eu]

	# write the table

	outputname = "cescutti2013.txt"

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


	# write yields


	for iM,M in enumerate(Ms):

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



	f.close()

Event Timeline