nugrid_set.py 文件源码

def final_bottom_envelope_set1(self):
        '''
            For paper1 marco routine:
             Numbers of remnant mass shell masses, exists also in mesa_set!
        '''
        inim=[]
        remnm=[]
        for i in range(len(self.runs_H5_surf)):
            m1p65_last=se(self.runs_H5_out[i])  
            mass_dummy=m1p65_last.se.get(m1p65_last.se.cycles[len(m1p65_last.se.cycles)-1],'mass')
            top_of_envelope=mass_dummy[len(mass_dummy)-1]
            h_dummy=m1p65_last.se.get(m1p65_last.se.cycles[len(m1p65_last.se.cycles)-1],'iso_massf','H-1')

            for j in range(len(mass_dummy)):
                    if h_dummy[j] > 0.05:
                            bottom_of_envelope = mass_dummy[j]
                            break
            inim.append(m1p65_last.get("mini"))
            remnm.append(bottom_of_envelope)
        print "M_initial | M_remn/bottom of envelope"
        for i in range(len(inim)):
            print inim[i],"|",remnm[i]