The results shown in the following were calculated using the drag coefficient due to the formula of Matthias and Godfrey
[16], bottom roughness parameter as 0.01 m suggested by Weatherly and Martin [2] and bottom coefficient from calculation
based on the logarithm law of the wall suggested by Mellor [3].