The paper presents a new variant of a nonsmooth continuation algorithm by means of which one can follow branches of solutions to 2D contact problems with Coulomb friction which are parameterized by the coefficient of friction. The algorithm is based on the predictor-corrector technique and uses the active set strategy implementation of the semismooth Newton method.