nugrid_set.py 文件源码

def read_yield_sn1a_tables(self,sn1a_table,isotopes):
        f1=open(sn1a_table)
        lines=f1.readlines()
        f1.close()
        iso_1a=[]
        yield_1a=[]
        for line in lines:
            #for header
            if '#' in line:
                continue
            iso_1a.append(line.split()[0])
            yield_1a.append(float(line.split()[1]))             


        yields=[]
        #fill up the missing isotope yields with zero
        for iso in isotopes:
            if iso in iso_1a:
                idx=iso_1a.index(iso)
                yields.append(yield_1a[idx])
            else:
                yields.append(0.)       
        return yields