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

	for Z in Zs:
	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)

	print Z,M,sum(Masses),M-sum(Masses)