In [30]:

    

from pyne import material
from pyne import data
from pyne import nucname

import operator


    

In [31]:

    

loc = 'tot.eq'


    

In [32]:

    

massvec= {}
molvec= {}
for line in file(loc):
    zaid, atoms = line.split()
    molvec[int(zaid)] = float(atoms)
    massvec[int(zaid)] = data.atomic_mass(int(zaid))*float(atoms)
    if atoms==None :
        raise BaseException("The file contains no material data.")
    
print molvec[551370]
print massvec[551370]

print sum([value for key, value in massvec.iteritems()])
print sum([value for key, value in massvec.iteritems() if key!=922350 and key!=922380])


    

    




23.642
3238.954
728900.942949
122292.042949

In [33]:

    

mat = material.Material()
mat.from_atom_frac(molvec)
#print mat
print "There are " + str(len(molvec)) + " radionuclides in this material composition."


    

    




There are 908 radionuclides in this material composition.

In [34]:

    

tot_file = open('tot_eq.csv','w')
tot_file.write("Z, A, Metastable?, ZZAASS, Name, mols, grams \n")
tot = 0
for iso, mols in molvec.iteritems():
    z= str(nucname.znum(iso))
    a= str(nucname.anum(iso))
    s= str(nucname.snum(iso))
    name= nucname.name(iso)
    grams = massvec[iso]
    tot_file.write(z + ', ' +
                   a + ', ' +
                   s + ', ' +
                   str(iso) + ', ' +
                   name + ', ' +
                   str(mols) + ', ' +
                   str(grams) + '\n' )
    tot+=grams


    

In [35]:

    

print tot


    

    




728900.942949

In [6]:

    

activity={}
activity_Ci={}

from scipy import constants

ci_dec_per_sec = 3.7E10 # conversion factor 1 curie = 3.7E10 decays/sec 

for iso, mols in molvec.iteritems():
    dec_per_sec = mols*constants.N_A*data.decay_const(nucname.id(iso))
    activity[nucname.name(iso)] = dec_per_sec
    activity_Ci[nucname.name(iso)]= dec_per_sec/ci_dec_per_sec
    
sorted_a = sorted(activity_Ci.iteritems(), key=operator.itemgetter(1), reverse=True)


    

In [1]:

    

for iso, val in sorted_a : 
    str(iso)+", "+'%.2e' % val


    

    




---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-1-f5ca346430f7> in <module>()
----> 1 tot = file(tot_eq.csv)
      2 for iso, val in sorted_a :
      3     print nucname,id(iso)+", "+str(iso)+", "+'%.2e' % val

NameError: name 'tot_eq' is not defined

In [54]:

    

# quick test - can we get the decay constant of Iodine 130, metastable ?
print nucname.name(531300001)
data.decay_const(531300001)


    

    




I130M






    
Out[54]:





0.0013068385757163373

In [55]:

    

from pyne.transmute.chainsolve import Transmuter


    

In [56]:

    

groups = {}
groups['Fission Products'] = nucname.fp
#groups['Minor Actinides'] = nucname.ma
#groups['Transuranics'] = nucname.tru
groups['Lanthanides'] = lans = nucname.lan
groups['Actinides'] = nucname.act
groups['Fission Products'] = nucname.fp - nucname.lan

elem_act_Ci={}
iso_act_Ci={}
for group, elems in groups.iteritems() :
    #print group
    elem_act_Ci={}
    iso_act_Ci={}
    for num in elems : 
        elem = nucname.name(num)
        elem_act_Ci[elem] = sum([value for key, value in activity_Ci.items() if num == nucname.znum(key)])
        for key, value in activity_Ci.items():
            if num== nucname.znum(key):
                iso_act_Ci[key] = value
    sorted_a = sorted(elem_act_Ci.iteritems(), key=operator.itemgetter(1),reverse=True)
    for entry in sorted_a :  
        elem, val = entry
        #if elem_act_Ci[elem] > 0:
            #if len(elem) == 1:
                #print elem + "      " + '%.2e' % elem_act_Ci[elem]
            #else:
                #print elem + "     " + '%.2e' % elem_act_Ci[elem]
    sorted_iso_a = sorted(iso_act_Ci.iteritems(), key=operator.itemgetter(1),reverse=True)
    i=0
    for entry in sorted_iso_a :  
        iso, val = entry
        #if iso_act_Ci[iso] > 0:
        nspaces=10-len(nucname.name(iso))
        i+=1
        print nucname.name(iso) + ","+" " + '%.2e' % iso_act_Ci[iso]
    print i


    

    




