docs: def set_plot_CC_T_rho_new(self,fig='CC evol',linestyle=['-'],burn_limit=0.997,color=['r'],marker=['o'],markevery=500): 
  		'''
  			Plots 
  			end_model - array, control how far in models a run is plottet, if -1 till end
  			symbs_1  - set symbols of runs
  		'''
  		if len(linestyle)==0:
  			linestyle=200*['-']
  		plt.figure(fig)
                  for i in range(len(self.runs_H5_surf)):
                          sefiles=se(self.runs_H5_out[i])
  			t1_model=-1
                          mini=sefiles.get('mini')
  			zini=sefiles.get('zini')
  			label=str(mini)+'$M_{\odot}$, Z='+str(zini)
                          model=sefiles.se.cycles
                          model_list=[]
                          for k in range(0,len(model),1):
                                  model_list.append(model[k])
  			print 'REad Rho,T, this might take a while...'
                          rho1=sefiles.get(model_list,'rho')   #[:(t1_model-t0_model)]
                          T1=sefiles.get(model_list,'temperature')#[:(t1_model-t0_model)]
  			print 'finished'
                          rho=[]
                          T=[]
  			T_unit=sefiles.get('temperature_unit')
  			labeldone=False
  			for k in range(len(model_list)):
  				#print 'test model ',model_list[k]
                                  t9max=max(np.array(T1[k])*T_unit/1.e9)
                                  #T.append(max(t9))
  				rho1max=max(rho1[k])
  				print 'model',model_list[k]
  				print 'maxT, maxrho'
  				print t9max,rho1max
  				if k==0:
  					t9_prev=t9max
  					rho1_prev=rho1max
  					idx_T=0
  					idx_rho=0
  					continue
  				if t9max>t9_prev:
  					idx_T=k
  					t9_prev=t9max
  				if rho1max>rho1_prev:
  					idx_rho=k
  					rho1_prev=rho1max
  			print 'found highest rho',idx_rho,max(np.array(T1[idx_rho])*T_unit/1.0e9),max(rho1[idx_rho]),model_list[idx_rho]
  			print 'found highest T',idx_T,max(np.array(T1[idx_T])*T_unit/1.0e9),max(rho1[idx_T]),model_list[idx_T]
  			if idx_T==idx_rho:
  				x=np.array(T1[idx_T])*T_unit/1e9
  				y=rho1[idx_T]
  				rho1=[]	
  				T1=[]
  				for k in range(len(x)):
  				    if not y[k]==1.0:
  			        	rho1.append(y[k])
          				T1.append(x[k])
  				x=T1
  				y=rho1
  				plt.plot(x,y,label=label,color=color[i],marker=marker[i],linestyle=linestyle[i],markevery=markevery)
                                  #rhoi.append(max(rho1[k]))	       
  			else:
  				#for max T
  				x=np.array(T1[idx_T])*T_unit/1e9
  				y=rho1[idx_T]
                                  rho_temp=[]
                                  T_temp=[]
                                  for k in range(len(x)):
                                      if not y[k]==1.0:
                                          rho_temp.append(y[k])
                                          T_temp.append(x[k])
                                  x=T_temp
                                  y=rho_temp
  				plt.plot(x,y,label=label,color=color[i],marker=marker[i],linestyle=linestyle[i],markevery=markevery)
  				#for max rho
          			x=np.array(T1[idx_rho])*T_unit/1e9
  				y=rho1[idx_rho]
                                  rho_temp=[]
                                  T_temp=[]
                                  for k in range(len(x)):
                                      if not y[k]==1.0:
                                          rho_temp.append(y[k])
                                          T_temp.append(x[k])
                                  x=T_temp
                                  y=rho_temp
  				plt.plot(x,y,color=color[i],marker=marker[i],linestyle=linestyle[i],markevery=markevery)