def multi_DUP(self,t0_model=[],h_core_mass=False,plot_fig=True,linestyle=[],marker=[],color=[]):
'''
Plot (!) DUP parameter lambda versus mass.
There is a algorithm finding the first TP (and set t0_model) for reach run dir
which worked fine for set1. But there could be a problem with SAGB stars
so use with caution! You can use a manual input (the array t0_model) to
skip this step.
It also returns peak lum and corresponding
model number of each TP.
h_core_mass - If True: plot dependence from h free core , else star mass
t0_model - Array of first TPs of models.
examples of set1:
t0_model paramter:
z1e-2:1.65-5: [13690,3120,3163,5306]
z2e-2:1.65-5: [16033,6214,3388,5368]
dirs=["M1.65Z2.0e-02/LOGS","M2.00Z2.0e-02/LOGS","M3.00Z2.0e-02/LOGS","M5.00Z2.0e-02/LOGS"]
dirs=["M1.65Z1.0e-02/LOGS","M2.00Z1.0e-02/LOGS","M3.00Z1.0e-02/LOGS","M5.00Z1.0e-02/LOGS"]
note: M3_1e-2 shows decrease in core mass in the end
set.multi_DUP(t0_model=[13690,3120,3163,5306])
set.multi_DUP(t0_model=[16033,6214,3388,5368])
'''
if len(t0_model)==0:
t0_model=self.set_find_first_TP()
dirs=self.run_LOGS
historydata=self.run_historydata
marker_type=marker
#line_style=['--','-','-.',':']
peak_lum_model_array=[]
h1_mass_min_DUP_model_array=[]
for i in range(len(dirs)):
#color,marker_type)
peak_lum_model,h1_mass_min_DUP_model = historydata[i].find_TP_attributes(0,t0_model[i],color[i],marker_type[i],h_core_mass,no_fig=plot_fig)
peak_lum_model_array.append(peak_lum_model)
h1_mass_min_DUP_model.append(h1_mass_min_DUP_model)
return peak_lum_model_array,h1_mass_min_DUP_model_array
###the following methods allow are parts of Falks vis3.py file, thanks Falk
评论列表
文章目录