Sb134M, 5.09e+15
As88, 1.63e+09
Y96, 2.67e+08
Nb103, 5.05e+07
Zr102, 3.17e+07
Sr96, 2.28e+07
Tc109, 2.08e+07
I134, 1.41e+07
I133, 1.25e+07
Xe133, 1.22e+07
Zr99, 1.20e+07
I135, 1.18e+07
Cs139, 1.16e+07
Ba139, 1.15e+07
Te134, 1.14e+07
Cs138, 1.13e+07
Xe137, 1.12e+07
Y95, 1.11e+07
Xe138, 1.11e+07
Mo99, 1.10e+07
Ba140, 1.10e+07
Zr97, 1.07e+07
Zr98, 1.07e+07
Nb97, 1.07e+07
Y94, 1.07e+07
Sb122M, 1.07e+07
Ba141, 1.06e+07
Nb98, 1.05e+07
Sr93, 1.05e+07
Y93, 1.05e+07
Zr95, 1.04e+07
Ba142, 1.03e+07
Nb99, 1.03e+07
Nb97M, 1.03e+07
Cs140, 1.03e+07
Nb95, 1.01e+07
Br95, 1.01e+07
Tc101, 1.00e+07
Mo101, 1.00e+07
Sr94, 9.77e+06
Zr100, 9.69e+06
Sr95, 9.68e+06
Sr92, 9.66e+06
Tc99M, 9.64e+06
Nb100, 9.50e+06
Nb101, 9.39e+06
Y92, 9.39e+06
Tc102, 8.94e+06
Sr91, 8.90e+06
Mo104, 8.81e+06
Zr101, 8.80e+06
Mo102, 8.78e+06
Y91, 8.73e+06
Xe139, 8.71e+06
Ba143, 8.57e+06
Rb91, 8.50e+06
I132, 8.42e+06
Te132, 8.34e+06
Cs141, 7.76e+06
Ru103, 7.50e+06
Rb92, 7.28e+06
Te133, 7.25e+06
Tc103, 7.23e+06
Rb89, 7.22e+06
Sb134, 7.01e+06
Sr89, 6.96e+06
Kr89, 6.90e+06
Mo103, 6.88e+06
Rb90, 6.87e+06
Kr90, 6.79e+06
Ba144, 6.79e+06
Rh103M, 6.75e+06
Xe140, 5.85e+06
Tc104, 5.80e+06
Te135, 5.73e+06
I137, 5.68e+06
I131, 5.66e+06
Rb88, 5.64e+06
I136, 5.54e+06
Kr88, 5.50e+06
Mo105, 5.37e+06
Rb93, 5.37e+06
Nb102, 5.34e+06
Y91M, 5.28e+06
Te131, 5.16e+06
Te133M, 5.06e+06
Kr91, 5.03e+06
Sb131, 4.94e+06
Ba145, 4.90e+06
Cs142, 4.81e+06
Nb100M, 4.79e+06
Tc105, 4.53e+06
Ru105, 4.32e+06
Kr87, 3.89e+06
Te136, 3.64e+06
Br88, 3.64e+06
Sb133, 3.60e+06
I136M, 3.45e+06
Rh105, 3.39e+06
Y98, 3.37e+06
Br87, 3.28e+06
Sb132, 3.02e+06
Tc106, 2.86e+06
I138, 2.85e+06
Rb94, 2.74e+06
Br89, 2.71e+06
Y97, 2.59e+06
Sb130M, 2.58e+06
Ru107, 2.48e+06
Cs143, 2.30e+06
Kr92, 2.20e+06
Rh107, 2.20e+06
Se86, 2.07e+06
Sb132M, 2.05e+06
Rh106, 2.04e+06
Sn131, 2.03e+06
Kr85M, 1.97e+06
Sn130, 1.93e+06
Br85, 1.93e+06
Xe141, 1.92e+06
Y99, 1.84e+06
Mo106, 1.82e+06
Sn129, 1.78e+06
Tc100, 1.62e+06
Te129, 1.60e+06
Sb129, 1.59e+06
Rb90M, 1.57e+06
Br84, 1.57e+06
Se84, 1.55e+06
Rh108, 1.51e+06
Y100, 1.50e+06
Ge88, 1.49e+06
Xe135, 1.49e+06
Ru108, 1.48e+06
Br90, 1.45e+06
Sr97, 1.45e+06
Sr98, 1.44e+06
Br86, 1.44e+06
Ru106, 1.40e+06
Zr103, 1.39e+06
Rb95, 1.29e+06
Rh111, 1.28e+06
Rh105M, 1.27e+06
I139, 1.26e+06
I134M, 1.20e+06
Mo107, 1.18e+06
Se85, 1.16e+06
Ba146, 1.15e+06
Ru111, 1.13e+06
Ag109M, 1.07e+06
Rh109, 1.05e+06
Pd109, 1.05e+06
Sn132, 9.96e+05
Te137, 9.91e+05
Y101, 9.51e+05
Ru109, 9.11e+05
Sb128M, 8.83e+05
Kr83M, 8.70e+05
Br83, 8.64e+05
Tc108, 8.27e+05
Sn128, 8.18e+05
Kr93, 7.39e+05
Te131M, 7.13e+05
As84, 6.46e+05
As83, 6.25e+05
Xe142, 6.02e+05
Br91, 5.96e+05
In128, 5.93e+05
Zr104, 5.80e+05
Cs144, 5.66e+05
Ba147, 5.56e+05
Sb130, 5.52e+05
Ru110, 5.30e+05
Se83M, 5.25e+05
Cs138M, 5.20e+05
Sr99, 5.09e+05
Se88, 5.04e+05
Pd109M, 4.70e+05
Cs134, 4.53e+05
Cd124, 4.51e+05
Sb127, 4.46e+05
Te127, 4.37e+05
Nb104, 4.12e+05
Nb105, 3.75e+05
Xe133M, 3.74e+05
Se81, 3.62e+05
As82, 3.60e+05
Ag110, 3.44e+05
In120, 3.35e+05
Se83, 3.35e+05
I140, 3.27e+05
As81, 3.27e+05
I133M, 3.26e+05
Rb96, 3.20e+05
Br96, 3.20e+05
Kr94, 3.16e+05
Sb135, 3.12e+05
Rh106M, 3.01e+05
In122, 2.99e+05
As85, 2.96e+05
Ge81, 2.91e+05
Sn127, 2.83e+05
Cs137, 2.81e+05
Sn133, 2.72e+05
Ba137M, 2.67e+05
Ba148, 2.60e+05
Sn129M, 2.52e+05
In130, 2.44e+05
As80, 2.35e+05
Te129M, 2.32e+05
Y90, 2.32e+05
Sr90, 2.27e+05
Ag111M, 2.27e+05
Ag111, 2.21e+05
Pd111, 2.13e+05
Mo108, 2.07e+05
Sr100, 1.99e+05
Ge82, 1.97e+05
Cs134M, 1.78e+05
Ge79, 1.75e+05
Rh104M, 1.69e+05
Te138, 1.64e+05
As86, 1.50e+05
Sn127M, 1.47e+05
Ge80, 1.45e+05
Cd123, 1.44e+05
Ge83, 1.43e+05
Cd126, 1.37e+05
Se89, 1.32e+05
In127, 1.28e+05
As82M, 1.26e+05
In129, 1.20e+05
Cs145, 1.19e+05
In126, 1.14e+05
Ag112, 1.10e+05
Cs136, 1.04e+05
Pd112, 1.04e+05
Nb106, 1.03e+05
Rh115, 9.60e+04
As79, 9.29e+04
Se79M, 9.24e+04
Xe143, 8.52e+04
Nb98M, 8.24e+04
Y102, 8.21e+04
Pd113, 8.12e+04
Cd119, 7.99e+04
Ag113, 7.46e+04
Kr95, 7.35e+04
Sn125M, 7.35e+04
Nb95M, 7.31e+04
Xe134M, 7.31e+04
I130, 7.13e+04
Ru113, 7.08e+04
Pd115, 7.08e+04
In127M, 6.91e+04
In124, 6.85e+04
Mo109, 6.43e+04
Tc110, 6.37e+04
Tc111, 6.36e+04
Xe131M, 6.31e+04
Nb107, 6.28e+04
Sb128, 6.25e+04
Cs135M, 5.89e+04
Ag114, 5.80e+04
Pd114, 5.74e+04
Rb97, 5.73e+04
Ru114, 5.69e+04
Cd125, 5.44e+04
Te127M, 5.41e+04
Sn121, 4.76e+04
In118, 4.66e+04
Cd118, 4.66e+04
As87, 4.66e+04
In121, 4.60e+04
Ag121, 4.57e+04
In131, 4.55e+04
Sn123M, 4.53e+04
Cd115, 4.52e+04
Pd116, 4.52e+04
Cd120, 4.51e+04
Sn125, 4.41e+04
In115M, 4.33e+04
Br84M, 4.27e+04
I141, 4.26e+04
Cd122, 4.19e+04
Cd121, 4.13e+04
In125, 4.02e+04
Pd117, 4.00e+04
As78, 3.86e+04
Ge78, 3.74e+04
Ge84, 3.64e+04
In123, 3.59e+04
Ag115, 3.56e+04
In119M, 3.39e+04
In117M, 3.37e+04
Cd119M, 3.32e+04
Ga80, 3.25e+04
In125M, 3.15e+04
I128, 3.00e+04
Cd117, 2.99e+04
Ag118, 2.97e+04
Rh114, 2.96e+04
Kr85, 2.88e+04
Pd119, 2.86e+04
I130M, 2.86e+04
Ga78, 2.77e+04
Pd118, 2.76e+04
In117, 2.71e+04
Ga79, 2.63e+04
Sb125, 2.61e+04
Ag118M, 2.61e+04
Br92, 2.56e+04
Pd120, 2.49e+04
Zr105, 2.30e+04
Rh117, 2.29e+04
Ag117, 2.17e+04
Ag117M, 2.15e+04
Cd128, 2.08e+04
Nb96, 2.08e+04
Rh113, 2.07e+04
Ru112, 2.05e+04
Cd127, 2.01e+04
Ag120, 2.01e+04
Mo110, 1.94e+04
In123M, 1.79e+04
As77, 1.74e+04
In116, 1.73e+04
Sn134, 1.70e+04
Ag116, 1.59e+04
Cd117M, 1.58e+04
Y103, 1.56e+04
Ge77M, 1.46e+04
Ga77, 1.46e+04
Kr96, 1.43e+04
Rh108M, 1.38e+04
In119, 1.37e+04
In116M, 1.35e+04
Ag115M, 1.33e+04
Br93, 1.33e+04
Sb136, 1.19e+04
Ag116M, 1.18e+04
Br82, 1.17e+04
Sn123, 1.12e+04
Se81M, 1.12e+04
Zn78, 1.03e+04
Sr101, 1.01e+04
Rb98, 1.01e+04
Ga81, 9.64e+03
Xe144, 9.36e+03
In121M, 8.83e+03
Ag110M, 8.46e+03
Cs146, 8.40e+03
Ag123, 8.30e+03
Ag113M, 8.21e+03
Rh116, 8.12e+03
Rb86, 7.82e+03
Tc102M, 7.52e+03
Ge77, 6.30e+03
Ga76, 6.29e+03
In132, 5.93e+03
Te125M, 4.85e+03
Zn76, 4.81e+03
Ba149, 4.78e+03
Br82M, 4.71e+03
Cd130, 4.52e+03
Cd115M, 4.29e+03
Ge85, 3.97e+03
Zn77, 3.86e+03
Sb126, 3.55e+03
Pd111M, 3.47e+03
I142, 3.45e+03
Rh110M, 3.34e+03
Se91, 3.02e+03
Tc112, 2.98e+03
Ge75, 2.86e+03
Ga75, 2.57e+03
Cs147, 2.56e+03
Sb122, 2.54e+03
Sb137, 2.27e+03
Nb108, 2.21e+03
Ag124, 2.12e+03
Zn75, 2.10e+03
Sb126M, 2.06e+03
Sr102, 1.99e+03
Xe145, 1.84e+03
Pd121, 1.80e+03
Sb124, 1.62e+03
In120M, 1.51e+03
Ru115, 1.46e+03
Zn80, 1.42e+03
Ga82, 1.40e+03
Rh118, 1.38e+03
H3, 1.37e+03
Rh119, 1.29e+03
Ga74, 1.23e+03
Zn74, 1.15e+03
Ag122, 1.15e+03
Sn135, 1.09e+03
Rb99, 9.82e+02
Ag125, 9.44e+02
Zn79, 9.40e+02
Rb86M, 9.29e+02
Br94, 7.59e+02
In122M, 7.24e+02
Ga83, 6.47e+02
Ge73M, 5.54e+02
Ba150, 5.42e+02
Ru118, 5.29e+02
Ga73, 5.24e+02
Y104, 5.01e+02
Zn73, 4.92e+02
In133, 4.45e+02
Sn119M, 3.84e+02
Zr107, 3.34e+02
Cu75, 2.75e+02
Cd131, 2.44e+02
Nb109, 2.32e+02
Ge87, 2.21e+02
Ag126, 2.17e+02
Cu73, 1.98e+02
Zn72, 1.90e+02
Ga72, 1.89e+02
Kr97, 1.87e+02
Pd107M, 1.60e+02
Sn136, 1.55e+02
Cu74, 1.49e+02
Pd124, 1.42e+02
Cu76, 1.41e+02
Cu72, 1.40e+02
Ge75M, 1.29e+02
Cu77, 1.09e+02
Ga84, 9.95e+01
Ni72, 9.44e+01
Cd113M, 7.92e+01
Zn81, 7.66e+01
Tc114, 7.64e+01
Ba135M, 7.01e+01
As76, 6.82e+01
Y90M, 6.50e+01
Rb100, 5.92e+01
Sn117M, 5.85e+01
Cs148, 5.06e+01
Se77M, 4.43e+01
Cd111M, 4.33e+01
Ni74, 3.70e+01
Sb124M, 3.51e+01
Tc99, 3.44e+01
In118M, 2.72e+01
Rh120, 2.52e+01
Cu78, 2.32e+01
Cd132, 1.56e+01
Sr87M, 1.53e+01
Ni73, 1.48e+01
Ag128, 1.13e+01
Te123M, 9.75e+00
In134, 8.00e+00
Xe129M, 6.68e+00
Cu79, 6.25e+00
Nb110, 5.67e+00
Ni76, 5.56e+00
Zr93, 5.18e+00
Ni75, 4.64e+00
Nb94M, 4.29e+00
Rh121, 4.00e+00
Co72, 3.72e+00
Br80, 3.47e+00
In114, 2.51e+00
Br80M, 2.10e+00
Y89M, 2.06e+00
In114M, 9.15e-01
Sn126, 6.11e-01
Sn121M, 5.23e-01
Cs135, 5.07e-01
Co74, 4.52e-01
Ag108, 4.41e-01
Se79, 2.36e-01
Pd107, 2.11e-01
Y107, 1.89e-01
Nb93M, 1.28e-01
Kr81M, 9.65e-02
Ni78, 9.58e-02
I129, 7.64e-02
Br79M, 6.22e-02
Cd109, 2.02e-02
Zr90M, 1.72e-02
Cu66, 8.11e-04
C14, 3.13e-04
Nb94, 2.73e-04
Rn220, 1.33e-04
Ra224, 1.33e-04
Rb87, 6.22e-05
In113M, 1.89e-05
Ag108M, 1.43e-05
Be10, 8.21e-06
Kr81, 6.77e-07
Ra223, 1.91e-08
Pb210, 1.04e-10
Bi210, 1.01e-10
Po210, 4.54e-11
In115, 3.28e-11
Ra226, 6.28e-12
Rn222, 6.02e-12
Ra228, 1.43e-12
Cd113, 9.11e-14
Te123, 1.51e-16
Sb123, 0.00e+00
Sb121, 0.00e+00
Cd108, 0.00e+00
Se78, 0.00e+00
Se77, 0.00e+00
Se76, 0.00e+00
Rh110, 0.00e+00
Rh112, 0.00e+00
Pd110, 0.00e+00
Xe129, 0.00e+00
Pb207, 0.00e+00
Pb206, 0.00e+00
Pb208, 0.00e+00
Y89, 0.00e+00
Tc107, 0.00e+00
Cu80, 0.00e+00
Cu81, 0.00e+00
I127, 0.00e+00
Co73, 0.00e+00
Nb99M, 0.00e+00
Ag107, 0.00e+00
Be9, 0.00e+00
Pd122, 0.00e+00
Pd123, 0.00e+00
Sn116, 0.00e+00
Ru120, 0.00e+00
Sn114, 0.00e+00
Pd126, 0.00e+00
Sn115, 0.00e+00
Xe128, 0.00e+00
Zn66, 0.00e+00
Zn67, 0.00e+00
Zn68, 0.00e+00
Ga69, 0.00e+00
Tc117, 0.00e+00
Te130, 0.00e+00
Te139, 0.00e+00
Li6, 0.00e+00
Li7, 0.00e+00
Ge86, 0.00e+00
Rb85, 0.00e+00
Br86M, 0.00e+00
Cd114, 0.00e+00
Cd116, 0.00e+00
Cd110, 0.00e+00
Cd111, 0.00e+00
Cd112, 0.00e+00
Rh104, 0.00e+00
Rh103, 0.00e+00
Mo113, 0.00e+00
Mo112, 0.00e+00
Mo111, 0.00e+00
Mo115, 0.00e+00
Mo114, 0.00e+00
Pd102, 0.00e+00
Pd105, 0.00e+00
Pd106, 0.00e+00
Pd108, 0.00e+00
Zr92, 0.00e+00
Zr91, 0.00e+00
Zr90, 0.00e+00
Zr96, 0.00e+00
Zr94, 0.00e+00
Se85M, 0.00e+00
Tc113, 0.00e+00
Tc116, 0.00e+00
Ru99, 0.00e+00
Tc115, 0.00e+00
Te142, 0.00e+00
Te141, 0.00e+00
Te140, 0.00e+00
Cs133, 0.00e+00
Nb111, 0.00e+00
Ba138, 0.00e+00
As75, 0.00e+00
Ba132, 0.00e+00
Ba134, 0.00e+00
Ba135, 0.00e+00
Ba136, 0.00e+00
Ba137, 0.00e+00
Xe147, 0.00e+00
Xe146, 0.00e+00
Cd129, 0.00e+00
Xe135M, 0.00e+00
Zn70, 0.00e+00
Sr104, 0.00e+00
Co75, 0.00e+00
Ga71, 0.00e+00
Sn117, 0.00e+00
Sn118, 0.00e+00
Sn119, 0.00e+00
Se87, 0.00e+00
Ru104, 0.00e+00
I143, 0.00e+00
I144, 0.00e+00
Xe126, 0.00e+00
Mo100, 0.00e+00
Y105, 0.00e+00
Se80, 0.00e+00
Se82, 0.00e+00
Sr103, 0.00e+00
Ru101, 0.00e+00
Ru100, 0.00e+00
Ru102, 0.00e+00
Xe132, 0.00e+00
Cs150, 0.00e+00
As89, 0.00e+00
Sb138, 0.00e+00
Kr98, 0.00e+00
Ba152, 0.00e+00
Rh109M, 0.00e+00
Zn82, 0.00e+00
Zn83, 0.00e+00
Sb139, 0.00e+00
Sr88, 0.00e+00
Sr87, 0.00e+00
Sr86, 0.00e+00
In113, 0.00e+00
Bi209, 0.00e+00
Ag109, 0.00e+00
Sn124, 0.00e+00
Sn122, 0.00e+00
Sn120, 0.00e+00
Xe136, 0.00e+00
Ge70, 0.00e+00
Xe134, 0.00e+00
Xe130, 0.00e+00
Xe131, 0.00e+00
Pd104, 0.00e+00
Ag119, 0.00e+00
Rh122, 0.00e+00
Rh123, 0.00e+00
Se92, 0.00e+00
Se90, 0.00e+00
Ru116, 0.00e+00
Ru117, 0.00e+00
Ni77, 0.00e+00
Mo97, 0.00e+00
Mo96, 0.00e+00
Mo95, 0.00e+00
Mo98, 0.00e+00
Te122, 0.00e+00
Te125, 0.00e+00
Te124, 0.00e+00
Te126, 0.00e+00
Te128, 0.00e+00
Nb93, 0.00e+00
Zr108, 0.00e+00
Zr109, 0.00e+00
Ge73, 0.00e+00
Ge72, 0.00e+00
Ge74, 0.00e+00
Ge76, 0.00e+00
Zr106, 0.00e+00
Br81, 0.00e+00
Kr86, 0.00e+00
Kr84, 0.00e+00
Kr83, 0.00e+00
Kr82, 0.00e+00
Kr80, 0.00e+00
Br79, 0.00e+00
666
Ho166M, 3.48e+11
Nd158, 1.04e+11
Ce154, 8.57e+09
La152, 2.35e+08
La140, 1.11e+07
La141, 1.06e+07
La142, 1.05e+07
Ce141, 1.01e+07
Ce143, 1.01e+07
La143, 1.01e+07
Pr143, 9.86e+06
La144, 8.97e+06
La145, 7.48e+06
Pr145, 6.84e+06
Ce145, 6.83e+06
Pr144, 5.70e+06
Ce146, 5.58e+06
Ce144, 5.53e+06
La146, 5.34e+06
Pr146, 5.17e+06
Ce147, 5.03e+06
La147, 4.65e+06
Nd147, 3.89e+06
Pr147, 3.72e+06
Pm149, 2.73e+06
Pr150, 2.58e+06
Ce148, 2.19e+06
Nd149, 2.18e+06
Pr149, 2.15e+06
Sm153, 1.63e+06
Nd151, 9.98e+05
Pm151, 9.84e+05
Nd153, 8.32e+05
Pm152, 6.98e+05
Nd152, 6.84e+05
La148, 6.63e+05
Pr152, 6.59e+05
Eu156, 6.31e+05
Pm147, 5.20e+05
Pm153, 4.54e+05
Pm148M, 4.08e+05
La149, 3.63e+05
Pm154, 3.62e+05
Ce152, 3.00e+05
Ce149, 2.83e+05
Nd154, 2.81e+05
Sm158, 2.65e+05
Nd156, 2.13e+05
Sm159, 2.03e+05
Pm157, 1.99e+05
Nd155, 1.90e+05
Sm160, 1.87e+05
Pr142, 1.73e+05
Sm155, 1.56e+05
Ce150, 1.54e+05
Pr153, 1.45e+05
Pr151, 1.44e+05
Pm155, 1.18e+05
Pm150, 1.14e+05
Eu157, 1.10e+05
Ce151, 9.42e+04
Sm156, 9.13e+04
Pr144M, 6.66e+04
Sm157, 5.82e+04
Pr142M, 3.52e+04
Eu158, 3.45e+04
Pm156, 3.33e+04
Pm154M, 2.55e+04
Gd159, 2.15e+04
Eu162, 2.15e+04
Eu159, 1.74e+04
La150, 1.53e+04
Pm152M, 1.22e+04
Eu154, 1.20e+04
Eu160, 1.06e+04
Pr154, 9.19e+03
Pm158, 7.85e+03
Eu155, 7.43e+03
Gd164, 6.44e+03
Pm159, 5.98e+03
Eu161, 5.20e+03
Tb161, 3.93e+03
Gd161, 3.79e+03
Sm161, 3.39e+03
Gd162, 1.97e+03
Tb162, 1.61e+03
Tb160, 1.53e+03
Sm151, 1.06e+03
Sm162, 1.00e+03
Gd163, 8.42e+02
Tb163, 6.73e+02
Dy165, 6.40e+02
Gd165, 6.31e+02
Eu163, 3.22e+02
Tb164, 2.89e+02
Dy165M, 2.76e+02
Eu152M, 2.35e+02
Ho166, 1.69e+02
Dy166, 8.36e+01
Gd153, 4.05e+01
Tb165, 3.10e+01
Eu164, 1.36e+01
Eu152, 1.03e+01
Eu165, 3.37e+00
Er167M, 2.15e+00
Tm170, 1.23e-02
Tm171, 1.70e-04
Sm147, 1.87e-06
Er171, 3.36e-08
La138, 2.94e-09
Nd144, 1.92e-09
Sm148, 1.08e-10
Gd152, 3.88e-13
Sm152, 0.00e+00
Sm150, 0.00e+00
La139, 0.00e+00
Sm154, 0.00e+00
Yb170, 0.00e+00
Dy162, 0.00e+00
Yb171, 0.00e+00
Dy164, 0.00e+00
Dy161, 0.00e+00
Gd158, 0.00e+00
Dy160, 0.00e+00
Gd154, 0.00e+00
Gd155, 0.00e+00
Gd156, 0.00e+00
Gd157, 0.00e+00
Sm164, 0.00e+00
Sm165, 0.00e+00
Er170, 0.00e+00
Sm163, 0.00e+00
Pr155, 0.00e+00
Pr157, 0.00e+00
Pr156, 0.00e+00
Pr159, 0.00e+00
Pr158, 0.00e+00
Nd159, 0.00e+00
Nd150, 0.00e+00
Nd157, 0.00e+00
Tm169, 0.00e+00
Tb162M, 0.00e+00
Pr141, 0.00e+00
Yb172, 0.00e+00
Pr148, 0.00e+00
Nd148, 0.00e+00
Nd142, 0.00e+00
Nd143, 0.00e+00
Nd146, 0.00e+00
Nd145, 0.00e+00
Sm149, 0.00e+00
Ce153, 0.00e+00
Ce157, 0.00e+00
Ce156, 0.00e+00
Ce155, 0.00e+00
Ho165, 0.00e+00
Pm161, 0.00e+00
Pm160, 0.00e+00
Pm162, 0.00e+00
Nd160, 0.00e+00
Nd161, 0.00e+00
Eu151, 0.00e+00
Eu153, 0.00e+00
Yb168, 0.00e+00
La151, 0.00e+00
La153, 0.00e+00
La154, 0.00e+00
Tb159, 0.00e+00
Sm145, 0.00e+00
Dy163, 0.00e+00
Gd160, 0.00e+00
Er168, 0.00e+00
Er167, 0.00e+00
Er166, 0.00e+00
Pm148, 0.00e+00
Ce140, 0.00e+00
Ce142, 0.00e+00
177
U239, 7.83e+07
Np239, 7.57e+07
U237, 3.91e+06
Np238, 1.15e+06
Pu243, 9.65e+05
Pu241, 2.73e+05
Am244M, 2.52e+05
Am242, 9.12e+04
Cm242, 3.81e+04
Am244, 1.30e+04
Pu238, 5.75e+03
Cm244, 5.70e+03
U240, 1.68e+03
Pu240, 6.75e+02
Pu239, 4.51e+02
Am241, 1.45e+02
Am243, 4.79e+01
Cm243, 1.24e+01
Am242M, 1.09e+01
Pu245, 9.44e+00
Pu242, 5.05e+00
U236, 9.67e-01
Pu237, 8.54e-01
Np237, 8.40e-01
Pa233, 7.47e-01
Cm245, 7.20e-01
Pu236, 1.86e-01
Pa234M, 1.84e-01
U238, 1.84e-01
Th234, 1.72e-01
U235, 1.30e-01
Cm246, 1.07e-01
Pa232, 4.38e-02
Cm249, 2.33e-02
Th233, 1.80e-02
Bk250, 1.38e-02
U234, 1.12e-02
Pa234, 1.03e-02
Bk249, 2.32e-03
Np235, 1.40e-03
U232, 1.20e-03
Th228, 1.38e-04
Cf252, 1.32e-04
Cf250, 1.17e-04
Cf253, 5.14e-06
U233, 4.05e-06
Es253, 4.01e-06
Pu244, 2.63e-06
Pa231, 1.48e-06
U231, 1.27e-06
Cf249, 1.23e-06
Np236, 1.14e-06
Cf251, 1.07e-06
Cm248, 8.55e-07
Cm247, 3.14e-07
Es254M, 1.42e-07
Cf254, 6.96e-08
Th230, 4.29e-08
Es254, 3.79e-08
Ac227, 2.36e-08
Th227, 2.05e-08
U230, 5.61e-09
Th229, 2.12e-09
Th232, 2.91e-11
Cm250, 3.44e-13
65

