def find_first_TP(self):
star_mass         = self.get('<STR_LIT>')<EOL>he_lumi           = self.get('<STR_LIT>')<EOL>h_lumi            = self.get('<STR_LIT>')<EOL>mx2_bot           = self.get('<STR_LIT>')*star_mass<EOL>try:<EOL><INDENT>h1_boundary_mass  = self.get('<STR_LIT>')<EOL>he4_boundary_mass = self.get('<STR_LIT>')<EOL><DEDENT>except:<EOL><INDENT>try:<EOL><INDENT>h1_boundary_mass  = self.get('<STR_LIT>')<EOL>he4_boundary_mass = self.get('<STR_LIT>')<EOL><DEDENT>except:<EOL><INDENT>pass<EOL><DEDENT><DEDENT>TP_bot=np.array(self.get('<STR_LIT>'))*np.array(self.get('<STR_LIT>'))<EOL>TP_top=np.array(self.get('<STR_LIT>'))*np.array(self.get('<STR_LIT>'))<EOL>lum_array=[]<EOL>activate=False<EOL>models=[]<EOL>pdcz_size=[]<EOL>for i in range(len(h1_boundary_mass)):<EOL><INDENT>if (h1_boundary_mass[i]-he4_boundary_mass[i] <<NUM_LIT>) and (he4_boundary_mass[i]><NUM_LIT>):<EOL><INDENT>if (mx2_bot[i]>he4_boundary_mass[i]) and (he_lumi[i]>h_lumi[i]):<EOL><INDENT>if TP_top[i]>he4_boundary_mass[i]:<EOL><INDENT>pdcz_size.append(TP_top[i]-TP_bot[i])<EOL>activate=True<EOL>lum_array.append(he_lumi[i])<EOL>models.append(i)<EOL><DEDENT><DEDENT>if (activate == True) and (he_lumi[i]<h_lumi[i]):<EOL><INDENT>if max(pdcz_size)<<NUM_LIT>:<EOL><INDENT>active=False<EOL>lum_array=[]<EOL>models=[]<EOL>print('<STR_LIT>')<EOL><DEDENT>else:<EOL><INDENT>break<EOL><DEDENT><DEDENT><DEDENT><DEDENT>t0_model = models[np.argmax(lum_array)]<EOL>return t0_model<EOL>
Find first TP of the TPAGB phase and returns the model
number at its LHe maximum.

Parameters
----------