Computing monodromy via continuation methods on random Riemann surfaces