In [57]:

    

for key, val in mat.comp.iteritems():
    print nucname.name(key)


    

    




H3
Li6
Li7
Be9
Be10
C14
Co72
Co73
Co74
Co75
Ni72
Ni73
Ni74
Ni75
Ni76
Ni77
Ni78
Cu66
Cu72
Cu73
Cu74
Cu75
Cu76
Cu77
Cu78
Cu79
Cu80
Cu81
Zn66
Zn67
Zn68
Zn70
Zn72
Zn73
Zn74
Zn75
Zn76
Zn77
Zn78
Zn79
Zn80
Zn81
Zn82
Zn83
Ga69
Ga71
Ga72
Ga73
Ga74
Ga75
Ga76
Ga77
Ga78
Ga79
Ga80
Ga81
Ga82
Ga83
Ga84
Ge70
Ge72
Ge73
Ge73M
Ge74
Ge75
Ge75M
Ge76
Ge77
Ge77M
Ge78
Ge79
Ge80
Ge81
Ge82
Ge83
Ge84
Ge85
Ge86
Ge87
Ge88
As75
As76
As77
As78
As79
As80
As81
As82
As82M
As83
As84
As85
As86
As87
As88
As89
Se76
Se77
Se77M
Se78
Se79
Se79M
Se80
Se81
Se81M
Se82
Se83
Se83M
Se84
Se85
Se85M
Se86
Se87
Se88
Se89
Se90
Se91
Se92
Br79
Br79M
Br80
Br80M
Br81
Br82
Br82M
Br83
Br84
Br84M
Br85
Br86
Br86M
Br87
Br88
Br89
Br90
Br91
Br92
Br93
Br94
Br95
Br96
Kr80
Kr81
Kr81M
Kr82
Kr83
Kr83M
Kr84
Kr85
Kr85M
Kr86
Kr87
Kr88
Kr89
Kr90
Kr91
Kr92
Kr93
Kr94
Kr95
Kr96
Kr97
Kr98
Rb85
Rb86
Rb86M
Rb87
Rb88
Rb89
Rb90
Rb90M
Rb91
Rb92
Rb93
Rb94
Rb95
Rb96
Rb97
Rb98
Rb99
Rb100
Sr86
Sr87
Sr87M
Sr88
Sr89
Sr90
Sr91
Sr92
Sr93
Sr94
Sr95
Sr96
Sr97
Sr98
Sr99
Sr100
Sr101
Sr102
Sr103
Sr104
Y89
Y89M
Y90
Y90M
Y91
Y91M
Y92
Y93
Y94
Y95
Y96
Y97
Y98
Y99
Y100
Y101
Y102
Y103
Y104
Y105
Y107
Zr90
Zr90M
Zr91
Zr92
Zr93
Zr94
Zr95
Zr96
Zr97
Zr98
Zr99
Zr100
Zr101
Zr102
Zr103
Zr104
Zr105
Zr106
Zr107
Zr108
Zr109
Nb93
Nb93M
Nb94
Nb94M
Nb95
Nb95M
Nb96
Nb97
Nb97M
Nb98
Nb98M
Nb99
Nb99M
Nb100
Nb100M
Nb101
Nb102
Nb103
Nb104
Nb105
Nb106
Nb107
Nb108
Nb109
Nb110
Nb111
Mo95
Mo96
Mo97
Mo98
Mo99
Mo100
Mo101
Mo102
Mo103
Mo104
Mo105
Mo106
Mo107
Mo108
Mo109
Mo110
Mo111
Mo112
Mo113
Mo114
Mo115
Tc99
Tc99M
Tc100
Tc101
Tc102
Tc102M
Tc103
Tc104
Tc105
Tc106
Tc107
Tc108
Tc109
Tc110
Tc111
Tc112
Tc113
Tc114
Tc115
Tc116
Tc117
Ru99
Ru100
Ru101
Ru102
Ru103
Ru104
Ru105
Ru106
Ru107
Ru108
Ru109
Ru110
Ru111
Ru112
Ru113
Ru114
Ru115
Ru116
Ru117
Ru118
Ru120
Rh103
Rh103M
Rh104
Rh104M
Rh105
Rh105M
Rh106
Rh106M
Rh107
Rh108
Rh108M
Rh109
Rh109M
Rh110
Rh110M
Rh111
Rh112
Rh113
Rh114
Rh115
Rh116
Rh117
Rh118
Rh119
Rh120
Rh121
Rh122
Rh123
Pd102
Pd104
Pd105
Pd106
Pd107
Pd107M
Pd108
Pd109
Pd109M
Pd110
Pd111
Pd111M
Pd112
Pd113
Pd114
Pd115
Pd116
Pd117
Pd118
Pd119
Pd120
Pd121
Pd122
Pd123
Pd124
Pd126
Ag107
Ag108
Ag108M
Ag109
Ag109M
Ag110
Ag110M
Ag111
Ag111M
Ag112
Ag113
Ag113M
Ag114
Ag115
Ag115M
Ag116
Ag116M
Ag117
Ag117M
Ag118
Ag118M
Ag119
Ag120
Ag121
Ag122
Ag123
Ag124
Ag125
Ag126
Ag128
Cd108
Cd109
Cd110
Cd111
Cd111M
Cd112
Cd113
Cd113M
Cd114
Cd115
Cd115M
Cd116
Cd117
Cd117M
Cd118
Cd119
Cd119M
Cd120
Cd121
Cd122
Cd123
Cd124
Cd125
Cd126
Cd127
Cd128
Cd129
Cd130
Cd131
Cd132
In113
In113M
In114
In114M
In115
In115M
In116
In116M
In117
In117M
In118
In118M
In119
In119M
In120
In120M
In121
In121M
In122
In122M
In123
In123M
In124
In125
In125M
In126
In127
In127M
In128
In129
In130
In131
In132
In133
In134
Sn114
Sn115
Sn116
Sn117
Sn117M
Sn118
Sn119
Sn119M
Sn120
Sn121
Sn121M
Sn122
Sn123
Sn123M
Sn124
Sn125
Sn125M
Sn126
Sn127
Sn127M
Sn128
Sn129
Sn129M
Sn130
Sn131
Sn132
Sn133
Sn134
Sn135
Sn136
Sb121
Sb122
Sb122M
Sb123
Sb124
Sb124M
Sb125
Sb126
Sb126M
Sb127
Sb128
Sb128M
Sb129
Sb130
Sb130M
Sb131
Sb132
Sb132M
Sb133
Sb134
Sb134M
Sb135
Sb136
Sb137
Sb138
Sb139
Te122
Te123
Te123M
Te124
Te125
Te125M
Te126
Te127
Te127M
Te128
Te129
Te129M
Te130
Te131
Te131M
Te132
Te133
Te133M
Te134
Te135
Te136
Te137
Te138
Te139
Te140
Te141
Te142
I127
I128
I129
I130
I130M
I131
I132
I133
I133M
I134
I134M
I135
I136
I136M
I137
I138
I139
I140
I141
I142
I143
I144
Xe126
Xe128
Xe129
Xe129M
Xe130
Xe131
Xe131M
Xe132
Xe133
Xe133M
Xe134
Xe134M
Xe135
Xe135M
Xe136
Xe137
Xe138
Xe139
Xe140
Xe141
Xe142
Xe143
Xe144
Xe145
Xe146
Xe147
Cs133
Cs134
Cs134M
Cs135
Cs135M
Cs136
Cs137
Cs138
Cs138M
Cs139
Cs140
Cs141
Cs142
Cs143
Cs144
Cs145
Cs146
Cs147
Cs148
Cs150
Ba132
Ba134
Ba135
Ba135M
Ba136
Ba137
Ba137M
Ba138
Ba139
Ba140
Ba141
Ba142
Ba143
Ba144
Ba145
Ba146
Ba147
Ba148
Ba149
Ba150
Ba152
La138
La139
La140
La141
La142
La143
La144
La145
La146
La147
La148
La149
La150
La151
La152
La153
La154
Ce140
Ce141
Ce142
Ce143
Ce144
Ce145
Ce146
Ce147
Ce148
Ce149
Ce150
Ce151
Ce152
Ce153
Ce154
Ce155
Ce156
Ce157
Pr141
Pr142
Pr142M
Pr143
Pr144
Pr144M
Pr145
Pr146
Pr147
Pr148
Pr149
Pr150
Pr151
Pr152
Pr153
Pr154
Pr155
Pr156
Pr157
Pr158
Pr159
Nd142
Nd143
Nd144
Nd145
Nd146
Nd147
Nd148
Nd149
Nd150
Nd151
Nd152
Nd153
Nd154
Nd155
Nd156
Nd157
Nd158
Nd159
Nd160
Nd161
Pm147
Pm148
Pm148M
Pm149
Pm150
Pm151
Pm152
Pm152M
Pm153
Pm154
Pm154M
Pm155
Pm156
Pm157
Pm158
Pm159
Pm160
Pm161
Pm162
Sm145
Sm147
Sm148
Sm149
Sm150
Sm151
Sm152
Sm153
Sm154
Sm155
Sm156
Sm157
Sm158
Sm159
Sm160
Sm161
Sm162
Sm163
Sm164
Sm165
Eu151
Eu152
Eu152M
Eu153
Eu154
Eu155
Eu156
Eu157
Eu158
Eu159
Eu160
Eu161
Eu162
Eu163
Eu164
Eu165
Gd152
Gd153
Gd154
Gd155
Gd156
Gd157
Gd158
Gd159
Gd160
Gd161
Gd162
Gd163
Gd164
Gd165
Tb159
Tb160
Tb161
Tb162
Tb162M
Tb163
Tb164
Tb165
Dy160
Dy161
Dy162
Dy163
Dy164
Dy165
Dy165M
Dy166
Ho165
Ho166
Ho166M
Er166
Er167
Er167M
Er168
Er170
Er171
Tm169
Tm170
Tm171
Yb168
Yb170
Yb171
Yb172
Pb206
Pb207
Pb208
Pb210
Bi209
Bi210
Po210
Rn220
Rn222
Ra223
Ra224
Ra226
Ra228
Ac227
Th227
Th228
Th229
Th230
Th232
Th233
Th234
Pa231
Pa232
Pa233
Pa234
Pa234M
U230
U231
U232
U233
U234
U235
U236
U237
U238
U239
U240
Np235
Np236
Np237
Np238
Np239
Pu236
Pu237
Pu238
Pu239
Pu240
Pu241
Pu242
Pu243
Pu244
Pu245
Am241
Am242
Am242M
Am243
Am244
Am244M
Cm242
Cm243
Cm244
Cm245
Cm246
Cm247
Cm248
Cm249
Cm250
Bk249
Bk250
Cf249
Cf250
Cf251
Cf252
Cf253
Cf254
Es253
Es254
Es254M

