docs: def final_bottom_envelope_set1(self):
  		'''
  			For paper1 marco routine:
  			 Numbers of remnant mass shell masses, exists also in mesa_set!
  		'''
  		inim=[]
  		remnm=[]
  		for i in range(len(self.runs_H5_surf)):
  			m1p65_last=se(self.runs_H5_out[i])	
  			mass_dummy=m1p65_last.se.get(m1p65_last.se.cycles[len(m1p65_last.se.cycles)-1],'mass')
  			top_of_envelope=mass_dummy[len(mass_dummy)-1]
  			h_dummy=m1p65_last.se.get(m1p65_last.se.cycles[len(m1p65_last.se.cycles)-1],'iso_massf','H-1')
  
  			for j in range(len(mass_dummy)):
          			if h_dummy[j] > 0.05:
                  			bottom_of_envelope = mass_dummy[j]
                  			break
  			inim.append(m1p65_last.get("mini"))
  			remnm.append(bottom_of_envelope)
  		print "M_initial | M_remn/bottom of envelope"
  		for i in range(len(inim)):
  			print inim[i],"|",remnm[i]