A one-dimensional model capable of simulating flood wave propagation in a river or network of channels is presented. The computer model is programmed to provide maximum flexibility in the adaptation of channel geometry, the specification of conveyance properties, and the treatment of boundary conditions. An equation transformation procedure is employed in the model to minimize computer storage and execution time requirements by reducing the order of the resultant coefficient matrices. Based on a four-point implicit finite-difference approximation of the governing, nonlinear, flow equations, the model can be used to simulate the wide range of flow conditions typically encountered in various natural waterbody systems. Two particular applications are presented to demonstrate the computational features and capabilities of the model in the simulation of flood wave propagation.