Computation of critical and normal depth in trapezoidal, triangular, or rectangular channels
(c) 2001-2003 V. M. Ponce, Version 3.4, October 11, 2003