docs: def set_plot_CC_T_rho(self,linestyle=[],burn_limit=0.997,color=['r'],marker=['o'],nolabelZ=False,markevery=500): 
  		'''
  			Plots HRDs
  			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('CC evol')
                  for i in range(len(self.runs_H5_surf)):
                          sefiles=se(self.runs_H5_out[i])
  			t1_model=-1
  			sefiles.get('temperature')
  			sefiles.get('density')
                          mini=sefiles.get('mini')
  			zini=sefiles.get('zini')
                          model=sefiles.se.cycles
                          model_list=[]
                          for k in range(0,len(model),1):
                                 model_list.append(model[k])
                          rho1=sefiles.get(model_list,'rho')   #[:(t1_model-t0_model)]
                          T1=sefiles.get(model_list,'temperature')#[:(t1_model-t0_model)]
                          rho=[]
                          T=[]
  			T_unit=sefiles.get('temperature_unit')
  			labeldone=False
                          for k in range(len(model_list)):
  				t9=np.array(T1[k])*T_unit/1e9
  				T_delrem=[]
  				rho_delrem=[]
  				if k ==0:
  					T.append(max(t9))
  					rho.append(max(rho1[k]))
  					for h in range(len(rho1[k])):
  						if rho1[k][h] < 1e3:
  							T_delrem.append(t9[h])
  							rho_delrem.append(rho1[k][h])
  					if nolabelZ==True:
  						plt.plot(T_delrem,rho_delrem,label=self.extra_label[i].split(',')[0],color=color[i],marker=marker[i],markevery=markevery)
  					else:
  						plt.plot(T_delrem,rho_delrem,label=self.extra_label[i],color=color[i],marker=marker[i],markevery=markevery)
  				else:
  					if (max(rho)<max(rho1[k]) or max(T)<max(t9)):
                                          	for h in range(len(rho1[k])):
  							if rho1[k][h] > 1e3:
                                                  		T_delrem.append(t9[h])
                                          	      		rho_delrem.append(rho1[k][h])
  						if labeldone==True:	
  							plt.plot(T_delrem,rho_delrem,color=color[i],marker=marker[i],markevery=markevery)
  						else:
  							label=str(mini)+'$M_{\odot}$, Z='+str(zini)
  							if nolabelZ==True:
  								plt.plot(T_delrem,rho_delrem,label=label.split(',')[0],color=color[i],marker=marker[i],markevery=markevery)
  							else:
  								plt.plot(T_delrem,rho_delrem,label=label,color=color[i],marker=marker[i],markevery=markevery)
  							labeldone=True
  						T.append(max(t9))
  						rho.append(max(rho1[k]))
  					else:
  						break
  					#else:
  						
  		plt.xlabel('$T_9 [GK]$',size=22)
  		plt.ylabel(r'$\rho [cm^{-3}]$',size=22)
  		plt.yscale('log')
  		plt.xscale('log')
  		plt.legend(loc=2)