Add discharge function to MHKiT

[jmcvey3]

Also capable of estimating fluid power and Reynolds number. The function natively fills in velocity between the riverbed and the lowest depth bin by assuming velocity at the seafloor is 0 and interpolating up to the lowest depth bin. Velocity at the rivers surface is assumed to be the same velocity as the uppermost depth bin.

Example notebook in progress.