Computation of flow velocity in gravel bed channels