In [53]:

    

tm = Transmuter()
t_decay = 5.*24.*60.*60. # 5 days in seconds 
obs = tm.transmute(mat, t=t_decay, phi=0.0, tol=1e-7)
print obs


    

    




---------------------------------------------------------------------------
RuntimeError                              Traceback (most recent call last)
<ipython-input-53-fc7cb2f669d9> in <module>()
      1 tm = Transmuter()
      2 t_decay = 5.*24.*60.*60. # 5 days in seconds
----> 3 obs = tm.transmute(mat, t=t_decay, phi=0.0, tol=1e-7)
      4 print obs

/Users/khuff/Library/Python/2.7/lib/python/site-packages/pyne/transmute/chainsolve.pyc in transmute(self, x, t, phi, tol, log)
    127             # Find output for root of unit density and scale all output by
    128             # actual nuclide density and add to final output.
--> 129             partial = self._transmute_partial(nuc)
    130             for part_nuc, part_adens in partial.items():
    131                 y_atoms[part_nuc] = part_adens * adens + y_atoms.get(part_nuc, 0.0)

/Users/khuff/Library/Python/2.7/lib/python/site-packages/pyne/transmute/chainsolve.pyc in _transmute_partial(self, nuc)
    161         rootval = np.exp(-dest * self.t)
    162         partial = {nuc: rootval}
