The time evolution of power systems is modeled by a system of differential and algebraic equations. The variables involved in the system may exhibit different time scales. In standard numerical time integration methods the most active variables impose the time step for the whole system. We present a strategy, which allows the use of different, local time steps over the variables. The partitioning of the components of the system in different classes of activity is performed automatically and is based on the topology of the power system.