We present a numerical and experimental study of the generation of harmonic mode locking in a silica toroid microcavity. We use a generalized mean-field Lugiato-Lefever equation and solve it with the split-step Fourier method. We found that a stable harmonic mode-locking regime can be accessed when we reduce the input power after strong pumping even if we do not carefully adjust the wavelength detuning. This is due to the bistable nature of the nonlinear cavity system. The experiment agrees well with the numerical analysis, where we obtain a low-noise Kerr comb spectrum with a narrow longitudinal mode spacing by gradually reducing the input pump power after strong pumping. This finding clarifies the procedure for generating harmonic mode locking in such high-Q microcavity systems.