--> 163         self._traversal(nuc, A, partial)
    164         return partial
    165 

/Users/khuff/Library/Python/2.7/lib/python/site-packages/pyne/transmute/chainsolve.pyc in _traversal(self, nuc, A, out, depth)
    269         for child in prod:
    270             # Grow matrix
--> 271             d = self._get_destruction(child)
    272             B = self._grow_matrix(A, prod[child], d)
    273             # Create initial density vector

/Users/khuff/Library/Python/2.7/lib/python/site-packages/pyne/transmute/chainsolve.pyc in _get_destruction(self, nuc, decay)
    182         """
    183         xs_cache = self.xs_cache
--> 184         sig_a = sigma_a(nuc, xs_cache=xs_cache)
    185         d = utils.from_barns(sig_a[0], 'cm2') * xs_cache['phi_g'][0]
    186         if decay:

/Users/khuff/Library/Python/2.7/lib/python/site-packages/pyne/xs/channels.pyc in sigma_a(nuc, temp, group_struct, phi_g, xs_cache)
    374     nuc = nucname.id(nuc)
    375     key = (nuc, rxname.id('absorption'), temp)
--> 376     return xs_cache[key]
    377 
    378 

/Users/khuff/Library/Python/2.7/lib/python/site-packages/pyne/xs/cache.pyc in __getitem__(self, key)
    100                 kw['dst_phi_g'] = self._cache['phi_g']
    101                 for ds in self.data_sources:
--> 102                     xsdata = ds.discretize(**kw)
    103                     if xsdata is not None:
    104                         self._cache[key] = xsdata

/Users/khuff/Library/Python/2.7/lib/python/site-packages/pyne/xs/data_source.pyc in discretize(self, nuc, rx, temp, src_phi_g, dst_phi_g)
    187         """
    188         src_phi_g = self.src_phi_g if src_phi_g is None else np.asarray(src_phi_g)
