compute_HeliumCore.py
View Options

			#!/usr/bin/env python


			import libcris
			import sys,os
			from numpy import *


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

			Ms.sort()
			Zs.sort()

			for Z in Zs:

			outputname = "cristallo-HeliumCore-%6.4f.txt"%Z
			f = open(outputname,'w')
			f.write("# Mass Residu Helium Core\n")

			for M in Ms:

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

			filename = "Total/yields_tot_%s_20130626_14736.txt"%(libcris.Mcode[M] + libcris.Zcode[Z])
			Elts,Masses = libcris.ReadTable(filename)

			Mres = M-sum(Masses)
			Malpha = M-sum(Masses)


			f.write( " %4.2f %13.7e %13.7e\n"%(M,Mres,Malpha) )

Gear/PyChem/shared/tables/AGBs/Cristallo/compute_HeliumCore.pyf7e36798d4e7master