Diffusion Gear (master) Browse

Gear/PyChem/shared/tables/Cescutti/compute_Yields.py
f7e36798d4e7master

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()