nugrid_set.py 文件源码

def initial_finall_mass_relation(self,marker='o',linestyle='--'):
        '''
            INtiial to final mass relation
        '''

        final_m=[]
                ini_m=[]
                for i in range(len(self.runs_H5_surf)):
            sefiles=se(self.runs_H5_out[i])
                        ini_m.append(sefiles.get("mini")) 
            h1=sefiles.get(int(sefiles.se.cycles[-2]),'H-1')
            mass=sefiles.get(int(sefiles.se.cycles[-2]),'mass')
            idx=-1
            for k in range(len(h1)):
                if h1[k]>0.1:
                    idx=k
                    break
            final_m.append(mass[idx])
        label='Z='+str(sefiles.get('zini'))
        plt.plot(ini_m,final_m,label=label,marker=marker,linestyle=linestyle)
        plt.xlabel('$M_{Initial} [M_{\odot}]$',size=23)
        plt.ylabel('$M_{Final} [M_{\odot}]$',size=23)