In this paper for simulating spatially and temporally varied moving storm in a watershed a distributed model was developed. The model is conducted by two major flow simulations which overland flow simulation and channel network flow simulation. Two dimensional continuity equation and momentum equation of kinematic approximation are used in the overland flow simulation. On the other hand, in the channel networks simulation two types of governing equations which are one dimensional continuity and momentum equations between two adjacent sections in a channel, and continuity and energy equations at a channel junction are applied. The finite element formulations were used in the overland flow simulation and the implicit finite difference formulations were used in the channel network simulation. The finite element formulations for the overland flow are analyzed by the Gauss elimination method and the finite difference formulations for the channel network flow are analyzed by the double sweep method having advantages of computational speed and reduced computer storages. Several recurrent coefficient equations for channel network simulation are suggested in the paper.