--> 189         src_sigma = self.reaction(nuc, rx, temp)
    190         dst_sigma = None if src_sigma is None else group_collapse(src_sigma,
    191                                                         src_phi_g, dst_phi_g,

/Users/khuff/Library/Python/2.7/lib/python/site-packages/pyne/xs/data_source.pyc in reaction(self, nuc, rx, temp)
    153 
    154         """
--> 155         nuc = nucname.id(nuc)
    156         rx = rxname.id(rx)
    157         rxkey = (nuc, rx, temp) if self._USES_TEMP else (nuc, rx)

/Users/khuff/Library/Python/2.7/lib/python/site-packages/pyne/nucname.so in pyne.nucname.id (/Users/khuff/repos/pyne/build/pyne/nucname.cxx:2701)()

RuntimeError: Indeterminate nuclide form: 100900000

In [ ]:

    

vecmass = obs.mult_by_mass()
vecmol={}
for key, val in vecmass.iteritems():
    vecmol[key] = val/data.atomic_mass(key)

activity_Ci={}
for iso, mols in vecmol.iteritems():
    dec_per_sec = mols*constants.N_A*data.decay_const(nucname.id(iso))
    activity_Ci[nucname.name(iso)]= dec_per_sec/ci_dec_per_sec
    

sorted_a = sorted(activity_Ci.iteritems(), key=operator.itemgetter(0))

elem_act_Ci={}
iso_act_Ci={}
for group, elems in groups.iteritems() :
    #print group
    elem_act_Ci={}
    iso_act_Ci={}
    for num in elems : 
        elem = nucname.name(num)
        elem_act_Ci[elem] = sum([value for key, value in activity_Ci.items() if num == nucname.znum(key)])
        for key, value in activity_Ci.items():
            if num== nucname.znum(key):
                iso_act_Ci[key] = value
    sorted_a = sorted(elem_act_Ci.iteritems(), key=operator.itemgetter(1),reverse=True)
    for entry in sorted_a :  
        elem, val = entry
        #if elem_act_Ci[elem] > 0:
            #if len(elem) == 1:
                #print elem + "      " + '%.2e' % elem_act_Ci[elem]
            #else:
                #print elem + "     " + '%.2e' % elem_act_Ci[elem]
    sorted_iso_a = sorted(iso_act_Ci.iteritems(), key=operator.itemgetter(1),reverse=True)
    for entry in sorted_iso_a :  
        iso, val = entry
        if iso_act_Ci[iso] > 0:
            nspaces=10-len(nucname.name(iso))
            print nucname.name(iso) + ","+" " + '%.2e' % iso_act_Ci[iso]


    

In [ ]: