nugrid_set.py 文件源码

def write_gce_input_lifetimes(self,file_name="isotopic_table.txt",final_cycles=[]):

        '''
            Calculates and writes lifetime in file_name
            X_carbon_center >  0.4

        '''

        ###Get lifetimes###

        print '##############GETTING THE LIFETIMES##############'
        time=[]
        model=[]
        mz=[]
        for i in range(len(self.runs_H5_surf)):
            sefiles=se(self.runs_H5_surf[i])
            mass=sefiles.get("mini")
            metallicity=sefiles.get("zini")
            mz1='(M='+str(round(mass,2))+',Z='+str(metallicity)+')'
            cycs=[]
            if final_cycles[i]>-1:
                final_cycle=final_cycles[i]
            else:
                final_cycle=int(sefiles.se.cycles[-1])
            time1=(sefiles.get(final_cycle,'age')*sefiles.get('age_unit'))/31557600.
            time.append(time1)
            mz.append(mz1)

        f1=open(file_name,'r')
        lines=f1.readlines()
        f1.close()
        i=-1
        line1=''
        print 'Adding for ',mz,'lifetime : ',time
        while (True):
            i+=1
            if i>len(lines)-1:
                break
            line=lines[i]   
            line1+=lines[i]
            for k in range(len(mz)):
                if mz[k] in lines[i]:
                    line1+=('H Lifetime: '+'{:.3E}'.format(time[k])+'\n')
                    break           

        f1=open(file_name,'w')
                f1.write(line1)
                f1.close()