	el_found, abundances = abunds_func(x_0)
	blobs_dtype = [('ep_r','f8'), ('rew_r','f8')]

	for el in el_found:
	    blobs_dtype.append((el, 'f8'))
	    blobs_dtype.append((el+'_sigma_mean', 'f8'))


	sampler = emcee.EnsembleSampler(n_walkers, n_dim, obj_func, blobs_dtype=blobs_dtype, args=(len(el_found),sun_el,sun_abs))
	results = sampler.run_mcmc(x_ball, n_steps, progress=True)

	flat_blob = sampler.get_blobs().reshape((n_walkers*n_steps))