Numerical solution of the time-dependent Ginzburg-Landau equation for a superconducting mesoscopic disk: Link variable method