AM_mcmc_refit function

Performs a Gibbs sampling reusing previous configuration