To develop the numerical model, we solve for the melt fraction and the velocity of each phase as a function of depth within the compacting medium. This is accomplished by solving statements of conservation of mass and momentum through finite volumes methods. We use an expression for the effective matrix viscosity that is present in the conservation of momentum statements derived for repacking. In this case the effective matrix viscosity is the sum of viscous contributions from hydrodynamic and particle-particle (frictional) interactions (Boyer et al. , 2011). Finally, we perform Markov-chain Monte-Carlo (MCMC) inversions to explore the best fitting solutions and the parameters for the effective matrix viscosity that result in the most likely fit to these sets of experiments.