The problem of finding the optimal control in numerical computer simulations of quantum control phenomena is usually addressed through the introduction of monotonically convergent algorithms that are guaranteed to improve the cost functional at each step. A recent extension of these algorithms implements a search for a control with given bounds. Within this context, this paper will present a generalisation of this algorithm to a more general class of local in time constraints. Numerical results that illustrate some of the properties of the algorithms introduced are given.