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)