Iteratively solves for critical depth in a circular culvert using the critical flow condition Q^2/g = A^3/T. Also computes critical velocity and Froude number.