Free surface models of partially molten rock with visco-elasto–plastic rheology: numerical solutions using a staggered-grid finite-difference discretisation