We give a combinatorial algorithm for equivariant embedded resolution of singularities of a toric variety defined over a perfect field. The algorithm is realized by a finite succession of blowings-up with smooth invariant centres that satisfy the normal flatness criterion of Hironaka. The results extend to more general varieties defined locally by binomial equations. Comment: 45 pages